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ABSTRACT 
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MODAL FLUTTER ANALYSIS OF ADVANCED TURBOPROPELLERS 


1 .1 Introduction 

An existing capability to conduct modal flutter analysis of tuned bladed- 
shrouded discs in NASTRAN (Ref. 1) has been modified to analyze the subsonic 
unstalled flutter characteristics of advanced turbopropellers in NASTRAN 
Level 17.7. 

The modifications pertain to the inclusion of oscillatory modal aerodynamic 
loads of blades with large (backward and forward) varying sweep. 

The following section summarizes the theoretical aspects of turboprop 
flutter analysis from Ref. 2. 

1 .2 Theory 

Multi-bladed advanced turbopropellers are geometrically cyclic structures 
with thin blades of low aspect ratio and varying sweep. The blades are mounted 
on a relatively rigid hub and, therefore, can be considered to be structurally 
independent. This permits modal analysis of only one root-fixed blade without 
recourse to special harmonic analysis techniques applicable to cyclic structures. 
From a flutter aerodynamics viewpoint, the estimation of the generalized oscil- 
latory aerodynamic loads on the propeller blades depends on the aerodynamic 
theory employed. In the present capability, the two-dimensional subsonic cascade 
unsteady aerodynamic theory of Jones and Rao (Ref. 3) is applied in a strip 
theory manner similar to that of Barmby et al (Ref. 4) with appropriate modifica- 
tions recognizing the variability of the blade sweep and chord with radius. 

To facilitate the use of the two-dimensional cascade theory, the aerodynamic 
model of the blade is based on a grid defined by the intersection of a series of 
chords and "computing stations" as shown by the thick solid lines in Figure 1. 

The chords are selected normal to any spanwise reference curve such as the blade 
leading edge. Due to its resemblance to the structural model of the blade, and the 


adequacy of a relatively coarse grid to describe the spa.^wise flow variations, 
the aerodynamic model is chosen as a subset of the structural model. 

The modified two-dimensional cascade theory is applied on each of these 
chords to determine the generalized aerodynamic forces acting on the associated 
strips. The strip results are added to obtain the blade aerodynamic matrix. 

An overall flowchart for modal flutter analysis of advanced turbopropellers 
is shown in Figure 2. 

The User's, Programmer's and Demonstration Manuals are presented in Sections 
2, 3 and 4, respectively . 
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Figure 2 . Overall Flowchart of Advanced Turbopropeller Modal Flutter 
Analysis (continued). 
























MODAL FLUTTER ANALYSIS OF ADVANCED TURBOPROPELLERS 


2.1 Introduction 

Subsonic unstalled (modal) flutter analysis of advanced turbopropellers 
can be conducted using this capability. The AERO APPROACH Rigid Format 9, 

Series R, in NASTRAN Level 17.7 (Section 2.6 ) with modified functional modules 
APDB, AMG and AMP forms the basis for modal flutter analysis. Section 4 
demonstrates the use of the capability. Section 1 presents the theoretical 
aspects. 

2.2 Solution Phases 

As illustrated by the example in the Demonstration Manual (Section 4), the 
flutter analysis can, in general, be conducted in three phases--phase 1 to 
generate differential stiffness, phase 2 to compute natural modes and frequencies 
and phase 3 to compute flutter eigenvalues. 

Phase 1 uses DISP RF4 while phases 2 and 3 use AERO RF9. 

2.3 NASTRAN Model 

With the assumption of tuned blades mounted on a relatively rigid hub, the 
user models one blade of the advanced turbopropel ler as shown in Figure 1. 

The structural model is prepared using the general modelling capabilities 
of NASTRhN. The basic coordinate system is fixed to the rotating propeller such 
that the X-axis always coincides with the axis of rotation and is directly 
opposing the flight direction. Location of the origin is arbitrary. 

The XZ plane is approximately located such as to contain the maximum pro- 
jected area of the blade being modelled. This orientation is consistent with 
the internally generated chordl ine coordinate systems for the unsteady aerodynamics. 

The aerodynamic model comprises a grid defined by the intersection of a 
series of chords and "computing stations" (Figure 1). The chords are selected 
normal to any spanwise reference curve such as the blade leading edge. The 


2.1 


choice of the numberand location of the chords and the competing stations is 
dictated by the expected variation of the relative flow properties across the 
blade span, and the complexity of the mode shapes exhibited by the propeller 
blade. Due to its resemblance to the structural model of the blade, and the 
adequacy of a relatively coarse grid to describe the spanwise flow variations, 
the aerodynamic model is chosen as a subset of the structural model as shown in 
Figure 1 . 

The aerodynamic grid is specified on STREAML1 bulk data cards. 

2.4 Remarks on the Use of Some Cards 

As noted in the Theoretical Manual (Section 1), the present capability to 
analyze turboprop flutter is derived by modifying the modal flutter analysis 
capability originally developed for bladed discs. These modifications as 
reflected in the user input preparation are discussed below. 

The NASTRAN card is required for flutter analysis only if sweep aerodynamic 
effects are to be included. It is placed preceding the Executive Control Deck, 
and is specified as 

NASTRAN SYSTEM (76) = 1 

Section 2.1 of the NASTRAN User's Manual, and Section 7.3.1 of the NASTRAN 
Programmer's Manual discuss this card in detail. 

In light of the assumption of structurally independent blades, 

a) CYJOIN bulk data cards are required merely for their presence in the 
Bulk Data Deck, 

b) PARAM KINDEX is set to zero to save computational time in the real 
eigenvalue extraction process and 

c) PARAM MTYPE is set to COSINE (default) because of KINDEX = 0. 

The STREAML2 bulk data card (see Section 2.5) has been modified to include 
the parameters associated with the swept blade aerodynamics. Figure 2 defines 
some of these parameters. In this figure, A_B_, AB and A + B + represent three 
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successive chords with points A's on the leading edge. For the chord AB, at 
any operating condition WA represents the absolute inflow velocity while AU 

(= x RA) is the blade (tangential) velocity. WA and AU uniquely define a 

plane in which the inflow properties are defined. 

In the plane WAU, VA - WA - AU represents the relative inflow velocity. 

C/t represents the chordwise, cascade relative inflow velocity (Field 2, con- 
tinuation card). Mach number in Field 8 is based on CA. AI is the line of 
intersection between the axial plane through point A and the plane WAU. Angle 
IAV defines the relative inflow angle 8 (shown positive). 

The angle of sweep A is defined as the angle of inclination of the chord 

BA with the plane WAU. A shown in Figure 2 is positive. 

AD is the projection of AC (BA extended to C) in the plane WAU. Angle IAD 

represents the stagger angle A, and is shown positive. 

A local coordinate system xyz is internally defined at the leading edge 
point A of the chord AB such that x is directed along AB. y is defined normal 
to the 'mean' surface containing the points A_, A, A + , B + , B and B_. The unit 

vector along y, for the sense of Q shown in Figure 2, is given by 

(A J} + ) x (AB) (AB) x (\B_) 

\(rf + ) x (AB)1 | (AB) x (A^f )| 

Modal translations along y and rotations about x are used in deriving the 
generalized airforce matrix. For the opposite sense of rotation, xyz is internally 
defined to be left handed with y reversing direction. The shaded area about the 
chord AB represents the strip of integration associated with AB. 

The FLFACJ, FLUTTER and MKAEROi bulk data cards (Section 2.5) have beers 
modified to specify the interblade phase angles. Referring to the sketch below, 
a positive interblade phase angle implies that blade 1 of the two-dimensional 
cascade leads the reference blade 0. 
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Blade 1 

Blade 0 (ref. ) 


Section 2.5 contains all the bulk data cards associated with aerodynamics 
data specification for the turboprop flutter analysis. Section 2.6 describes 
all the PARAMeters. 
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2.5 BULK DATA DECK 
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Input Data Card AER0 Aerodynamic Physical Data 

Description : Gives basic aerodynamic parameters. 


Format and Examples - 


1 23456789 10 


AER0 

ACSID 

VELOCITY 

REFC 

RH0REF 

SYMXZ 

SYMXY 




AERO 

3 

1.3+4 

100. 

1.-5 


1 





Field 


Contents 


ACSID 

VELOCITY 

REFC 

RH0REF 

SYMXZ 

SYMXY 


Aerodynamic coordinate system Identification (Integer 1 0). See Remark 2. 
Velocity (Real). 

Reference length (for reduced frequency) (Real). 

Reference density (Real). 

Syimetry key for aero coordinate x-z plane (Integer) (+1 for sym, »0 for no sym, 
-1 for anti -sym). 

Symmetry key for aero coordinate x-y olane can be used to simulate ground 
effects (Integer), same code as SYMXZ. 


Remarks : 1. This card Is required for aerodynamic respoise problems. Only one AER0 card is 

allowed. 

2. The ACSID must be a rectangular coordinate system. Flow Is In the positive x 
dl recti on. 


1. i Ze$«se~c*. eeFc/2, 

C k- v) 


2.7 



OttlGttiftL PAGE m 

OF POOR QUALITY 


BULK DATA DECK 


Input Data Card FLFACT Aerodynamic Physical Data 


IL e i . cr j - PAlPf 1 : Used to specify densities, Mach numbers or Interblade phase angles 

and reduced frequencies for flutter analysis. 


Format and Example: 


1 , l , ? , 4 , 5 6 , Z , 8 , 9 , lH 


FLFACT 

SID 

FI 

F2 

F3 

F4 

F5 

F6 

F7 

ABC 

FLFACT 

97 

.3 

.7 

3.5 





abc 

+BC 

F8 

F9 

--etc.— 


















Alternate Form : 


FLFACT 

SID 

FI 

THRU 

FNF 

NF 

FMID 

XI 

x 


FLFACT 

201 

.200 

THRU 

.100 

11 

.133333 





Field 


C ontents 


SID Set identification number {Unique Integer > 0). 

Fi Aerodynamic factor (Real). 


Remarks : 1. These factors must be selected by a FLUTTER data card to be used by NASTRAN. 

2. Imbedded blank fields are forbidden. 

3. Parameters must be listed in the order in which they are to be used within the 
looping of flutter analysis. 

4. For the alternate form, NF must be greater than 1. F . . must lie between F. and 

F nf , otherwise F mi(j will be set to (F 1 + F Np )/2. Then 0 1 


_ V F NF-W<" M > + 

1 ' ^NF-midH^-D^F^yd-l) 


i = 1,2 NF 


The use of F m - d (middle factor selection) allows unequal spacing of the factors. 
F mid = 2fr ] F NF^ F l +F NF^ 9ives e ^ ual values to increments of the reciprocal of Fj. 




BULK DATA DECK 

Input Data Card FLUTTER Aerodynamic Flutter Data 


omomL uses 

0F p 00fi QUAtfTY 


Description : Defines data needed to perform flutter analysis. 


Format and Example : 


1 

2 1 

3 

4 

5 

6 

7 

6 

9 

10 

FLUTTER 

SID 

METH0D 

DENS 

MACH 

RFREQ 

IMETH 

NVALUE 

EPS 


FLUTTER 

19 

K 

119 

219 

319 

S 

5 

1.-4 



Field 


Contents 


SID 

METH0D 

DENS 

MACH 

RFREQ (or VEL) 

IMETH 

NVALUE 

EPS 


Set identification number (Unique Integer > 0). 

Flutter analysis method, "K" for K-method, "PK" for P-K method, "KE" for the 
K-method restricted for efficiency. 

Identification number of an FLFACT data card specifying density ratios to be 
used in flutter analysis (Integer > 0). 

Identification number of an FLFACT data card specifying 
MACH numbers or interblade phase angles (m) to be used in 
f.l utter analysis (integer > 0). 

Identification number of an FLFACT data card specifying reduced frequencies (k) 
to be used in flutter analysis (Integer > 0); for the p-k method, the velocity. 

Choice of interpolation method for matrix interpolation (BCD: L ° linear, 

S B surface). 

Number of eigenvalues for output and plots (Integer > 0). 

Convergence parameter for k; used in the P-K method ( Real )(default « 1 0“ 3 ) . 


Remarks : 1. The FLUTTER data card must be selected in Case Control Deck (FMETH0D = SID). 


2. The density is given by DENS ■ RH0REF, where RH0REF is the reference value given on 
the AER0 data card. 


3. The reduced frequency is given by k = (REFC*w/2> V) , where REFC is given on the AER0 
data card, u is the circular frequency and V is the velocity. 

4. An eigenvalue is accepted in the P-K method when |k - k estimate | < EPS. 
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BULK DATA DECK 
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Input Data Card MKAER01 Mach Number - Frequency Table 

Descrl pti on : Provides a table of Mach numbers or interblade phase angles (ro) 

and reduced frequencies (k) for aerodynamic matrix calculation. 

Format and Example : 


1 2 3 4 5 6 7 8 9 10 


MKAER01 

m 1 

1 

m 3 

EH 


m 6 

EME 

mg 

A8C 

MKAER01 

.1 

.7 







+ABC 

+BC 

k i 

k 2 

k 3 

k 4 

k 5 

k 6 

Jiz 

ie 


+BC 

.3 

.6 

1.0 





j 


Field 


Contents 


k 3 


List of Mach nurribers (Real; 1 < 1 < B). 

List of reduced frequencies (Real > 0.0, 1 i j iB), 


Remarks : 1 . 

2 . 

3. 

4. 


Blank fields end the list, and thus cannot be used for 0.0. 

All combinations of (m,k) will be used. 

The continuation card is required. 

Since 0.0 is not allowed, it may be simulated with a very small number 
such as 0.0001. 


Mach numbers are input for wing flutter and interblade phase angles for 
blade flutter. 
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BULK DATA DECK 


Input Data Card MKAER02 Mach Number - Frequency Table 

Oescrl ption: Provides a list of Mach numbers or interbladc phase angles (m) and 

reduced frequencies (k) for aerodynamic matrix calculation. 


Format and Example : 


123456789 10 


MKAER02 

mm 

m 

m 2 

k 2 

m 3 

k 3 

m 4 

Uia 


MKAER02 

.10 

— 

.30 

.10 

.60 

.70 

.30 

.70 

1.0 



Field 


Contents 


m. 

k i 


List of Mach numbers (Real > 0.0). 

List of reduced frequencies (Real > 0.0). 


Remarks : 1. This card will cause the aerodynamic matrices to be computed for a set of parameter 

pairs. * 

2. Several MKAER02 cards may be in the deck. 

3. Imbedded blank pairs are skipped. 

4. Mach numbers are input for wing flutter and interblade phase angle for blade 
flutter. 
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BULK DATA DECK 


Input Data Card STREAML 1 Blade Streamline Data 

Descrl ptlon : Defines grid points on the blade streamline from blade leading 

edge to blade trailing edge. 

Format and Example : 


1 23456789 10 


| STREAML 1 

SLN 

G 1 

G2 

G 3 

G4 

G5 

G6 

G7 

+ABC 

1 STREAHLl 

3 

2 

4 

6 

8 

10 








♦ABC 

G8 

G9 

-etc- 







+ABC 










Alternate Form: 

STREAML 1 

SLN 

GID1 

"THRU" 

GID2 

X 

X 

X 

X 


STREAML 1 

5 

6 

THRU 

12 







Field Contents 

SLN Streamline number (Integer > 0). 

G 1 , G1D1 Grid point Identification numbers (Integer > 0). 


Remarks : 

1. This card is required for blade steady aeroelastic and blade flutter 
problems. 

2. There must be one STREAML1 card for each streamline on the blade. 

For blade flutter problems, there must be an equal number of STREAHLl 
and STREAML2 cards. 

3. The streamline numbers, SLN, must increase with Increasing radial 
distance of the blade section from the axis of rotation. The 
lowest and the highest SLN, respectively, will be assumed to 
represent the blade sections closest to and farthest from the axis 
of rotation. 

4. All grid points should be unique. 

5. All grid points referenced by GID1 through G!D2 must exist. 

6. Each ST RE AML 1 card must have the same number of grid points. The 
nodes must be input from the blade leading edge to the blade 
trailing edge in the correct positional order. 





BULK DATA DECK 


Input Data Card STREAHL2 Blade Streamline Data 

Description : Defines aerodynamic data for a blade streamline. 

Format and Example : 


STREAML2 

SLN 

NSTNS 

STAGGER 

CHORD 

RADIUS/ 

DCBDZB 

BSPACE 

MACH 

DEN 

+abc 

STREAML2 

2 

3 

23.5 

1.85 

6.07 

.3B6 

.934 

.066 



+abc 

VEL 

FLOWA/ 

SWEEP 






— 


+A8C 

1014.2 

55.12 









Field 

SLM 

NSTNS 

STAGGER 

CHORD 

RADIUS/DCBDZB 


BSPACE 

MACH 

DEN 

VEL 

FLOWA/SWEEP 


Contents 

Streamline number {Integer >0) 

Number of computing stations on the blade streamline. 

(3 ^ NSTNS < 10, Integer) 

Blade stagger angle (-90.0 <stagger <90.0, degrees) 

Blade chord (real >0.0) 

Radius of streamline for flutter analysis without sweep effects 
(real >0.0) or^ 

3C/3Z for flutter analysis with sweep effects. C Is the swept 
chord and 1 is the (local) spanwise reference direction (real) 
Blade spacing (real >0.0) 

Relative flow mach number at blade leading edge (real >0.0) 

Gas density at blade leading edge (real >0.0) 

Relative flow velocity at blade leading edge (real >0.0) 

Relative flow angle at blade leading edge for flutter analysis 
without sweep effects (-90.0 <FL0WA <90.0 degrees) 0 £ 

Blade sweep angle for flutter analysis with sweep effects 
(-90.0 <SWEEP <90.0 degrees) 


2 




Remarks: 


1. At least three (3) and no more than fifty (50) STREAML2 cards are required 
for a blade flutter analysis. 

2. The streamline number, SLN, must be the same as Its corresponding SLN on 
a STREAML1 card. There must be a STREAML1 card far each STREAML2 card, 

3. It is not required that all streamlines be used to define the aerodynamic 
matrices used In blade flutter analysis. 

4. For flutter analysis with sweep effects, the use of the NASTRAN card Is 
required as follows: 

NASTRAN SYSTEM (76) = 1 

Refer to Section 2.1 of the User's Manual and Section 6.3.1 of the Programmer' 
Manual for description and placement in the Executive Control Deck. 
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COMPRESSOR BLADE CYCLIC MODAL FLUTTER ANALYSIS 

2.6 COMPRESSOR BLADE CYCLIC MODAL FLUTTER ANALYSIS 

2 . 6.1 DMAP Sequence For Compressor Olade Cyclic Modal Flutter Analysis 

RIGID FORMAT DMAP LISTING 

aero approach, rigid FORMAT 9 

LEVEL 2.0 NASTRAN DMAP COMPILER - SOURCE I. * STING 


OPTIONS IN EFFECT * GO E RR =2 NOLIST NOOECK NOREF NDUSCAR 

AERO NO. 9 COMPRE S SUK BLADE CYCLIC MODAL FLUTTER ANALYSIS 8 

PHlHL=APPEND/A J JL =AP PENJ /F SA VE= APPEND/C A S£YY= APPEND/ CL AM AL ° 
AHP END/U VG=A? pi; ND /k)HHL=APPEND 6 

CEDI l « GEUM2 , /GPL ,EOE XI N ,G PD T ,C S T M, 8G POT , S I L/ V , N . LUS ET/ V,N» 
NOGPOT $ 

LUSET.NOGPOT 8 

ERROR 1 .NUGPDT 8 

GPL.EOEXIN.CPDT.C i TM.BGPDT.SIL 8 
U1JE.D2JE/N03 JE 8 
GEtJ M 2 , fc U E X 1 N/ E C T 8 
CCT 8 

: : '.JM3(E0EX1N»GE0M2/,GPTT /V,N,NOGRAV 8 
GPTT 8 

ECT , E P T , BOPDT , SIL ,GP TT,C STM/EST.GEl ,GPEC Tr«/ V , N« LUSET/ V 0 N , 

NOS l MP /C ,N, 1/ V,N, NOSE NL/V, N, GENE L 8 

NUGENL iNUSMP (GENEL 8 

ERR JK l , NO SIMP 8 

UGPST/GENEL 8 

EST.GPECT.GEI .UGPST 8 

//C.N,ADD/V,M tNUKGvi X /C »N»1 /C ,N,0 8 

//C,N , ADD/V, V ,NOMGG/C ,N,I/C»N,0 8 

// C.N.NOP / V, Y« KGG I N a "i 8 

JMPKGGlN ,KGGI N 8 

//Cifi i ADD / V* N , NOKGG X /C,N,-1 /C.N.O 8 
/KTOT-Afcj, > » » /I , V.LOCA T| UN=-l /C , V ,1 NPTUNI T=D 8 


1 

BEGIN 

2 

FILE 

3 


A 

SAVE 

5 

(COND ^ 

6 

CWtPNT 

7 

PURGE 

e 

<gE> 

9 

CHXPNT 

10 


1 1 

ChKPNT 

12 

( 

13 

SAVE 

1 A 

<CCND JD 

15 

PURGE 

16 

cmkp.nt 

17 

PARAM 

18 

P ARAM 

19 

PAR AM 

20 


21 

(Param3 

22 

^NPUTf^ 
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ORIGINAL RAGE F3 
OF POOR QUALHT 


. RIGID FORMATS 


RIGID FORMAT OMAP LISTING 

AERO APPRUACH, RIGID FURM4 T 9 

LEVEL 2.0 NA S IK AN DKAP COMPILER - SOURCE LISTING 


23 

24 

25 



*6 



27 

2fl 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 


KtoMU,RGGX 8 
KGGX 4 
JKPKCGIN 8 

E ST * C STM ,MP T, 13 I T ,GE0M2 , /Kt LM, KOI CT , HELM, MDI CT , , / V , N , NCKGGX/ V, 
N (NU.1GG/C iN ,/C ,N, /C ,N. /C , 3',C DUP MASS /C , V * CP OAR/C , V „C PRUO/ C,Y, 
CP'JU4J1/C,Y,: PJUAD2/C, Y.CPTRI At /C»Y,CPTRJ A2/C,V ,CPIUBE/ CiV , 
CPQDPLT/C.Y.CPTRPL I /C , V ,C P TRB SC 4 

NOkGGX,NUMGG s 

KELM ,KJICT,MELM,MD1C T 8 

JMPK GCX , ND KG G X 8 

t 

GPEC r.KJIC T,XELM/KuGX,C/ST 8 
KGGX.GPST 4 
JUPKGGX 4 
ERROR l'NUMGG 4 

GPECT ,MDICT,HELM/MGG ,/C ,N,~1 /C, V ,WT MASS' 1 .0 8 
MGG 8 

L GPWG.GKOPNT 8 

BGPDT.C STM.E3EKIN ,HGG/UGPWG/V ,Y ,GRDPNT=-1/C«Y,WTHASS 8 
0 GP WG, ,»,,// 8 
LGPWG 8 

K CGX i KGG/NOGE NL 8 
KGG 8 

LBL 1 1 , NOGENL 8 

GEI , KGGX /KGG/ V,N,LU$E?/V,N,NOGENL/V,N,NOSlHP 8 
KGG 8 
LBL 11 8 

C AS EC C »GE0M4, EQE X I N» f,POT ,BGPDT,CSTf)/RG. .USET, ASET/ V,N, 
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ORIGINAL PfiOE 
OF POOR QUALHY 


COMPRESSOR BLADE CYCLIC MODAL FLUTTER ANALYSIS 
RIGID FORMAT MAP LISTING 


AERO APPROACH, RIGID FORMAT 9 

LEVEL 2.0 NAS TRAN OMAP COMPILER - SOURCE LISTING 


w 

SAVE 

48 

P ARAM 

49 


50 

PURGE 

51 


52 

SAVE 

53 

CFKPWT 
__ ■■ 

54 

QUNuJ} 

55 

<^ONtP> 

56 

CgTs'p 3 

57 

SAVE 

58 

(Tung 3 

59 

(5fp^> 

60 

LABEL 

61 

dpufy^ 

62 

CEKPNT 

63 

{CONG 2 ) 

64 

<$CE1 2 ) 

65 

ChKPN T 

66 

<gCF2 ^ 

67 

CmPiMT 

68 

LABEL 

69 

(Two i Q 


LUSET/V.N.MPIF l/V,N,MPCF2/V,h,Sl NGLE/ VoN jOMIT/ V »N, REACT/ C,N«0/ 
V,\l , A CPE A f/V»NiNUSET/V»N, NOL/ V , N , NQA/C , V ,SU81D 5 

RPC FI .SINGLE, OMIT ,RE AC T, NOSET »MPCF2 , REPEAT, NCL» NO A 8 

//C,N,NJT/V,N ,re acga ta /V,N, REACT 8 

ERROR S,KEACL>A TA 8 

GM, GMD/MPCF 1 /GO ,G OU/OMI T/KFS cQPC/St NGLE 8 
G &JM 4 , tO E X I N, USET /CTCO/ V,Y,CTYPE / V,N,NUGO B 
NUGU 8 
CYCO 8 

ERR JR 6 , .IUGQ 4 
LBL 4 , GE.MEL 8 

GPL , OP ST , USE T , SI L /USP ST / V, N, NOG PST B 
N06P ST 8 
LUL4,N0GPST 4 
OGPST,,,,.// 8 
L BL 4 8 

KGG.KNN/MPCFl /MGG ,MNN/MPCF1 8 

KNN.MNN 8 

LBL2.MPCF1 8 

US6T.RG/CM 8 

GM 8 

USLT,CM,KGG,MGG,,/KNN,MNN, . 8 
KNN , MNN $ 

LBL2 8 

RNN.KFF/SINGLE/MNN.MFF/SINGLE 8 


70 


CMCPNT 


KFF.MFF t 


RIGID FORMATS 


ORlQfim Pt&E 53 

OF POOR QUALITY 


RIGID FORMAT O.MAP LISTING 

Ac«(J APPROACH, RIGID FORMAT 9 

LEVEL 2.0 NASTRAN DMAP COMPILER - SOORCE LISTING 


71 

(£pNlT^> 

72 


73 

CHKPNT 

*7A 

LABEL 

75 

(Tom 

76 

IHXPNT 

77 

(£ONO 2) 

78 

<|mpT^) 

79 

CHKPNT 

80 


at 

CHKPNT 

82 

LABEL 

83 

Cupc 3 



B8 SAVt 
B9 CHKPNT 

90 (CUNtT^ 

91 C ^E AO^ ) 


92 SAVE 

93 CHKPNT 


LOLi, SINGLE 8 

USET.KNN ,MNN, , /KFF ,KF S, ,MFF , , 8 
KFF.KFSiMFF 8 
LUL 3 8 

KFF ,KAA /OMI T/ MFF , HA A /GMI T 8 
K AA , M AA 8 
L0L5.OM1T 8 

USET ,KfF,,,/GO,KAA,KCO,LOCJ,,,,, 8 
GU,KAA 8 

Li SET , GO ,MFF /MAA 8 
M AA 8 
LBL5 8 

OWAM ICS.GPL, SIL, USE T/GPLD, St LO , USE TO , TF POUL, , , , , , EEO, ECDVN/V 
N ,LUSET/V,N .LUSETD/V ,N,NUN L/ V , H ,NOOL T/ V ,N , NUPS DL/V ,N,NOFRL/V 
N.NONLF f/V,N,NUTRL/V,N,N(JEtD/C,fWV,N,NJUE 8 

LUSETO, NUUE ,NOEED 8 

ERROR 2, NOEEO 8 

GO,GJU/NOUE/3H,GMD/NCUE 8 

CYCDiKAA ,MAA, ,, /KKK ,MKK , , , /C,N,FORE / V,V,NSEGS=-1 /V,V, 
KIN0EX=-1 t V, Y,C YCS£0=-1 / C,N,1 / V ,N, N0G0 8 

NOGO 8 

KkK,MKK 8 

ERROR 6 ,N0G0 8 

KKK «MKK , , ,EED , iCA SECC / LA MK , PH I K, ,OEIGS / CiH , MOOES /V,N, 
NEJGV 8 

NtlGV 8 

LAMK ,PHIK» 02 IG S 8 


V 


8 


r% 

OF POOR QUALTTY 

COMPRESSOR BLADE CYCLIC MODAL FLUTTER ANALYSIS 

RIUO FORMAT OMAP LISTING 

AERO APPROACH. RIGID FUR Hit T <3 

LEVEL 2.0 NASTRAN OMAP COMPILER - SOURCE LISTING 


94 

95 

96 

97 
SB 


PARAM 



99 

100 

101 

102 

103 


SAVE 

CHKPNT 



104 

105 

106 



107 
10S 
104 
1 10 
111 


SAVE 

CHKPNT 



112 

113 

114 
1 IS 
116 



/ /C t M jNP Y / V .N .C AKO NO / C.N.O / C.N.O S 
0 Ei GS «L AMR) r a i // V» NgCARDNO 8 
CARDNO 8 
ERROR AfNEI GV 5 

CVCO) i > fPHJKiLAMR /, , , PH l A, LAMA f C.N.B4CK / V.Y, MSECS /V,T, 
KINOEX / V.Y.CYCSEO / C.N.l / V > N « NUGO 8 

N0G3 8 

LAMA tPril A 8 

ERROR 6 .NUGO 6 

USET« fPrilAt tiCUtGMiiKFS) • / PHIG.i / C.N.l t C.N.RC1G 8 

CASECCiCSTMjMPT ,D I TfEQEXINgSl L g t .8G POT .LAMA. . PH 1G. EST . « / n 
OPHI Of i i / CiNfREIG i 

OPHI G it) it// ViNjC ARONO 8 

CAR UNO 8 

EOI.USETiBGPOT.CSTP.EQEXI N,GM,GO / AE RO, ACPT »FL 1ST , GT KA, PV ECT / 
V,N,NK/V,N,NJ/V,Y»MUMACH/V,Y,MAXMACH/V,Y,1REF/V»V,MTYPE/V,N, 
NElGv/V, Y.K INOE K=- 1 8 

NK.NJ 1 

AERO. ACP T.FLi STiGTKA iPVECT 8 
PlilAiPVECT) / PHI AX. i « / C . N . 1 8 

P Hi AX i MAA i PHI A X 1 1 » / MI / C . N .3 /C i N . 1 /C ■ N. 1 1 C. N >0/C i N. 1 8 

CASECC.MATPOOL .EQDYN . , TF P0CL/K2 PP » M2 PP , B2PP/V. N.LUS ETD/V.Ni 
NUK2PP/ViNiN0M2PP/ViN.NUb2PP $ 

NUK2PP jNOMZPP .NU82PP 8 

X 200/NUK.2PP/M2UD/NUM2PP/B200/NOB2PP 8 

M2PP .M2UO/NJSET/B2PP »B2J0 /NGSE T /K2 PP.K2DG/ NOSET 8 

X2PP .M2PP.U2PP .K20U.M2D0 .8200 8 

USETO. GMiGO .. t. . K2PP .M2PP.B2PP/ . . oGMO.GOD.R200. M2 00. B2D0/C. N. 
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OBK#**: 

OF P° 0R 


QUfcUTV 


RIGID FORMATS 


RICIO FORMAT DMAP LISTING 

AERO APPROACH, RIGID FORMAT 9 

LEVEL 2.0 NASTRAN JMAP COMPILER - SOURCE LISTING 


117 

CHKPOT 

i 18 

AN 3 


119 

SAVE 

120 

CHKPNT 

121 

PAKAML 

122 

PURGE 

123 

QunjJ) 

12 A 

<£l_ise£) 

125 

SAVE 

126 

(^rtmT^) 

127 

PARAM 

128 

P Afi AM 

129 

t£QNQ ^ 

1 30 

< guTP Z ) 

131 

SAVE 

132 

<£iTTMs£) 

133 

LABEL 

13 A 

(CuND 

135 

PARAM 

136 

Camg 


CNPLEV/C ,N,OI SP/C ,N,MOOAL/C oN, 0.0/C,N,0.0/C,N,0 „Q/V,N,N0K2PP/V, 
N,N0M2PP/V,N.NU82PP/V,N,MPCFl/V»N,SI NG LE/V , N.ON IT/ V , N, NOUE/ C, 
N.-l/C.N.-l/G ,N,-1/C iN,-l S 

K2DU,H2DD,B23Q »GOD ,GMO 6 

USETJ.PHJAX.MI .LAMK.Ul T, rt2DD,B2D0,K2UD,CASECC / MHH, 8HH, KHH. 

PHI OH / V.N.NUUC/C .Y,LMUL>ES“999999/C.¥,LFREQ=0.0/C( Y,HPRE0=0.0/ 
V,9,NCJM2PP/V,N,NO02PP/V,N ( NOK2PP/V,N,NONCOP/V,H.FMOOE/C«V, 
K0AMP=-1 & 

NUNC UP , FMJOE S 

N HH , 8 HH » KHH , P H l DH $ 

PCOB//C iN,PR£ S/C,N,/C.N,/C.N,/V,N,NUPCOJ 6 
PLTSE IX , PL TPAR .GP SE T SiE L SE IS / NOPCOfl a 
p 2 , nop coa a 

PCOd, CaOYN.e; T / PLTSElX.PLTPAft.GPSETSfELSETS / V.N.NSUI /V.N, 
JUMPPLU T =- 1 S 

NSIL 1, JUMPPLO T & 

PLTSETA // * 

//C ,N ,MP Y/V,N .PLTFLG/C ,N ,1 /C ,N,1 $ 
//C.N.MPY/V.N.PFILE/C.N.O/C.N.O 6 
P 2, J UMPPLO T a 

PL TP AR , uP SETS .ELSE TS .C ASECC . 8 GPD T , EQO YN , . . , ,/PLCT X l/Vo N» NS I L 1 / 
V.V.LUSF J/V,U (JUMPPLUT/V.N.PLTFLG/V.N.PFILE a 

jumpplo t.pltf lg „pf ile t 

plotxi // a 

p 2 a 

ERROR 2.N0EED a 

//C,.I,AUD/V,N .OESTRY/C ,N,0/C .Npl 8 

AERU.ACP T/A JJL .SKJ.DI JK.02 JK/V.N ,NK/V ,N» NJ/V.N, OESf BY S 


1J7 SAVE 


OESTRY a 


OREQHMAi; P&3E 

OF POOR tfUALRY 


COMPRESSOR BLADE CYCLIC MODAL flUTTER ANALYSIS 


RIGID FORMAT DM AP LISTING 


AERO APPROACH. RIGID FORMA T 9 

LEVEL 2.0 NASTKAN DMAP COMPILER - SOURCE LISTING 


138 

139 

140 



141 LABEL 
14 2 PAR AM 
143 < $MP~ D 


144 SAVE 

145 CHKPNT 
1'46 PAR AM 
147 PARAM 
14b PARAM 




159 

160 



A J JL iSKJ.Dl JR [D 2 J K 8 
NODJ E «NUDJE $ 

/01JE ,D2JE, /C « Y »PUSl Tl UN=-1 /C , Y,UNI TNUM=1 1/ C , Y gUSRl ABEL“ 

TAPE ID $ 

NOUJG 8 

//C*N»AOO/ViN|SQHHL/C ,N,1 /C,N,0 8 

A J J L , SK.J .01 JR ,0 2 JR ,G TKA , PH I OH ,01 J£ , 02 JC , US ET 0, AERO/ QHHL 7V, 
N.NOUE/V.N, KJHHL & 

XO'tHL 8 

QHHL & 

//C, N.MP Y/V.N ,NUP/C,<V,-l /C,N,1 8 
//C, \!,MPY/V,N tNOP/CjN.l/C.Ntl 8 
//L.N.MP Y/V.N .NOH/C ,N,0/C,N»1 8 
//C,,M,MPY/V,N ,FL00P/V,Y,NUDJE“-1/C,N,0 6 
LOOP TUP 8 
LOUP TOP 8 

K.HH,dHH,MHH,JHIIL .CASECC ,PLI ST/FSAVE ,KKHH , BKHH, MKHH/ V , N, FLOOP/V, 
N.T START 8 

FLCJUP , T STAR T 8 

KKMH, BXHH, MXHH.EEO, CASECC /PH 1H,C LAMA, OCE IGy V.N.EIGVS 8 
E1GVS 8 

LULZAP.EIGVS 8 
LUL 16.NUH 8 

CASECC, EQOYN, USE TO , PHI H ,C LAMA , , /CPH1H, /C ,N» CE I GEN/ C, N, MODAL /C 5 
N, 123/V,N,NUH/V,N,NUP/V,N,FM00E $ 

NOHgNOP 8 

LBL 16 jNtlH 8 
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RIGID FORMATS 


RIGID FURHAT DHAP LISTING 

AERO APPROACH, RIGID FORMAT <3 

LEVEL 2.0 NAS IRAN OMAP COMPILER - SOURCE LISTING 


161 

3 

0PH 1 H , , i , , / /V » N »C ARD NO 4 

162 

SAVE 

caronu $ 

163 

LAUEL 

L8L 16 4 

16 A 

ann^> 

PHIHtCLAMA ,FSAVE/PHI HL,CLA KAL ,C A SE YV » OVG/V , N, TS T ART / C« V jVREFb 
1.0/C ,Y,PrtlNT»YES8 4 

165 

SAVE 

T STAR T $ 

166 

CMPNT 

PhlHL .CLAf'.ALtCASE WoOVG 5 

167 

(SoNU 

CONTINUE.TSTART 4 

168 

LABEL 

L BL ZAP 4 

169 


CONTINUE .FLOOP 4 

170 

REP1 

LOOP TOP , 100 4 

171 

(JumT^ 

ERROR 3 S 

172 

LABa 

CONTINUE 4 

173 

CHKPNT 

OVG 4 

176 

PAR AML 

XYCUB//C ,N, PRES/C ,N» /C ,N , /C » N ,/ V ,N t NOX YCOB 6 

17!i 

(Tund 

NOXYUUT.NOXYCDb 4 

176 

(XYTR AnT) 

XYCJBfOVG, ,, , /KYPLTCE/C.N, vg /C , N ,PSET/ V, n, pf ile/v , n , c arono 4 

177 

SAVE 

P F IL E jCARUNO 6 

178 

(aVpljt) 

XYPl ICE// 4 

179 

label 

NOXYOLT 4 

180 

PAR AM 

//C 5 N,ANO/V t N,PJUMP/V,N,NOP°-l/V,N, JJMPPLOT 4 

101 

(£pND 

FIVIStPJUMP 4 

182 

<5pOACC) 

CASE YY,CLAMAL ,PHl HL ,CASECC , , /CL AMAU ,CPHIH1,CAS E22, »/C,N» 
CEIGN 4 

183 

(Tcrl'3 

CPHIH1,PHI0H/CPHID 4 

186 

ChKPNT 

CPHIO 4 
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OF POOR QUALITY 


COMPRESSOR BLADE CYCLIC MODAL FLUTTER ANALYSIS 
RIGID FORMAT UMAP LISTING 


AERO APPROACH, RIGID FORMAT 9 

LEVEL 2.0 NAS TRAN l)NAP COMPILER “ SOURCE LISTING 


185 

<gm£} 

186 

(£UNU 3 

167 

CscTO 

IBS 

LABEL 

189 

CHXPNT 

190 

(equiO 

191 

CHonu D 

192 

<20 

193 

<£ARlO 

199 

LABEL ' 

195 

dcR2 Z) 

196 

CHKPNT 

197 

(GfP 3) 

198 


199 

(£LOT 3 

200 

Cprt^sg) 

201 

LABEL 

202 

(Tump 

203 

LABEL 

209 

^RTP AMft 

205 

LABEL 

206 

^RTPAKjj) 

207 

LABEL 

208 

(gjUPARjfo 


CPHIU tCPMIP /NOA 8 
LULi9„.\IOA $ 

USETD , bCPHID, , ,GOO,GMQ, , Kf S, , /C PHI P » , QPC/C, N,l/ C» N» OY NAM ICS 6 

LBL 19 8 

CPHIP,UPC 8 

CPHIU,CPHIA/NOUE s 

LBLNUE ,NUUE 8 

U$ETO/RP/C,N»D/C,N,A/C,N,E 8 
CPHIU, , RP/CPH IA ,, ,/C ,N, 1/C ,Ns3 8 
LBLNUE 8 

CASEZ2,CSTM,MPT,DI T.EQOVN.SI CD, , , BGPOT , CL A MALI , 0PC , CPHI P, ES T, , / 
,UgPCI,UCPHIP ,UE SCI ,UEFul , PC PHI P/C,N,CFIGN 8 

PCPH1P 8 

UCP H 1 P , UQPC 1, OE SC 1 ,UEFC1 , , // V ,N ,CARDNO 8 
P 3, J UMPPLOT 8 

PLTPARcGP SETS ,ELSE IS, CASE 22, BGPOT, EOOVN, SILO,, PCPHI P, ,/ PLOT S 3/ 

v,n,nsili/v,\,luset/v,n,jumpplot /v,n,pltflg/v,n,pfile 8 

PLDTX3// 4 
P 3 8 
FINIS 8 
ERROR 1 8 

//c ,m »- i /c , n, f sub son 8 

ERROR 2 8 

//C,No-2/C,N'FSUaSCN 8 
ERROR 3 8 

//CoNi-3/C,N,FSUBSCN 8 
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ORfQfNAl’ & 
OF POOR QUALITY 


RIGID FORMATS 


RIGID FORMAT OMAR LISTING 


AERO APPROACH, RIGID FORMAT 9 



LEVEL 2 . 

0 NASTRAN OMAP COMPILER - 

SOURCE 

209 

LABEL 

ERROR *> 5 


210 

(gjuPAHJ*) 

//C i H (- 6/C jN, F SUB SON 8 


211 

LABEL 

ERROR 5 8 


212 

(gRTPAKiP 

n c,n,-a / : ,n,:.yc modes 

8 

213 

LABEL 

ERROR 6 5 


2H 

(gwTPAR^ 

n C.Ut-5 / : ( Ni,C VC MODE s 

$ 

2 15 

LABEL 

FINIS 8 


216 

ENC 

1 



9®MCI ERRORS FOUNO - EXECUTE NASTRAN PROGRAM®® 
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COMPRESSOR BLADE CYCLIC MODAL FLUTTER ANALYSIS 


2 . 6.2 Description of DMAP Operations for Compressor Blade Cyclic Modal 
Mutter Analysts ~ — ' 

3. GP1 generates coordinate system transformation matrices, tables of grid 
point locations, and tables for relating Internal and external grid 
point numbers. 

5. Go to OMAP No. 203 and print error message if no grid points are present, 

8. GP2 generates Element Connection Table with internal Indices. 

10. GP3 generates Static Loads Table and Grid Point Temperature Table. 

12. TA1 generates element tables for use in matrix assembly and stress 

recovery. 

14. Go to DMAP No. 203 and print error message if no elements hove been 
defined. 

20. Go to DMAP No. 25 if stiffness matrix is not user input. 

21. Set parameter NOKGGX ° -1 so that the stiffness matrix will not 

be generated in DMAP No. 26. 

22. INPUTT1 reads the user supplied stiffness matrix from tape (GINO 
file INPT ) . 

23. Equivalence [K* g ] to Ug g ]. 

26. EMG generates structural element matrix tables and dictionaries for 
later assembly. 

29. Go to DMAP No. 32 If no stiffness matrix is to be assembled. 

30. EMA assembles stiffness matrix [K gg ] and Grid Point Singularity Table. 

33. Go to DMAP No. 203 and print error message if no mass matrix exists. 

34. EMA assembles mass matrix [M gg ]. 

36. Go to DMAP No. 39 if no weight and balance request. 

37. GPUG generates weight and balance information. 

38. 0FP formats weight and balance information and places it on the 
system output file for printing. 

40. Equivalence [K* 3 to [K ] if no general elements. 

42. Go to DMAP No. 45 if no general elements. 

43, SMA3 adds general elements to [K gg ] to obtain stiffness matrix [ 3 . 

46. GP4 generates flags defining members of various displacement sets 

(USET), forms multipoint constraint equations {R 0 ]{u } » 0. 

9 9 

49. Go to DMAP No, 211 and print error message if free-body supports 
are present. 

51. 6PCYC prepares segment boundary table. 

54. Go to DMAP No. 213 and print error message if C V JO I N data is Inconsistent. 
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ORIGINAL PASS E8 
OF POOR QUALITY 


55. 


56. 


58. 


59. 


61. 


63. 


64. 


66 . 


69. 


71. 


72. 


75. 


77. 


RIGID FORMATS 


Go to DMAP No. 60 If general elements present. 

GPSP determines If possible grid point singularities remain. 

6o to OMAP. No. 60 if no grid point singularities remain, 

0FP formats the table of possible grid point singularities and places 
it on the system output file for printing. 

Equivalence [K ] to £k ] and [«_„] to [M nn ] if no multipoint 
constraints. 99 99 n 


Go to DMAP No. 68 if MCE1 and MCE2 have already been executed for 
current set of multipoint constraints. 

MCE 1 partitions multipoint constraint equations [R 1 = [R R ] and 

9 m ^ n 

solves for multipoint constraint transformation matrix [G^] » 

MCE2 partitions stiffness and mass matrices 

i 1 


[> W ■ 


nn I nm 


mn 


mm 


and 


[ V 


nn 


mn 


nm 


mm 


and performs matrix reductions 


l‘,.3 • I R „n3 * 

t" BB ) ■ [»„„) * 

Equivalence [K nf) ] to [K^] end [M nn ] to if no single-point 

constraints. 

Go to DMAP No. 74 If no single-point constraints. 

SCEl partitions out single-point constraints. 


' K ff • K f S ' 


M ff j M fs 

i 


* 

i 

and [M nn ] - 

« 

V { K ss. 


M sf { «ss 


Equivalence [K^j to [lt aa ] and [M^J to lM aa J if r omitted degrees 


of freedom. 

Go to DMAP No. 82 if no omitted coordinates. 
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78 . 


80 . 


83. 


85. 


86 . 


87. 


90 . 


91 . 


95 . 


97. 


98. 


SMP I partitions constrained stiffness matrix 




aa 


ao 


i. - - 

l K «« ' K a* 

oa i oo 


and solves for transformation matrix ° -[K I'Vk ] 

v 0 00 Ofl 

and performs matrix reduction D [tT 1 + [K T 1[G ]. 

a a Afl - *- o a - o 


SMP 2 partitions constrained mass matrix 


fa 


Ltoff] 


aa i 

* 


oa { 


M 


ao 


00 


and performs matrix reduction 


W..1 ■ [".,3 1 K.JIG,] * IgJHH ][G„] * £sj3CH„ 0 J. 


OPD generates flags defining members of various displacement sets used 
in dynamic analysis (USETD), tables relating internal and external 
grid point numbers, including extra points introduced for dynamic 
analysis, and prepares Transfer Function Pool and Eigenvalue Extraction 
Data. 


Go to DMAP No. 205 and print error message if no Eigenvalue Extraction 
Data. 

Equivalence [G Q 3 to [G^] and [G^J to [Gjjj] if no extra points introduced 
for dynamic analysis. 

CYCT2 transforms matrices . rcm symmetric components to solution set. 

Go to DMAP No. 213 and print error message if CYCT2 error was found. 
READ extracts real eigenvalues from the equation 

i‘n - ■ 0 • 

and normalizes eigenvectors according to one of the following user 
requests: 

1) Unit value of selected coordinate 
21 Unit value of largest components 
3) Unit value of generalized mass. 

0FP formats eigenvalues and summary of eigenvalue extraction infor- 
mation and places them on the system output file for printing. 

Go to DMAP No. 209 and exit if no eigenvalues found. 

CYCT2 finds symmetric components of eigenvectors from solution set 
eigenvectors . 
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Go to DMAP No. 213 and print error message If CYCT2 error was found 
SDR 1 recovers dependent components of the eigenvectors 



104. 

106. 


109. 

no. 


1 

I 


111 . 

114. 


116 . 


SDR2 prepares eigenvectors for output (0PHIG). 

0FP formats tables prepared by SDR2 and places them on the system 
output file for printing. 

APDB processes the aerodynamic data cards from EDT. AERO and ACPT 
reflect the aerodynamic parameters. PVECT is a partitioning vector 
and GTKA is a transformation matrix between aerodynamic (K) and 
structural (a) degrees of freedom. 

PARTN partitions the eigenvector into all sine or all cosine 
components . 

SMPYAD calculates modal mass matrix 

w * CO T r« aa ) u* a ] 

MTRXIfl selects the direct Input matrices [K pp ], E M p p 3 . and [C pp ]. 

Equivalence [H* p ] to [M^], [E* p J to [0^] and [K* p ] to [K* d J If no 
no constraints applied. 

2 2 

GKAD applies constraints to direct input matrices [K pp ]„ [M pp J, and 

and [B dd 3 (see Section 9.3.3 of the Theoretical Manual) and 
forms [G md ] and [G fl .]. 


1 


118 . 


123. 

124. 

126. 

129. 

130. 
132. 

134. 

136. 

139. 

140. 
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GKAM selects eigenvectors to form U dh ] and assembles stiffness, matrices 
and damping matrices In modal coordinates: 


CK hh J 

"W 

f'hiA 





where 

KDAMP n 1 

° modal masses 
b f » m i 2 tt f t g ( f ^ ) 
kj c m^ 4 ti* f^ 


KDAMP ° -1 (default ) 
m^ = modal masses 

bj ° 0 

k 1 c (1+fg(f 1 )> 4T! S f*m j 


Go to DMA? No. 133 if no plot package Is present. 

PLTSET transforms user input into a form used to drive structure 
plotter, 

PRTMSG prints error messages associated with structure plotter. 

GO to DMAP No. 133 if no undeformed aerodynamic structure plot request. 

PL 0T generates all requested undeformed structure plots. 

PRTMSG prints plotter data and engineering data for each undeformed 
aerodynamic plot generated. 

Go to DMAP No. 205 and print error message if no Eigenvalue 
Extraction Data. 

AMG forms the aerodynamic materix list [Ajj], the area matrix 


and the downwash coefficients [Dj k 3 and [D^j. 

Go to DMAP No. 141 if no user- supp 1 1 ed downwash coefficients. 


I N PUTT 2 provides the user-suppl led downwash factors due to extra 
points , C Oj e 3 ) . 
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154. CEAD extracts complex eigenvalues from the equation 

K»i> ! * 0 ♦ K S h 3<» h > ■ 0 

and normalizes eigenvectors to unit magnitude of largest component. 

156. Go to DMAP No. 168 If no complex eigenvalues found. 

157. Go to DMAP No. 163 If no output request for the extra points Intro- 
duced for dynamic analysis or modal coordinates. 

158. VDR prepares eigenvectors for output, using only the extra points 
Introduced for dynamic analysis and modal coordinates. 

160. Go to DMAP No. 163 If no output request for the extra points 
Introduced for dynamic analysis or modal coordinates. 

161 0FP formats eigenvectors for extra points Introduced for dynamic 

analysis and modal coordinates and places them on the system output 
file for printing. 

164. FA2 appends eigenvectors to PH1KL, eigenvalues to CLAHAL, Case Control 
to CASEYY , and V-g plot data to 0 V G . 

167, Go to DMAP No. 172 if there Is insufficient time for another 
flutter loop. 

169. Go to DMAP No. 172 If flutter loop complete. 

171. Go to DMAP No. 207 for additional aerodynamic configuration triplet 

values . 

175, Go to DMAP No. 179 If no X-Y plot package is present. 

176, XYTRAN prepares the input for requested X-Y plots. 

178. XYPL0T prepares requester X-Y plots of displacements, velocities, 

accelerations , forces , y cresses , loads or single-point forces 
Of constraint vs. t . . 

181. Go to DMAP No. 215 if no output requests Involve dependent degrees 
of freedom or forces and stresses. 

182. M0DACC selects a list of eigenvalues and vectors whose Imaginary parts 
(velocity in Input units) are close to a user Input list. 

183. DDR1 transforms the complex eigenvectors from modal to physical 
coordinates 

° I*dh^*h^ * 

185. Equivalence to [$'] If no constraints applied. 

186 Go to DMAP No. 188 if no constraints applied. 
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190, 

191, 

192. 

193. 


195. 

197. 

198. 

199. 

200 . 

202 . 

204. 

208. 
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SDR1 recovers dependent components of eigenvectors 


(♦J) ■ lejluj) 


° Uf ♦ 


►5 ♦ # c 
f . _®. 

* c 


{< + ♦£) . 


u c } « rG d j<4> c ♦ * c ) 

1 n» J v n ' 


!♦? ♦ - uf) 


and recovers single-point forces of constraint { q ) o 


Cicf 6 3U f ), 


fQp> 


Equivalence [$|j] to [$£3 If no extra points Introduced for dynamic 
analysis . 


Go to DMAP No. 194 If no extra points present. 


VEC generates a d-size partitioning vector (RP) for the a and e sets, 
PARTN performs partition of [$£] using RP. 



SDR2 calculates element forces and stresses (0EFC1, 0ESC1) and 
prepares eigenvectors and slnqle-polnt forces of constraint for 
output (0CPIILP, 0QPC1). It also prepares PCPHIP for deformed plotting. 

0FP formats tables prepared by SDR2 and places them on the system 
output file for printing. 

Go to DMAP No. 194 If no deformed structure plots are requested. 

PL0T prepares all deformed structure plots, 

PRTMSG prints plotter data and engineering data for each deformed 
plot generated. 

Go to DMAP No. 215 and make normal exit. 


M0DAL COMPLEX EIGENVALUE ANALYSIS ERR0R MESSAGE N0. 1 - MASS MATRIX 
REQUIRED F0R M0OAL F0RMULATI0N. 


W0DAL C0MPLEX EIGENVALUE ANALYSIS ERR0R MESSAGE N0. 2 - EIGENVALUE 
EXTRACTION DATA REQUIRED FOR REAL EIGENVALUE ANALYSIS. 
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208. M0DAI. COMPLEX EIGENVALUE ANALYSIS ERR0R MESSAGE N0. 3 - ATTEMPT TO 
EXECUTE MORE THAN 100 LOOPS. 

210. M0DAL COMPLEX EIGENVALUE ANALYSIS ERROR MESSAGE N0 . 4 - REAL EIGEN- 
VALUES REQUIRED F0R MODAL FORMULATION. 

212. NORMAL MODES WITH CYCLIC SYMMETRY ERROR MESSAGE NO . 4 - FREE BODY 
SUPPORTS NOT ALLOWED. 

214. NORMAL MODES WITH CYCLIC SYMMETRY ERROR MESSAGE N0. 5 - CYCLIC 
SYMMETRY DATA ERROR. 
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2.6.3 Output for Compressor Blade Modal Flutter Analysts 

The Real Eigen value Summary Table and the Real Eigenvalue Analysis 
summary, as described under Normal Mode Analysis, ore automatically printeo. 

All real eigenvalues are Included even though all may not be used In the 
modal formulation. 

The grid point singularities from the structural model are also output, 

A flutter summary for each value of the configuration parameters Is 

printed out If PR INT = YESB . This Shows p, k, l/k,<r, a «v V a and f 

sound * ’ s ^ na T 

for each complex eigenvalue. 

V-g and V-f plots may be requested by the XYflUT control cards bv 
specifying the curve type as VG. The "points” are loop numbers and the 
“components" are G or F. 

Printed output of the following types, sorted by complex eigenvalue 
root number (S0RT1 ) and (m, k, p) may be requested for all complex eigenvalues 
kept, as either real and imaginary parts or magnitude and phase angle 
(0° - 360° lead): 

1. The eigenvector for a list of PHYSICAL points (grid points, 
extra points) or S0LUTION points (modal coordinates and extra points). 

2. Nonzero components of the single-point forces of constraint for 
a list of PHYSICAL points. 

3. Complex stresses and forces In selected elements. 

The ^FREQUENCY case control card can select a subset of the complex eigenvectors 
for data recovery. In addition, undeformed and deformed shapes mav be reouested. 
Undeformed shapes may Include only structural elements. 
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2 . 6.4 Case Control Deck, and Parameters for Compressor Blade Cyclic 


Modal flutter Analysts 


1. Only one subcase is allowed 


2. Desired direct input matrices for stiffness t^ nD ], boss t» 2 B J. ^d 

rr 


damping £8^ J must be selected via the keywords K2PP, M2PP, or 

rr 


B2PP . 


3. CMETH0O must be used to select an EIGC card from the Bulk Data Deck. 


FMETH0D must be used to select a FLUTTER card from the Bulk Data Deck, 


METH0D must be used to select an EIGR card that exists in the Bulk 


Data Deck. 


6. SCAMPING must be used to select a TABDMP1 table if structural damoina 


is desired. 


7. An SPC set must be selected unless the model Is a free bodv or all 


constraints are specified on GRID cards, Scalar Connection Cards or 


with General Elements. 


Each NASTRAN run calculates modes for only one symnetry index, It. 


The following user parameters are used in Compressor Blade Cyclic Modal 
Flutter Analysis. 

1. GRDPNT - optional - A positive integer value of this parameter will 
cause the Grid Point Weight Generator to be executed and the result- 
ing weight and balance Information to be printed, All fluid related 


masses are Ignored, 


2. WTMASS - optional - The terms of the structural mass matrix are 
multiplied by the real value of this parameter when they ore 
generated in SMA2. Not recommended for use in hvdroelastic 
problems . 

C0UPMASS - CPBAR . CPR0D. CP0UAD1 , CP0UAD2 . CPTRIA1 , CPTRIA2 , 


optional - These parameters 


will cause the generation of coupled mass matrices rather than 
lumped mass matrices for all bar elements, rod elements, and plate 
elements that include bending stiffness. 
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LFREQ and HFREQ » required unless LM0DES Is used. The real values 
of these parameters give the frequency range (LFREQ Is lower 
limit and HFREQ Is upper limit) of the modes to be used In the 
modal formulation. To use this option, LM0DES must be set £o 0. 
LMBDES - used unless set to 0. The Integer value of this parameter 
Is the number of lowest modes to be used In the modal formulation. 
The defalult value will request all modes to be used. 

NSDJE - optional In modal flutter analysis. A positive Integer 
of t h 1 !• parameter Indicates that user supplied downwash matrices 
due to entra points are to be read from tape via the INPwTTE 
module In the rigid format. The default value Is -1. 

PI , P2 and P3 - required In modal flutter analysis when using 
N0DJE parameter. See Section 5.3.2 for tape operation parameters 
required by INPUTT2 module. The defaults for PI, P2, and P3 are 
-1, 11 and TAPEID , respectl ve ly . 

VREF - optional In modal flutter analysis. Velocities are 
divided by the real value of this parameter to convert unit? 
or to compute flutter indices. The default value Is 1.0. 

PRINT - optional In modal flutter analysis. The BCD value N@, 
of this parameter will suppress the automatic printing of the 
flutter summary for the k method. The flutter summary table 
will be printed if the BCD value Is YES for wing flutter, or 
YESB for blade flutter. The default is YES. 

CTYPE - required - the BCD value of this parameter defines 
the type of cyclic sy«metry as follows: 

(1) R0T - rotational symmetry 

(2) DRL - dihedral symmetry, using right and left halves 

(3) DSA - dihedral symmetry, using symmetric and anti- 
symmetric components 

filSEGS - required - the Integer value of this parameter is the 


number of Identical segments In the structural model. 


COMPRESSOR BLADE CYCLIC MODAL FLUTTER ANALYSIS 


CVCSEQ - optional - the Integer value of this parameter specifies 
the procedure for sequencing the equations In the solution set. 

A value of +1 specifies that all cosine terms should be sequenced 
before all sine terms, and a value of *1 for alternating the 
cosine and sine terms. The default value Is -l, 

KINDEX - required In compressor blade cyclic modal flutter analysis. 
The Integer value of this parameter specifies a single value of 
the harmonic Index. 

MINMACH - optional in blade flutter analysis. This is the minimum 
Mach number above which the supersonic unsteady cascade theory Is 
valid. The default Is 1.01. 

MAXMA d - optional in blade flutter analysis. This is the maximum 
Mach number below which the subsonic unsteady cascade theory Is 
valid. The default value Is 0.80. 

IREF - optional In blade flutter analysis. This defines the 
reference streamline number. IREF must be equal to a SLN on a 
STREANL2 built data card. The default value, -1, represents the 
streamsurface at the blade tip. If IREF does not correspond to 
a SLN, then the default will be taken. 

MTYPS - optional in cyclic modal blade flutter analysis. This 
controls J hich components of the cyclic modes are to be used in 
the modal formulation. MTYPE ° SINE for sine components and 
MTYPE ° COSINE for cosine components. The default BCD value is 
COSINE. 

KGGIN - optional in blade flutter analysis. A positive Integer 
of this parameter indicates that the user supplied stiffness 
matrix Is to be read from tape (6IN0 file INPT) via the 
INPUTT1 module in the rigid format. The default is -1. 
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Aw 


3>l.l EDT (TABLE) 


fj 

n 


i 

•< 




Card Types and Header Information : 


Card Type 

Header Word 1 
Card Type 

Header Word 2 
Trailer Bit Position 

Header Word 3 
Internal Card Number 

AEFACT 

4002 

40 

273 

AER0 

3202 

32 

265 

CAER01 

3002 

30 

263 

CAER02 

4301 

43 

301 

CAER03 

4401 

44 

302 

CAER04 

4501 

45 

303 

CAER05 

. 5001 

50 

309 

DEF0RM 

- 104 

1 

81 

FLFACT 

4102 

41 

274 

FLUTTER 

3902 

39 

272 

MKAER01 

3802 

38 

271 

MKAER02 

3702 

37 

270 

PAER01 

3102 

31 

264 

PAER02 

4601 

46 

304 

PAER03 

4701 

47 

305 

PAER04 

4801 

48 

306 

PAER05 

5101 

51 

310 

SET1 

3502 

35 

268 

SET2 

3602 

36 

269 

SPLINE! 

3302 

33 

266 

SPUNE2 

3402 

34 

267 

■SPLINE3 

4901 

49 

307 

STREAHL 1 

3292 

92 

292 

STREAML2 

3293 

93 

293 

VARIAN 

4202 

42 

290 


Card Type Formats : 


AEFACT (Open Ended) 

SID 

FI 

F2 


etc. 

-1 


AERO (6 words) 

ACS ID 

VS0UND 

BREF 


RH0REF 

SVMXZ 

SYMXX 

CAER01 (16 words) 

PID 

CP 

NSPAN 


NCH0RD 

LSPAN 

LCH0RD 


0 

XI 

Y1 


Z1 

XI 2 

X4 


Y4 

Z4 

X43 

CAER02 (16 words) 

EID 

PID 

CP 


NSB 

MINT 

LSB 


LONT 

IGID 

XI 


Y1 

Z1 

XI 2 


CAER03 (16 words) 

EID 

PID 

LP 


LISTW 

LISTC1 

LISTC2 


» , , 

♦ . • 

XI 


Y1 

Z1 

XI 2 


X4 

Y4 

Z4 


X43 
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Card Type Formats (Cont.): 


SET2 (8 words) 

SID 

EID 

SP1 


SP2 

CHI 

CH2 


Z1 

Z2 


SPLINE! (6 words) 

EID 

CAER0 

B0X1 

B0X2 

SETG 

OZ 

SPLINE2 (10 words) 

EID 

CAER0 

B0X1 

B0X2 

SETG 

DZ 


DT0R 

CIO 

DTHX 


DTHV 



SPLINE3 (Open Ended) 

SID 

CAER0 

UFID 


C0MP 

G1 

Cl 


A1 

• • • 

GM 


CM 

AM 

-1 

STREAML 1 (open ended) 

2LN 

G1 

G2 

G3 

G4 

G5 


G6 

, , , 

Gn 


-1 



STREAML2 (10 words) 

2LN 

N0TN0 

STAGGER 


CH0RD 

RADIUS/DCBDZ8 

B0PACE 


MACH 

FLOWA/ SWEEP 

DEN 

VEL 

VAR IAN (Open Ended) 

DB l z 

QBl 2 

etc. 
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j?, i .Jl. Data Blocks Output from Module APDB 
2. I .1 .1 AERO (Table) 

Description 

See description and format of AERO table 
2, hi -2 FUST (Table) 

Descrl ptlon 

See description and format of FUST table 
2,1 „ 2. .3 GTKA (Matrix) 

Description 

See description and format of GTKA matrix 
J. I.X.4 PVECT (Matrix) 

Descrl ptlon 

{ PVECT } - Partitioning vector for cyclic modes. 


Matrix Trailer 

Number of columns = 1 

Number of rows = NEIGV (for KINDEX > 0, 2 ■ NEIGV) 

form ° rectagular 

Type ° real-single precision 


preceding page blank 


not filmed 
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3- 1.3 ACPT (Table) 


Description 

Aerodynamic connection and property table for compressor or turboprop blades. Contains 
one record for each compressor/or turboprop blade. 


Table Format 


Record 

0 

1 


or 1 


Word Type 
1-2 B 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 


I 

1 

R 

R 

I 

I 

r 

i 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

n 

12 

13 

14 

15 

16 

17 

18 
19 


I 

I 

R 

R 

I 

I 

I 

I 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 


2 


Table Trailer 


Item 

Data block name (ACPT) 


Key word, 6 for compressor blades 
IREF parameter 
MINMACH parameter 
MAXMACH parameter 

Number of blade streamlines, NLINES 

Number of stations on blade, NSTNS 

Streamline number, SLN 

Number of stations on streamline, NSTNSX 

Stagger angle, STAGGER 

Chord length, CHORD 

Radius of streamline, RADIUS 

Blade spacing, BSPACE 

Mach number, MACH 

Gas density, DEN 

Flow v locity, VEL 

Flow angle, FLOWA 

X-coordinate, basic] REPEAT 

Y-coordinate, basic \ NSTNS 

Z-coordinate, basic j TIMES 


REPEAT 
> NLINES 
I TIMES 


Key word, 7 for turboprop blades 
IREF parameter 
(1INMACH parameter 
MAXMACH parameter 

Number of blade streamlines, NLINES 

Number of stations on blade, NSTNS 

Streamline number, SLN 

Number of stations on streamline, NSTNSX 

Stagger angle, STAGGER 

Chord length, CHORD _ _ 

Radius of streamline .or 3C/3Z, RADIUS/DCBDZB 
Blade spacing, BSPACE 
Mach number, MACH 
Gas density, DEN 
Flow velocity, VEL 

Flow Angle or Sweep Angle, FLOWA/ SWEEP 
X-coordinate, basic 1 REPEAT 
Y-coordinate, basic > NSTNS 
Z-coordinate, basic J TIMES 



Additional records for other blade 


Word 1 = 1 

Word 2-6 = zero 
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3. i .1 ACPT (Table) (Cont'd.) 


Notes 

1. Words 7-19 are repeated for each streamline. There are NUNES streamlines 
and they are from the blade root to the blade tip. These data Items are 
taken from the STREAML2 bulk data cards. 

1 . Words 17-19 are repeated for each node on the streamline. There are 
NSTNS triplets (X, Y, Z). They are from the blade leading edge to the 
blade trailing edge. 

3. All records 1 to N will be all Keyword 6 or Keyword 7 and may not be mixed. 
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S.X FUNCTIONAL MODULE AMG (AERODYNAMIC MATRIX GENERATOR) 
y.i.l Entry Point : AMG 
3,2 .2 Purpose 

To generate a list of aerodynamic influence matrices (AJJL) and the transformation matrices 
needed to convert these to the interpolated structural system (SKJ, D1JK, D2JK). 

2, X -3 DMAP Calling Sequence 

AMG AER0,ACPT/AJJL,SKJ,S1JK,DZOK/V , N ,NK/V ,N,NJ $ 

j.i.4 Input Data Blocks 

AER0 - Aerodynamic matrix generation data 
ACPT - Aero connection and property data. 

Note: Neither AER0 or ACPT may be purgeo. 

j , L .S Output Data Blocks 

AJJL - Aerodynamic influence matrix list 
SKJ - Integration matrix list 
D1JK - Real part of downwash matrix 
D2JK ■ Imaginary part of downwash matrix 

3. X .6 Parameters 

NK - Input - integer - no default, number of degrees of freedom in k-set 
NJ - Input - integer - no default, 

N sb N 1b 

3 D i«i 1 j»i 1 

where N b = number of aero boxes 

N sb = number of slender bodies 
N.. = number of interference bodies 

ID 

f^ = degrees of freedom for each slender body 
g 1 = degrees of freedom for each interference body 
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.7 Method 

Module AMG Is broken into two sections. Section one outputs AJJL and SKJ and Section two 
outputs D1JK and D2JIC. Each section has a branch on method to output columns of the matrices. 
Each section has a common block to communicate between the module driver and the method dependent 
code. (See description of design requirements.) 
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MODULE FUNCTIONAL DESCRIPTIONS 

The flow for Section one is as follows: Four buffers are allocated from the bottom of core 
and ACPT, AER0, AJJL, and SKJ are opened. Record 1 of AER0 Is read Into /AMGMN/ at ND and the 
header and trailer for AJJL are initialized by passing the ACPT and counting the method groups. 

The rest of Section one Involves a loop where pairs of m and k are read into /AMGMN/ and 
then each record in the ACPT is processed. For each record of ACPT processed, one word (the 
method) is read and then a branch on method is taken. Each method is expected to: 

1. Read its record of the ACPT and leave it positioned to read the next record. 

2- Outout columns of AJJL and SKJ of proper size. 

The row number to start at is NR0W+1 for AJJL and ISK for SKJ. 

3. Always increment NK0W, ISK and NSK by the number of rows added. 

Once an end of file is reached in AER0, AER0, AJJL and SKJ are closed and Section two starts. 

The flow in Section two is as fellows: Three buffers are allocated and DlJK and D2JK are 
opened (ACPT is left open). The trailers are initialized and the method of ACPT is read. The 
program branches on method and then loops until all records on the ACPT have been processed. When 
an E0F is reached the files are closed and their trailers written. 

J. 1-.7.1 Doublet Lattice Method without Bodies 

The flow for Section one of the Doublet Lattice method is as follows: Subroutine DLAMG is 

the driver for this method. It reads in the ACPT record for this method and sets up the pointers 
to the various arrays in common DLC0M. Columns of SKJ are output. If there is enough core avail- 
able, GEND is called to output one matrix of the AJJL list. When GEND is through, DLAMG bumps 
NR0W and returns. 

Subroutine GEND outputs a row of the AJJL matrix for each box on the CAER01 element. To do 

this, it picks up the proper strip and panel data and calls DPPS once for each box. A row is 

packed out after each call to DPPS. 

Subroutine DPPS is also in a loop. There is one row element for each box and DPPS prepares 

the variables necessary for the computation of each element and calls SUBP to calculate the 

element. When the row is done DPPS returns to GEND. 

Subroutine SUBP computes the downwash factor elements by calling subroutines SNPDF and INCR0 
to compute the indicated components that make up the element. SNPDF computes the steady downwash 
factors and INCR0 computes the unsteady downwash factors for one receiving-sending point combination. 
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Each combination has four influence quadrants (upper left, upper right, lower left, lower right), 
so these routines must be called four times for each element and then the result summed before 
SUBP returns. Subroutine 1NCR0 uses subroutines TKER, IDF1, and IDF2 to compute the final result. 

The flow for Section two of the Doublet Lattice method is as follows: Subroutine DLPT2 
prepares all the computations necessary. DLPT2 reads the record of ACPT and then loops through 
each box packing out a column of D1JK and D2JK for each box. 

The row position of each pair of values for a column is 2*(box number-1) + 1 . Successive 
rows of SKJ (output in Section one) have the following form: 


SKJ ■ 


2.0 * EE str1p * DELX bM 


EE strtp * “ EL1 W / 2 ° j 


Successive rows of D1JK have the following form: 


DIJKh 


0.0 

1.0 


Successive rows of D2JK have the following form: 


D2JK-o 


-2.0/REFC 


DELX faox /2.0*REFC 


( 1 ) 


( 2 ) 


(3) 


3, X .7.2 Doublet Lattice Method with Bodies 

The flow for Section one of the Doublet Lattice method with bodies is as follows: Subroutine 

DLAMBG is the driver for this method. It reads the ACPT record and sets up the pointers to the 
various arrays in common DLBDY . The flow then depends on the type of problem submitted. If 
panels exist then GENDSB is called to build part of AJJL on a scratch file. If panels and slender 
bodies are used, AMGR0D is called to build this combination on a scratch file. Then AMGSBA is 
called to output this method's part of AJJL. AMGBFS is called to build SKJ for this method. 


Up to four additional buffers may be needed by this method. 

Figure one shows the subroutines and scratch files used to put AJJL together. 
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NTP 

NTZ 

NTY 

NTZS 

NTYS 


SCR5 


DPPSB 

DZPY 

DYPZ 

DZYMAT 

SCRl 

SCR2 

SCR3 

SCRl 

SCR2 

DPZY 





SCRl 





DPZY 





SCR4 





0 

0 

0 

AMGSBA 

0 

0 

0 




AJJL ' 


NTP, NTZ, NTY, NTZS, and NTYS 
are lengths. 

SCR1 - SCR5 are scratch files. 
Other names are subroutines. 


Figure 1. Subroutines and Scratch files for AJJL 


Subroutine GENDSB calls DPPSB once for each row on SCR1 (NTP times), and DPZY once for each 
NTZ and NTY rows on SCR1 and SCR4. Then DZPY once for each row on SCR2 and DYPZ once for each 
row on SCR3. Once this process is done, GENDSB builds the data collected from SCR1 , SCR2, SCR3 
and SCR4 on SCR5, 

Subroutine AMGR0D calls OZYMAT once to build SCR1 and once to build SCR2. Then AMGSBA is 
called to put SCR1 , SCR2 and SCR5 together plus AMGSBA's part of AJJL onto AJJL. 

Matrix SKJ is built by AMGBFS. It calls BFSMAT, builds some data on SCR1 , then 
SKJ is built and packed out. (See the Theoretical Manual for description of SKJ.) 

Figure 2 shows the various calls that may take place during Section one of this method. 

For Section two DLBPT2 is called and the ACPT is read. Then the pointers to the necessary 
arrays are computed and D1JK and D2JK are packed out. D1JK and D2JK look like the Doublet Lattice 
Method without Bodies except D1JK has a -1. instead of a 1.0 for D1JK on y-slender elements, 
and D2JK has a 0.0 for slender elements instead of DELX/2*REFC. 
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] 

j 


Figure 2. Section One Calls for the Doublet Lattice Method with Bodies 


3 , 1..7.3 Mach Box Method 


The flow for Section one of the Mach Box Method is as follows: Subroutine MBAMG is called to 
read the ACPT and set up pointers to arrays in common MBAMGX. MBAMG then makes calls to various 
subroutines to get AJJL built. Most of the arrays used by the Mach Box method are generated by 
subroutine MBREG, so, in general, fixed dimensions are used. One additional buffer is used. 


Once the pointers, common MB0XN and common MB0XL are set up, MBAMG calls M8GE0D to set up the 
regions and geometry for the problem in common MB0XA. Than MBAMG calls MBREG to generate the boxes. 
MBREG fills in the box arrays based on the area to cover, Mach number, and number of boxes 
requested. MBCTR1 and MBCTR2 can be called to make box arrays if the control surface exists. 

M3PL0T is then called to print a picture of the planform regions. 


MBAMG then calls MBM0DE to generate mode-like data on SCR2. The SSPLIN routine is used to 
spline from the Mach box points to the input control points for the wing, control one and control 
two separately. 
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The Influence coefficients are computed by MBCAP, then SKJ (identity) is output, and finally 
MBDHDH is called to compute and output the AJJL contribution. 

Section two is a call to STPPT2 with outputs DlJK (Identity) and D2JK (Null). 

7, 2.7.4 Strip Theory Method 

Section one of the Strip Theory Method is driven by Subroutine STPDA. STPDA reads the ACPT, 
fills incormon STRIPL, and sets up pointers to common STRIPX where the various arrays will be 
stored. After all the input arrays have been set up an SKJ (Identity) matrix isbuilt. 

STPDA then calls: 5TPBG to build a BM and GM matrix for each strip; STPPHI to build the PHI 
functions for each strip; and finally STPA1C to combine these matrices and build AJJL. 

Section two is a call to STPPT2 which output DlJK (Identity) and D2JK (Null). 

3. 2.. 7.5 Piston Theory Method 

Section one of the Piston Theory method is driven by subroutine PSTAMG. PSTAMG reads the 
ACPI and sets up the core pointer to the arrays. Then SKJ (Identity) is output and PSTA is 
called to build AJJL. 

Section two is a call to STPPT2 with outputs DlJK (Identity) and D2JK (Null), 
j , 2 . 7 . $ Compressor Blade Method 

The flow for Section one of the compressor blade method is as follows. 
Subroutine AMGB1 is the driver for this method. It reads in the ACPT record for 
this method and locates reference parameters from the reference streamline on the 
blade. If there is enough core available, it calls AMGB1A to output one matrix 
of the AJJL list. When AMGB1A is through, AMGB1 bumps NRfSW and returns. Sutr»uT,u4. 

An&BtS is +0 coLunMi <s( 

Subroutine AMGB1A outputs a portion of the AJJL matrix for each streamline 
on the compressor blade. Each streamline may be subsonic, transonic or supersonic, 
depending on the Mach number for that streamline. Subroutine AMGB1B calculates 
terms for subsonic streamlines. Subroutine AMGB1C calculates terms for supersonic 
streamlines and subroutine AMGB1D calculates terms for transonic streamlines. 
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Each subnutrlx cf AJJL corresponds tc a blade streamline and Is of order 
NSTNS X NSTNS , where NSTNS Is the number of computing stations on the blade. The 
submatrices are located along the diagonal of AJJL. ■ 

The flow for Section two of the compressor blade method is as follows. 
Subroutine AMGB2 prepares all the computations necessary. It reads the ACPT 
record and locates the reference streamline parameters. Subroutine AMGB2A Is 
called to calculate ■ matrix [F'^J for each streamline. AMGB2 

outputs the NSTNS X NSTNS submatrix for each streamline to [D1JK],, 

Each submatrix of [SKJ] and [DIJtC] has the following form: 

[SKJ] = U • CF" 1 

and 

[D1JK] = [F-1] T 


The [D2JK] matrix is null. 

j?,3- .7.7 Swept Turboprop Blade Method 

The flow for Section one of the swept turboprop blade method is as follows: 
Subroutine AMGT1 is the driver for this method. It reads in the ACPT record for this 
method and locates reference parameters from the reference streamline on the blade. It 
calls AMGT1A to output one matrix of the AJJL list. Subroutine AMGT1S is called to 
output matrices used in SKJ. 

Subroutine AMGT1A outputs a portion of the AJJL matrix for each streamline on the 
swept turboprop blade. Each streamline must be subsonic based on the Mach number for 
that streamline. Subroutine AMGT1T calculates constants used by AMGT1B which computes 
the terms for the subsonic streamlines. AMGT1C and AMGT1D cause error exits with 
diagnostic printout when the streamlines are supersonic or transonic, respectively. 

Each submatrix of AJJL corresponds to a blade streamline and is of the order 
NSTNS *2 X NSTNS*2, where NSTNS is the number of computing stations on the blade. The 
submatrices are located along the diagonal of AJJL. 
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The flow for Section two of the swept turboprop blade method is as follows: 
Subroutine AMGT2 reads the ACPT record and locates the reference streamline parameters. 
Subroutine AMGT2A is called to calculate the matrix [F’^j for each streamline. AMGT2 
then outputs the transpose of this NSTNS X NSTNS submatrix twice along the diagonal 
for each streamline to [D1JK], Each submatrix of [SKJ] and [D1JK] has the following 
form: 


[SKJ] 


'i[f s rl)I; ’ 


[DIJKj = [SKJ] 


The [D2JKJ matrix is null. 


i .8 Subroutines 

Besides the module driver AMG, the subroutines of Section one ere divided into groups by 
method: 

i 

For the Doublet Lattice Methods five subroutines are shared: 

SNPDF, INCRfl), TKER, I0F1, and IDF2 

The Doublet Lattice Method without Bodies also uses: 

DLAMG, GEND, DPPS, and SUSP 

The Doublet Lattice Method with Bodies also uses: 

DLAMBY, SUBI, AMGBFS, FZY2, FWMW, BFSMAT, AMGR0D, AMGSBA, GENDSB, DPPSB, DPZY, DYPZ, 
OZPY, SUBB, SUBPB, DZY, FLLD, TV0R, DZYMAT, and P.jJWDZY 
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The Mach Box Method uses: 

MBAMG, MBPRIT, MBGEtfD, MBREG, MBCTR1 , MBCTR2, MBPL0T, MBM0DE, MBCAP, MBBSLJ, 06 , ZJ, 
MBDPOH, MBGAE , MBGAW , M6GATE , SUMPHI, and TRAILE 

The Strip Theory Method uses: 

STPDA , STPBG, STPPHI, STPAIC, STPK, STPBSO, and STPBS1 

The Piston Theory Method uses: 

PSTAMG and PSTA 

For Section two the subroutines are DLPT2, STPPT2 and DLBPT2. 

Subroutine Name: AMG 

1 . Entry Point: AMG 

2. Purpose: Module driver for AMG - see description above. 

3. Calling Sequence: CALL AMG 

3,j~.8.2 Subroutine Name: DLAMG 

1. Entry Point: DLAMG 

2. Purpose: Output AJJL and SJK parts for Doublet Lattice without bodies, 

3. Calling Sequence: CALL DLAMG (ACPI, AJJL, SKJ) 

ACPT - GIN0 file number of ACPT. 

AJJL - GIN0 file number of AJJL. 

SKJ - GIN0 file number of SKJ. 

4. Core Requirement: Core needed is four buffers plus record of ACPT plus 2*NJ, 

3 * }-.8.3 Subroutine Name: GEND 

1. Entry Point: GEND 

2. Purpose: Output all the columns of AJJL associated with a record on ACPT. 

3. Calling Sequence: CALL GEND ( NCARAY , NBARAY , YS , ZS , SG , CG , DT , W0RK , MAT0UT ) 

NCARAY, NBARAY, YS, ZS, SG, and CG are the locations of these arrays from the ACPT record. 
DT - location to put column of AJJL - complex 
W0RK - start of open core 
MAT0UT - GIN0 file number of AJJL. 
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31 -8.4 Subroutine Name: DPPS 

1. Entry Point: DPPS 

2. Purpose: Compute the elements in a row of AJJl. 

3. Calling Sequence: CALL DPPS (KS,I ,Jl ,J2,SGR,CGR,YS,ZS ,NBARAY,NCARAY,DT,W0RK) 

KS - Strip number in which receiver point I lies 

I - 3ox number of receiver point 

J1 - 1 

J2 - Number of boxes 

SGR - Sine of dihedral angle of receiver strip (from SG array) 

CGR - Cosine of dihedral angle of receiver strip (from CG array) 

YS, ZS, NBARAY, NCARAY - location of these arrays from ACPT record 
DT - location to start putting elements of column - complex 
W0RK - start of open core 

3.1 .8.5 Subroutine Name: SUBP 

1 . Entry Point: SUBP 

2. Purpose: Compute downwash factor element. 

3. Calling Sequence: CALL SUBP(I,L,LS,J,SGR,CGR,YREC,ZREC,SUM,XIC,DELX,EE,XLAM,SG,CG,YS,ZS) 

I - Box number of receiving point 

L - Panel number in which sending point J lies 

LS - Strip number in which sending point J lies 

J - Box number of sending point (also row number of output column) 

SGR - Sine (see DPPS) 

CGR - Cosine (see DPPS) 

YREC - YS(KS) - y coordinate from ACPT array 

ZREC - ZS(KS) - z coordinate from ACPT array 

SUM - Output element - complex 

XIC ,DELX ,EE ,XLAM,SG,CG,YS ,ZS - locations of these arrays for ACPT record. 

3.1 .8.6 Subroutine Name: SNPDF 

1 . Entry Point: SNPDF 

2. Purpose: Compute the steady downwash factors for one receiving-sending point combination. 
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3. 


Calling 

SL 

CL 

TL 

SGS - 
CGS - 
SGR - 
CGR - 
XO 
YO 
ZO 
ES 

DIJ - 
BETA - 
CV 


Sequence: CALL SNPDF (SL,CL,TL,SG$,CGS,SGR,CGR,XO,YO,ZO,ES,DIJ,BETA T CV) 

Sine of sweep angle of sending box 

:or.’ne of sweep angle of sending box 

Tanganc of sweep angle of sending box (from ACPT) 

Sine of dihedral angle of sending point 
Cosine of dihedral angle of sending point 
Sine of dihedral angle of receiving point 
Cosine of dihedral angle of receiving point 

X coordinate of receiving point, X coordinate of "center" of sending point 

Y coordinate of receiving point, Y coordinate of "center" of sending point 

Z cuordinate of receiving point, Z coordinate of "center" of sending point 

Sending point strip half width 

Steady contribution to downwash - output 

Square root of 1.0-M^ 

Chord of sending point 


3 , 1 . 8 . 7 Subroutine Name: INCR0 


1. Entry Point: INCR0 

2 . Purpose: Computes the unsteady downwash factor for one receiving-sending point combination 

3. Calling Sequence: CALL I NCRtir ( AX , AY ,AZ,AX 1 ,AY1 ,AZ1 , AX2,AY2,AZ2, SGR, CGR, SGS, CGS, KR,RL, BETA, 

SDELX.DELY.DELR.DELI) 


AX 

AY 

AZ 

AX1 

AY1 

AZ1 

AX2 

AYl' 

AZ2 

SGR 

CGR 

SGS 

CGS 


- XO 

- YO 

- ZO 

- X0+ES*TL 

- Y0+ES*CGS 

- Z0-ES*SGS 

- X0-ES*TL 

- Y0-ES*CGS 

- Z0-ES*SGS 


See definitions for SNPDF (Section 4.114.8.6). 


t 


5 
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4 . 


KR - Reduced frequency 
RL - REFC 

BETA - Square root of 1.0-M 2 

SOELX - Box chord of sending point 

DELY - 2.0* sending point strip half width 

DELR • - Output - real part of downwash factor 

DELI - Output - imaginary part of downwash factor 

Method: INCR0 calls TKER for three points for each receiving-sending box combination (the 

center, the inboard point, and the outboard point). Then INCR0 calls I0F1 and IDF2 to 
perform the integration of the kernels. 


?. > .8,8 Subroutine Name: TKER 

1. Entry Point: TKER 

2. Purpose: Compute incremental oscillating kernel 

3. Calling Sequence: CALL TKER(X,Y,Z,KR,BR,SGR,CGR,SGS ,CGS,T1 ,T2,M) 

X - AX, AX1 , or AX2 for center, inboard, outboard \ 

Y - AY, AY1 , or AY2 for center, inboard, outboard ( see INCR0 (see Section 4.114.8.7 

Z - AZ, AZ1 , or AZ2 for center, inboard, outboard ) 

KR - Reduced frequency 

BR - REFC/2.0 

SGR, CGR, SGS, CGS - See SNPDF (section 4.114.8.6) 

T1 - Output - cosine (v r - Y s )y - dihedral angle (receiving (r) or sending { s ) ) 

T2 - Output - [(Z cos Y r - Y sin y r ) x (Z cos Y s - Y sin y $ )]/(BR/M) 2 

M - Mach number 

4. Method: Kernel components are returned in common DLM. 

3 . 2 - ,8.9 Subroutine Name: IDF1 


1 . Entry Point: IDF1 

2. Purpose: Integration of the planar parts of the kernels 

3. Calling Sequence: CALL I DFT ( EE , E2 , ETA01 .ZET01 ,ARE,AIM,BRE,BIM,CRE,CIM,R1SQX,XIIJR,XIIJI) 

EE - Sending strip half width 

E2 - EE 2 
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ETAOl - AY cos Y r + AZ sin Y r (AY and AX, see I NCR0 ) ( y see TKER) (Sections 4.114.7, 4.114.8) 
ZETOl - AZ cos Y r - AV sin Y r 

ARE, AIM, BRE, BIM, vRE, CIM - coefficients of the parabola for planar part 
R1SQX •• AY 2 + A 7 ? 

XI I JR - output - real part of planar integral contribution 
XIIJI - output - Imaginary part of planar integral contribution 

3.2. .8.10 Subroutine Name: IDF2 

1. Entry Point: IDF2 

2. Purpose: Integration of the nonplanar parts of kernels 

3. Calling Sequence: CALL IDF2(EE, E2, ETAOl .ZETOl ,A2R,A2I,B2R,B2I .C2R.C2I ,R1SQX,DIIJR,DI IJI) 

EE, E2, ETAOl, ZETOl - same as IDF1 (Section 4.114.8.9) 

A2R, A2I, B2R, B2I , C2R, C2I - coefficients of the parabola for the nonplanar part 
R1SQX - See IOF'1 (Section 4.114.8.9) 

DIIJR - output - real part of nonplanar integral contribution 
DIIJI - output - imaginary part of nonplanar integral contribution 

3, .1.8. 11 Subroutine Name: DLPT2 

1. Entry Point: DLPT2 

2. Purpose: To output the Doublet Lattice without Bodies parts for matrices D1JK and D2JK. 

3. Calling Sequence: CALL DLPT2(ACPT,D1JK,D2JK) 

ACPT - GIN0 number 

D1JK - GIN0 number 
D2JK - GIN0 number 

j.jL.8-12 Subroutine Name: DLAMBY 

1 . Entry Pol nt: DLAMBY 

2. Purpose: Output AJJL and SKJ oarts for Doublet Lattice with Bodies 

3. Calling Sequence: CALL OLAMBY(ACPT ,AJJL ,SKJ) 

ACPT, AJJL, and SKJ are GIN0 file numbers 

4. Core Requirements: Four buffers plus record of ACPT plus 4*NJ. 
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j. 2.8.13 Subroutine Name: GENDS 8 

1. Entry Point: GENDSB 

2. Purpose: Generate part of the AJJL Influence coefficient matrix 

3. Calling Sequence: CALL GEND$Fi(NCARAY,NBARAY,SG,CG,NFL ,NBEA1 ,NBEA2,IFLA1 ,IFLA2,DT,0PY) 

NCARAY to IFLA2 - the locations of these arrays from ACPT record 

DT - storage for 2*NJ words 

OPY - storage for 2*NJ words 

4. Core Requirements: Up to 4 buffers may be used (2 for Y bodies, 1 for Z bodies, and 1 

for panels) . 

3. i .8.14 Subroutine Name: DPPSB 

1. Entry Point: DPPSB 

2. Purpose: Compute the element in a panel on panel row of AJJL. 

3. Calling Sequence: CALL DPPSB (KS, I ,J1 ,J2,SGR,CGR,YS,ZS,NBARAY, NCARAY, DT.W0RK) 

Same as DPPS 

J.i.8.15 Subroutine Name: DPZY 

1. Entry Point: DPZY 

2. Purpose: Compute the elements in an interference element on a panel in AJJL 

3. Calling Sequence: CALL DPZY(KB,IZ,I ,J1 ,J2,IFIRST,ILAST,Y8,ZB,AVR,ARB,TH1A,TH2A,NT121 , 

NT1 22 .NBA RAY .NCARAY ,NZYK 8 ,DPZ ,DPY) 

KB - Body number in which receiving point I lies 

IZ - Body element number of body KB in which I lies 

I - Receiving point 

J1 - Starting element number 

J2 - End'' 3 element number 

IFIRST- 0 ^ starting element 

ILA5T - 0-| ending element 

YB to NCARAY - locations of arrays in ACPT record 

NZYKB - Z-Y flag 

DPZ - Storage for row of AJJL 

DPY - Storage for row of AJJL 
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3 t X. .8.16 ' ^routine Name: DZPY 

1. Entry Point: DZPY 

2. Purpose: Compute the elements in a column of AJJL for Z interference elements 

3. Calling Sequence: CALL DZPY(KB,KS ,LS,I ,J1 ,J2 ,NYFLAG,SGR,CGR,FMACH,ARB,NBEA1 ,DT) 

KB - See DPZY (Section 4.114.8.15) 

KS - index of receiving point Y-Z coordinates 

LS - strip number 

I,J1,J2 - See DPZY (Section 4.114.8.15) 

NYFLAG- Type to build 

SGR.CGR - See DPPS (Section 4.114.3.4 

FMACH - Mach number 

ARB.NBEAl - location of arrays in ACPT record 
DT - Storage for "ow of AJJL 

3 . A .8.17 Subroutine Name: DYPZ 

1. Entry Point: DYPZ 

2. Purpose: Compute the elements in a column of AJJL for Y-interference elements 

3. Calling Sequence: CALL DYPZ(KB,KS ,LS, I ,J1 ,J2, NYFLAG, SGR.CGR, FMACH, ARB, NBEA1 .IBO.LSO, JBO.DT) 

KB to N8EA1 - See DZPY (Section 4.114.8.16) 

LBO - first body with Y orientation 

LSO - Z-Y coordinate index for fi 'St element of LBO 

JBO - Sending point index for first Y oriented body element 

DT - Storage for row of AJJL 

3iwL.8.18 Subroutine Name: SUBPB 

1. Entry Point: SUBPB 

2. Purpose: Compute downwash facto" elements on panels 

3. Calling Sequence: CALL SUBPB (I ,L,LS,J,SGR,CGR,YREC,ZR£C,SUM,XIC,DELX,EE,XLAM,SG,CG,YS,ZS, 

NAS , NASB , AVR , ZB . YB . ARB ,XLE , XTE , X , NB ) 

I to SUM - See SUBP (Section 4.114.8.5} 

XIC to X - location of arrays from ACPT record 
NB - number of bodies 
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3.1 .8.19 Subroutine Name: SU8B 

1. Entry Point: SU8B 

2. Purpose: Compute downwash factor elements on bodies 

3. Calling Sequence: CALL SUBB(KB,KS,I,J,JB,LB,LS,NDY,NYFL,PI,EPS,$GR,CGR,AR,BETA,SUM,RIA, 

DELX ,YB,ZB ,YS,ZS,X) 

KB - Index of receiving body 

KS - strip number of receiving point 

I - receiving point index 

J - sending point index 

JB - sending point index 

LB - body number of sending point 

NDY - Z-Y flag 

NYFL - type to build 

PI - TT 

EPS - .00001 

SGR - CGR - See DPPS (Section 4.114.8.4) 

AR - aspect ratio of body 

BETA - See SNPDF (Section 4.114.8.6) 

SUM - Output 

RIA-X - locations of arrays in ACPT record 
3 - Z .8.20 Subroutine Name: SUBI 

1 . Entry Point: SUBI 

2. Purpose: Compute the image point coordinates inside associated bodies on MU-Z and MU-Y. 

3. Calling Sequence: CALL SUBI(DA,DZB,DYB,DAR,DETA,OZETA,DCGAM,OSGAM,DEE,DXI ,TL,DETAI , 

DZETAI .DLGAMI .DSGAMI ,DEEI ,DTLAMI ,DMUY ,DMUZ, INFL , IOTFL) 

See Reference 1 (Section 4.114,11) for argument list. 

J? , .8.21 Subroutine Name: DZY 

1. Entry Point: DZY 

2. Purpose: Calculated effect of slender body element on a panel element 

3. Calling Sequence: CALL DZY(X,Y,Z,SGR,CGR,SI1 ,XI2, ETA, ZETA.AR.AO.KR, REFC, BETA, FMACH>LNS, 

IDZDY,DZDYR,DZDYI) 

See Reference 1 (Section 4.114.11) for argument list. 
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3 . X .8.22 Subroutine Name: TV0R 

1. Entry Point: TV0R 

2. Purpose: Calculate norinalwash at a point due to a trapezoidal unsteady Index ring 

3. Calling Sequence: CALL TV0R(SL1 ,CL1 ,TL1 ,SL2,CL2,TL2,$G$,CGS,SGR,CGR,X01 ,X02,YO,ZO»E, 

BETA,REFC,FMACH,KR,BRE,BIM) 

See Reference 1 (Section 4.114.11) for argument list. 

3 . 2 . .8,23 Subroutine Name: FLLD 

1. Entry Point: FLLD 

2. Purpose: Calculate the velocity normal to a surface due to a finite length line doubled. 

3. Calling Sequence: CALL FLLD(X01 ,X02,Y0,Z0,SGR,CGR,SGS,CGS,KR,REFC,FMACH,E t L,KDlR.KDl I , 

KD2R.KD2I) 

See Reference 1 (Section 4.114.11) for argument list. 

3,2 .8,24 Subroutine Name: AMGR0D 

1. Entry Point: AMGR0D 

2. Purpose: Calculate normalwash at panels and Interference elements due to slender elements 

3. Calling Sequence: CALL AMGR0O(D,BETA) 

D - storage for a row of AJJL 

2 

BETA - square root of l.-M 
3. ^.8. 25 Subroutine Name: DZYMAT 

1. Entry Point: DZYMAT 

2. Purpose: Calculate a slender element column of AJJL 

3. Calling Sequence: CALL DZYMAT(D,NFB,NLB,NTZYS,IDZDY,NTAPE ,X,BETA,IPRT,NS,NC,YS,ZS,SG,CG, 

YB,ZB,NBEA1 ) 

D - storage for a row of AJJL 

NFB - number of first body 

NLB - number of last body 

NTZYS - number of slender elements 

IDZDY - Z-Y flag 

NTAPE - GIN0 file for output 

X ■* location of array from ACPT record 
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BETA - see AMGR0D 
IPRT - print flag 

NS to NBEA1 - locations of array from ACPT record 
J.X .8.26 Subroutine Name: R0WDZY 

1. Entry Point: R0WDZY 

2. Purpose: To set up call to OZY 

3. Calling Sequence: CALL R0WDZY(NFB,NLB,R0W,NTZY5,D,DX,DY,DZ, BETA, IDZDY.NTAPE.SG.CG, IPRT, 

Y8,ZB,ARB,NSBEA,XIS1 ,XIS2,A0) 

NFB, NLB.NTZYS.D, BETA, IIDZD1 ,NTAPE, IPRT - same as DZYMAT 
R0W « row position of answer 
DX,DY,DZ - X, Y, Z of receiving point 

SG,CG ,YB,ZB ,ARB ,NSBEA ,XIS1 ,XIS2,A0 - locations of arrays from ACPT record 
.8.27 Subroutine Name: AMGSBA 

1. Entry Point: AMGSBA 

2. Purpose: Add slender body terms and pack out final AJJL for bodies 

3. Calling Sequence: CALL AMGSBA{ AJJL , A0 ,AR ,NSBE ,A) 

AJJL - GIN0 file number 

A0,AR - locations of arrays from ACPT record 
NSBE - number of slender body elements 
A - storage for a row of AJJL 

3,X .8-?8 Subroutine Name: AMGBFS 

1. Entry Point: AMGBFS 

2. Purpose: Build the SKJ matrix for bodies 

3. Calling Sequence: CALL AMGBFS (SKJ ,EE ,DELX .NCARAY .NBARAY ,XIS2 ,XIS1 ,A0 ,A0P ,NSBE) 

SKJ - GIN0 file number 

EE to A0P - locations of arrays from ACPT record 
NSBE - number of slender body elements 
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3.1 .8.29 Subroutine Name: BFSMAT 

1. Entry Point: BFSMAT 

2. Purpose: Form force matrices for slender elements 

3. Calling Sequence: CALL BFSMAT(ND,NE,NB, NP.NTP, LENGTH, NT0.SCR1 ,JF,JL, NAS.FMACH.YB, ZB, YS, 

ZS,X,DELX,EE,XIC,SG,CG,AR,RIA,NBEA1 ,NBEA2,NASB,NSARAY, 
NCARAY,BFS,AVR,REFC,A0,XIS1 ,XlS2,KR,NSBEA,NTlJ) 

See Reference 1 (Section 4.114.11) for argument list. 

5,2, ,8.30 Subroutine Name: FWMW 

1. Entry Point: FWMW 

2. Purpose: Add in images, symmetry, plane and ground effects 

3. Calling Sequence: CALL FWMW(ND,NE,SGS,CGS,IRB,A(J,ARB,XBLE,XBTE,YB,ZB,X5,YS,ZS,NAS,NAS3, 

KR,B£TA2 ,REFC,AVR,FWZ,FWY) 

See Reference 1 (Section 4.114.11) for argument list. 

3.2 .8.31 Subroutine Name: FZY2 

1. Entry Point: FZY2 

2. Purpose: Calculate the force numbers for FWMW 

3. Calling Sequence: CALL FZY2(XIJ,X1 ,X2,ETA,ZETA,YB,ZB,A,BETA2,REFC,KR,FZZR,FZZI ,FZYR, 

FZYI ,FYZR,FYZI ,FYYR,FYYI ) 

See Reference 1 (Section 4.114.11) for argument list. 

V.2-. 8.32 Subroutine Name: DLBPT2 

1. Entry Point: DLBPT2 

2. Purpose: Output the Doublet Lattice with Bodies parts of D1JK, D2JK 

3. Calling Sequence: CALL DLBPT2(ACPT,D1JK,D2JK) 

ACPT, D1JK, D2JK • GIN0 file numbers 

8.33 Subroutine Name: MBAMG 

1. En::ry Point: MBAMG 

2. Purpose: Driver for Mach Box Method 

3. Calling Sequence: CALL MBAMG ( ACPT ,AJJL,SKJ) 

ACPT, AJJL , SKJ - GIN0 file numbers 
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3. i ,8.34 Subroutine Name: MBPRIT 

1. Entry Point: MBPRIT 

2. Purpose: Print geometry data 

3. Calling Sequence: CALL MBPRIT(AW,AC,AT) 

AW - area of wing 

AC - area of control one 

AT - area of control two 

3 , 2 . .8.35 Subroutine Name: MBGE0D 

1. Entry Point: MBGE0D 

2. Purpose: Compute the geometry of the planform 

3. Calling Sequence: MBGE0O 

.i..2..8.36 Subroutine Name: MBREG 

1. Entry Point: MBREG 

2. Purpose: Compute the limits of the region and the percentage of box in each 

3. Calling Sequence: CALL MBREG( I REG , NW1 ,NWJ ,NC21 ,NC2N , NC 1 ,NCN,ND1 ,NDN,XK,YK,XK1 ,YK1 ,XK2, 

YK2 ,XWTE , YWTE ,KTE ,KTE1 ,KTE2 .PAREA ) 

IREF - flag for MBREG success - 2 = fail 

NW1 - PAREA - location of arrays which MBREG is to build 

1,1 .8.37 Subroutine Name: MBCTR1 

1. Entry Point: MBCTR1 

2. Purpose: Compute the region calculations for control one 

3. Calling Sequence: CALL MBCTR1 { ICl ,IRl .NCN.NC1 .NWN.NW1 .PAREA) 

IC1 - starting box number for control one 

IR1 - ending box number for control one 

NCN to PAREA - locations of arrays which M8CTR1 is to build 

3 , 1 - 8.38 Subroutine Name: MBCTR2 


1. Entry Point: MBCTR2 

2. Purpose: Compute the region calculations for control two 
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3. Calling Sequence: CALL MBCTR2( IL2,IR2,NC2N,NC21 .NWN.NWl ,PAREA) 

IL2 - startini box for control two 
1R2 - ending box for control two 

NC2N to PAREA - location of arrays which M8CTR2 Is to build 
3.‘i.8.39 Subroutine Name: MBPL0T 

1. Entry Point: MBPL0T 

2. Purpose: Print a representation of the planform 

3. Calling Sequence: CALL MBPL0T (NW1 ,ND1 ,NWN,NC21 ,NC2N,NC1 ,NCN,NDN) 

NW1 - NON - locations of arrays which define planform boxes 

.8.40 Subroutine Name: MBM0DE 

1. Entry Point: MDM0DE 

2. Purpose: Build the mode-like data from surface interpolation 

3. Calling Sequence: CALL MBM0DE(ACPT,SCR2,IC0R,NC|&R,Z,NI ,ND,XD,YD,1S,CR) 

ACPT, SCR 2 - GIN0 file numbers 

IC0R - first available location in MBAMGX 

NC0R - last available location in MBAMGX 

Z - start of open core 

NI - number of independent points 

ND - number of dependent points 

XD - X location of dependent points 

YD - Y location of dependent points 

IS - singularity flag 

CR - non-dimensionalizing number 

3 .8,41 Subroutine Name: M8CAP 

1. Entry Point: M8CAP 

2. Purpose: Compute the velocity potential Influence coefficients 

3. Calling Sequence: CALL MBCAP(NPNI,CAPPNI) 

NPNI - number of coefficients computed 
CAPPNI - location to store coefficients 
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J.JL *8.42 Subroutine Name: MBBSLJ 

1. Entry Point: MBBSLJ 

2. Purpose: Compute even-ordered Bessel Functions 

3. Calling Sequence: CALL MBBSLJ (ARG.N.BSL) 

ARG - input argument 
N - order 

BSL - storage for answers {length N) 
J.^.8.43 Subroutine Name: ZJ 



*•* 


E 


1 . Entry Point: ZJ 

2. Zero order Bessel Function 

3. Calling Sequence: X=ZJ(ARG) 

ARG - input argument 

3 . X -8.44 Subroutine Name: G 0 

1. Entry Point: G0 

2. Purpose: Evaluate an Expression n u ) - ip{$7, n^)] 

3. Calling Sequence: ANS=G0{R,ETAR,ETAL,EKM) 

R = ¥ 

ETAR = n u 
ETAL = n^ 

EKM = SI 

3, il. 8.45 Subroutine Name: M8DPDH 

1. Entry Point: MBDPDH 

2. Purpose: Driver for computing and outputing the terms of AJJL for Mach Box Method 

3. Calling Sequence: CALL MBDPDH {AJJL,F,DF, FI ,DF1 ,F2 ,DF2 ,XWTE,YWTE,PAREA,CAPPNI .DPNITE, 

DSS.Q.Ql ,Q2,NDN,ND1 ,NW1 ,NWN,KTE,KTE1 ,KTE2,NTE,NNC8, 
NNSBD, IW1 7 , IBUF,A) 

AJJL - GIN0 file number 
F - NTE - locations of array for Mach box 
NNCB - number or chordwise boxes 


I ? jE i! "•* ' T «— * * 
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NNSBD - number of spanwise boxes 
IW17 - GIN0 file number 
IBUF - pointer to a buffer 
A - storage for a row of AJJL 

], % .8.46 Subroutine Name: MbCAE 

1. Entry Point: MBGAE 

2. Purpose: Final calculation and output for AJJL 

3. Calling Sequence: CALL MBGAE (AJJL ,IN17,A ,F ,DF,F1 ,DF1 ,F2,DF2 ,Q ,Q1 ,Q2 >M00D) 

AJJL.IN17 - G1N0 file numbers 

F to Q2 - locations of Mach box arrays 
M00D - row number' of AJJL 

3 . 2- .8.47 Subroutine Name: MBGATE 

1. Entry Point: M8GATE 

2. Purpose: Compute sum on trailing edge 

3. Calling Sequence: CALL MBGATE ( NT0TE , DPH I TE , N , YWTE ,Q , Q1 , Q2 , KTE , KTE1 , KTE2 ) 

NT0TE - number of trailing edge terms 

OPHITE - KTE2 - locations of Mach Box arrays 

3. JL-.8.48 Subroutine Name: MBGAW 

1. Entry Point: MBGAW 

2. Purpose: Compute sum on wing 

3. Calling Sequence: CALL MBGAW(B0XL ,DPHI ,WS,PAW,PAF1 ,PAF2,Q, Q1 ,Q2,M,KC ,KC1 ,KC2) 

B0XL - box length 

DPHI - Q2 - location of Mach Box arrays 
M - KC2 - indexes to arrays 

p.i_.8.49 Complex Function Name: SUMPHI 

1. Entry Point: SUMPHI 

2. Purpose: Compute sum of (N*£K) on the wing 

3. Calling Sequence: SUM=SUMPHI(IXR,IYR,NDl ,NDN,CAPPNI ,DSS,N,M,ASYM) 


IXR, IYR,N,M,ASYM - index and flags 
ND1 ,NDN,CAPPNI ,D5S - location of arrays 
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3. X- 8. SO Complex Function Name: TRAILE 


1 . Entry Point: TRAILE 

2. Purpose: Compute sum of (N*AH) on tip 

3. Calling Sequence: SUM=TRAILE(X,J,H,P,M,B0XL) 

J,M,N - pointers 

X - values 

P - location of array 
B0XL - box length 

3 . X .8.51 Subroutine Name: STPDA 

1. Entry Point: STPDA 

2. Purpose: Driver for Section one of Strip Theory 

3. Calling Sequence: CALL STPDA(ACPT,AJJL,SKJ) 

ACPT.AJJL ,SKJ - GIM0 file numbers 

5.X .8.52 Subroutine Name: STPBG 

1 . Entry Point: STPBG 

2. Purpose: Builds two Intermediate matrices for Strip Theory calculations (BM and GM) 

3. Calling Sequence: CALL STPBG(BM,GM,NS,BL0C,D,CA,NSIZE) 

BM - storage for BM matrix 

GM - storage for GM matrix 

NS - number of strips 

BL0C - array of semi-chord lengths for strips 
D - array of hinge line lengths 

CA - array of control surface chords 
NSIZE - array of strip types 


J.A.8.53 Subroutine Name: STPPHI 

1. Entry Point: STPPHI 

2. Purpose: Calculate the functions 

3. Calling Sequence: CALL STPPHI(CA,BL0C,PM,NS) 

CA,BL0C,NS - See STPBG 

DM - Storage for <j> functions 
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3. 2 . .8.54 Subroutine Name: STPAIC 

1. Entry Point: STPAIC 

2. Purpose: Calculate and output AJJL for Strip Theory 

3. Calling Sequence: CALL STPAIC(BL0C,OY,NS1ZE,GAP ) BM,GM,PM,(1S,CLA,AJJL) 

BL0C,NSIZE,NS,BM,GM - See STPBG 

GAP - array of control surface gap 
PM - See STPPHI 

DY - array of Strip widths 

CLA - array of lift curve slopes 
AJJL - GIN0 flie numbers 

3, X. 8. 5 5 Subroutine Name: STPK 

1. Entry Point: STPK 

2. Purpose: Calculate the K-matrix for Strip Theory 

3. Calling Sequence: CALL STPK(£K,N,NST0P,N0PEN,NSTED,TSR,PM,CR,CI ,IM,J1 ) 
EK - modified reduced frequency 

N - strip number 

NST0P - strip type 

N0PEN - control surface flag 

NSTED - reduced frequency flag 

TSR - . 5*GAP/BL0C 

PM - 4 

CR - Theodorsen Function 

Cl - 0. 

IM - k-size 

Jl - J -size 

3.A..8.56 Subroutine Name: STPBSO 

1. Entry Point: STPBSO 

2. Purpose: J and Y Bessell functions of order zero 

3. Calling Sequence: CALL STPBSO(X,NC0DE,BJO,BYO) 

X - input argument 
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NC0DE - flag 
BJO - J Bessel 

BYO - Y Bessel 

3.JL.8.57 Subroutine Name: STPBS1 

1. Entry Point: STPBS1 

2. Purpose: J and Y Bessel Function of first order 

3. Calling Sequence: CALL STPBS1 {X,NC0DE,BJ1 ,BY1 ) 

X - input argument 

NC0DE - flag 
BJ1 - J Bessel 

BY1 - Y Bessel 

3, X .8.58 Subroutine Name: STPPT2 

1. Entry Point: STPPT2 

2. Purpose: Output D1JK and D2JK for Strip Theory, Mach Box and Piston Theory 

3. Calling Sequence: CALL STPPT2 ( ACPT , D1 JK,D2JK) 

ACPT.Dl JK.D2JK - GIN0 file numbers 

3. % .8.59 Subroutine Name: PSTAMG 

1. Entry Point: PSTAMG 

2. Purpose: Driver for Section one of Piston Theory 

3. Calling Sequence: CALL PSTAMG(ACPT,AJJL,SKO) 

ACPT,AJJL ,SKJ - GIN0 file numbers 

3,2 .8.60 Subroutine Name: PSTA 

1. Entry Point: PSTA 

2. Purpose: Calculate and output AJJL for Piston Theory 

3. Calling Sequence: CALL PSTA{DELTY ,BI ,CA,ALPH,THI ,AJJL) 

DELTY - array of strip width 

BI - array of semi-chord lengths for strips 

CA - array of chord lengths of each strip 
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ALPH - Alpha array (angle of attack) 

THI - Theta array (thickness ratio) 

AJJL - G1N0 file number 

j.JL .8.6/ Subroutine Name: AMG81 

1. Entry Point: AMG81 

2. Purpose: Driver for the compressor blade method tor A^L o.uA 

3. Calling Sequence: CALL AMGB1 ( INPUT , MAT0UT , ) 

INPUT = GI N0 file number for ACPT 

HAT0UT = GIN0 file number for AJJL 
SKJ = G1N0 file number for SKJ 
■j . X .8.£l Subroutine Name: AMGB1A 

1 . Entry Point: AMGB1A 

2. Purpose: Output all the columns of AJJL associated with A record 

of ACPT. 
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3. Calling Sequence: CALL AMGB1A (INPUT, MAT0UT , AJJ, AJJT, T20NX , 

TAMACH, TREFD) 

INPUT » GIN0 file number of ACPT 

MAT0UT o GIN0 file number of AJJL 

AJJ ° Storage for AJJL submatrices - complex 

AJJT = Storage for one column of AJJL 

T?0NX = Stores position of transonic submatrix In AJJL for a particular 
transonic streamline 

TAMACH = Stores Mach numbers of transonic streamlines 
TREFD = Stores reduced frequencies of transonic streamlines 

3 , i. .8.63 Subroutine Name: AMGB18 

1. Entry Point: AMGB1B 

2. Purpose: Calculates AJJL terms for subsonic streamlines. 

3. Calling Sequence: CALL AMGB1B (AJJL) 

AJJL ■» Location to put subsonic AJJL submatrix for this streamline 
J. 2 .B. 6 Y Subroutine Name: AMGB1C 

1 . Entry Point: AMGB1C 

2. Purpose: Calculates AJJL terms for supersonic streamlines. 

3. Calling Sequence: CALL AMGB1C (AJJL) 

AJJL n Location to put supersonic AJJL submatrix for this streamline 
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3 .^ . 8 .^ Subroutine Name: AMGB1D 

1. Entry Point: AMGB1D 

2. Purpose: Calculates AJJL terms for transonic streamlines. 

3. Calling Sequence: CALL AMGBlD (AJJL. T?0NX, TAMACH, TREDF) 

AJJL = AJJL submatrices for all subsonic and supersonic streamlines. 

It also contains space for transonic submatrices. 

T?0NX = (integer) - vector - non-zero Indicates transonic streamline 
zero if known streamline 
TAMACH « Vector of streamline Mach numbers 
TREDF = Vector of streamline reduced frequencies 

3.X • 8 . 6 fc Subroutine Name: INTERT 

1. Entry Point: INTERT 

2. Purpose: To linearly Interpolate by Mach number a transonic general 

Air Force matrix given two known streamline matrices. 

3. Calling Sequence: CALL INTERT (NL. NL1 , NL2 , NM, AJJ, TA) 

NL = Streamline number of unknown transonic 

NL1, NL2 = Two known streamlines 

NM « Size of matrix in AJJ = 2 * NSTNS * NSTNS 

AJJ ° Contains all generalized Air Force matrices for all 

s treaml ines 

TA = Vector of streamline Mach numbers 

Subroutine Names: SU 8 A, SUBBS^SUBC, SUBP, ALAMDA , AKP2. AKAPPA, 

DLKAPM, ASYCON, AKAPM, DRKAPM 

1. Entry Points: Tlu* &am«t as name 

2. Purpose: Called by AMG81C 
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3-2~ .8.4? Subroutine Name: GAUSS 

1. Entry Point: GAUSS 

2. Purpose: Equation Solver used by AMGB1B. 

3. Calling Sequence: CALL GAUSS (A, N » NL) 

3.X. .B.C*) Subroutine Name: AMGB2 


1. 

Entry Point: AMGB2 


2. 

Purpose: To output the compressor blade 

parts for matrices 


and 02JK. 



3. 

Calling Sequence: 

CALL AMGB2 (INPUT. 

WTJK, W2JK) 


INPUT - GIN0 file 

number for ACTP 



WlJtC «• GIN0 file 

number for DlJK 



W2JK = GIN0 file 

number for D2JX 



J,i .8. To Subroutine Name: AMGB2A 

1. Entry Point: AMGB2A 

2. Purpose: Calculate matrix used In the generation of &JUK.. 


3. Calling Sequence: CALL AMGB2A (INPUT. FMAT , XYHB, INDEX ^ 


INPUT = GIN0 file number of ACPT 
FMAT » Location for f F“ ^ J matrix 

XY2B « Location for basic coordinates of nodes on streamline 
INDEX “ Work storage for INVERS 


1 
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'),! -8.71 Subroutine Name: AMGB3S 

1. Entry Point: AMGBlS 

2 . Purpose: Calculate [F"1] matrix and W factor used In the generation of 

5KJ . 

3. Calling Sequence: CALL AMGBlS (INPUT, FMAT , XYZB, INDEX, RAOII , UFACT, 

NLINE) 

INPUT = GIN0 file number of ACPT 
FMAT c Location for [F"l] matrix 

XYZB c Location for basic coordinates of nodes on streamline 
INDEX « Work storage for INVERS 
WFACT ° Factor for output 
NLINE c Number of streamlines 
RADII = Streamline radius 
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3,2 .8.72 Subroutine Name: AMGT1 

1. Entry Point: AMGT1 

2. Purpose: Driver for the swept turboprop blade method for AJJL and SKJ 

generation. 

3. Calling Sequence: CALL AMGT1 (INPUT, MAT0UT, SKJ) 

INPUT = GIN0 file number for ACPT 

MAT0UT = GIN0 file number for AJJL 

SKJ = GIN0 file number for SKJ 


l 

I 
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3 . 1.8.73 Subroutine Name: AMGT1A 

1. Entry Point: AGMT1A 

2. Purpose: Output all the columns of AJJL associated with a record of ACPT. 

3. Calling Sequence: CALL AMGT1A (INPUT, MAT0UT, AJJ, AJJT, T?0f)X, TAMACH, 

TREFD, N2TN?2). 

INPUT = GIN0 file number of ACPT 

MAT0UT = GIN0 file number of AJJL 

AJJ = Storage of AJJL submatrices - complex 

AJJT = Storage for one column of AJJL 

T20NX = ] 

TAMACH = k Not used. 

TREFD = J 

N0TN02 * 2* no. of stations on a streamline. 


K 

I 

1 
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.8.74 Subroutine Name: AMGT13 

1. Entry Point: AMGT1B 

2. Purpose: Calculates AJJL terms for subsonic streamlines. 

3. Calling Sequence: CALL AMGT1B (AJJL, NSTNS2, C3, C4) 

AJJL = Location to put subsonic AJJL submatrix for this streamline 
NSTNS2 = No. of stations on a streamline X 2 
C3, C4 = Input constants 

3 , X ..8.75 Subroutine Name: AMGT1C 

1. Entry Point: AMGT1C 

2. Purpose: Writes error message for supersonic streamlines for turboprop blades. 

3. Calling Sequence: CALL AMGT1C (AJJL, NSTNS2) 

3,2, .8.76 Subroutine Name: AMGT10 

1. Entry Point: AMGT1D 

2. Purpose: Writes error message for transonic streamlines for turboprop blades. 

3. Calling Sequence: CALL AMGT1D (AJJL, TSONX, TAMACH, TRDF, NSTNS’) 

3.X .8.77 Subroutine Name: AMGT1S 

1. Entry Point: AMGT1S 

-It 

2. Purpose: Calculate [F J matrix used in the generation of SKJ. 

3. Calling Sequence: CALL AMGT1S (INPUT, FMAT, XYZB, INDEX). 

INPUT = GINO file number of ACPT 
FMAT = Location for [F _1 ] matrix 

XYZB = Location for basic coordinates of nodes on streamline 
INDEX = Work storage for INVERS 
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S-X ,8.78 Subroutine Name: AMGT1T 

7. Entry Point: AMGT7T 

2. Purpose: Calculate constants C3, C4 for subspnlc streamline on swept turboprop 
blade. 

3. Calling Sequence: CALL AMGT1T (NLINES, NLINE, INPfJT, NSTNS , C3, C4) 

NLINES = Total no. of streamlines on blade 
NLINE = Streamline being considered 
INPUT « GINO file number for ACPT 
NSTNS = No. of stations on a streamline. 

C3, C4 = Output constants 

3.X -8.79 Subroutine Name: AMGT2 

1. Entry Point: AMGT2 

2. Purpose: To output the swept turboprop blade parts for matrices D1JK and D20K. 

3. Calling Sequence: CALL AMGT2 (INPUT, D1JK, D2JK) 

INPUT = GINO file number for ACTP 

D1JK = GINO file number for DUX 

D2JK = GINO file number for D2JK 
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7ri.8.80 Subroutine Name: AMGT2A 

1. Entry Point: AMGT2A 

2. Purpose: Calculate [F -1 J matrix used In the generation of D1JK. 

3. Calling Sequence: CALL AMGT2A (INPUT, FMAT, XYZB, INDEX) 

INPUT = GINO file number of ACPT 
FMAT = Location for [F’^J matrix 

XYZB » Location for basic coordinates of nodes on streamline 
INDEX = Work storage for INVERS 
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3-2 .9 Design Requirements 


For Section one, four buffers are allocated at the bottom of core. For Section two, three 
buffers are allocated at the bottom of core. Each method may have its own open core common block 
but they must not overlap these buffers. 

y.l .9.1 Common Blocks 


AMGMN - 

Doublet Lattice without Bodies Communication 

Words 

1-7 

MCB - Trailer for AJJL 


8 

NR0W - Last row number output for any method on AJJL 

9 

ND - Y-syrrmetry flag | 


10 

NE - Z-syranetry flag ) 

1 record of AER0 Data Blpck 

11 

, REFC - Reference card ] 


12 

FMACH - Mach number (M) ] 

1 

| Pairs from 2 record pf AER0 Data Block 

13 

RFK - Reduced frequency > 

14-20 

75KJ - Trailer for SKJ 


21 

1SK - Row number to start building on SKJ 

22 

NSK - Last row number output for any method on SKJ 

AMGP2 - 

Section Two Communication 


Words 

1-7 

TW1JK - trailer for 01 JK 


8-14 

TWZJK - trailer for D2JK 


DLC0M - 

Doublet Lattice without Bodies Comnunication 

Words 

1 

NP - number of panels 


2 

N5TRIP- number of strips 
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DLC0M (Cont'd.) 

Words 

3 

NTP - 

4 

F 

5 

NJJ - 

6 

NEXT - 

7 

LENGTH- 

8 

INC - 

9 

IN8 - 

10 

IYS - 

11 

IZS - 

12 

IEE - 

13 

ISG - 

14 

ICG - 

15 

IX IC - 

16 

IDELX - 

17 

IXLAM - ; 

18 

IDT - ■ 

OLM - 

Both Doublet 


number of boxes 


NJ (Input parameter) 


pointer to XIC array 


ords 

1 

2 

3 

4 

5 

6 

7 

8 


KIO - Planar part of steady contribution to the kernel 

K20 - Nonplanar part of steady contribution to the kernel 

K1RT1 - 


K1IT1 - 
K2RT2P- 
K2IT2P- 

K10T1 - K10*T1 
K20T2P- K20*T2 


Unsteady parts of modified kernel 


T1 and T2 are defined under TKER 
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DLAXX - Open Core for Doublet Lattice without Bodies 

f — 

Words 6 
to 

End of Record on ACPT 


Column of AJJL 
2*NJ 


Free 


4 Buffers 


SKJ 

AJJT 

SerT 

mr 


DLP2X - Open core for Section two 


Record of ACPT 


Free 


D2JK 

3 Buffers dTJK 

ME! 


KPS - doth Doublet Lattice Methods 


Words 

1 

2 

3 

4 

5 


IHD - 0 = total kernel, 1 = incremental part only 

KD1R - real part of 

KD1 1 - Imaginary part of k-| 

KD2R - real part of k£ 

KD2I - imaginary part of k. 


(Doublet lattice Arrays) 
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PLBDY - Doublet Lattice with Bodies Communication 
Words 

1-12 Words 2-13 of ACPT record 
13-51 pointers into DLBXX for arrays on ACPT 

52 EC0RE - end of core in DLBXX 

53 NEXT - next available location in DLBXX 

54-58 SCR1-SCR5 - G I N0 file numbers for scratch files 
59 NTBE - number of columns to add to AJJL 


DLBXX - Open core for Doublet Lattice with Bodies 



(Dout let Lattice Arrays) 


MB0XA - Mach Box Wing Definitions 


Words 



1-12 

X 

X locations of wing 

13-24 

Y 

Y locations of wing 

25-34 

TANG - 

Tangents of wing sweep angles 

35-44 

ANG - 

Sweep angles of wing 

45-54 

C0TANG- 

Cotangents of wing sweep angli 


FUNCTIONAL MODULE AMG (AERODYNAMIC MATRIX GENERATOR) 


MB0XC - Mach Box Communications 


Words 

1-9 Words 2-10 of ACPT record 

10-30 Intercommunication between Mach Box subroutines 
MBAMGX - Open core for Mach Box Method 



Mach Box Arrays 


STRIPC - Strip Theory Communications 
Words 


1 

NS 

- number of strips 

2 

BREF 

- reference chord/ 2,0 

3 

CLAM 

- cosine of sweep angle 

4 

FM 

- Mach number 

5 

NCIRC 

- Theadorsen function selection 

6 

NNCIRC 

- NCIRC+1 

7 

EKR 

- reduced frequency 

8 

9-12 

Not Used 
BB(u) 

- b's for approximate function 


13-16 BEl'A(u) - B's for approximate function 

17-48 EKM(u,u) complex - storage for STPK output (k matrix) 
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STRIPX - Strip Theory Open Core 


Strip Theory Arrays 


from ACP| 


FREE 


SKJ 

4 Buffers 

AJJL 

AERO 


P ST0NC - Piston Theory Communication 


Words 


1-9 Words 2-10 of ACPT record 


PST0NX - Piston Theory Open Core 


Piston Theory Arrays 
from ACPT 


FREE 


4 Buffers 


SKJ 

a3jl 

TEW 

5cpt 
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BAHG1L and BAMG2L - Common Blocks for Compressor Blade Method 
Mords: 


1 

IREF 

- 

Reference streamline number 

2 

HiNMAC 

- 

Parameter M1NMACH 

3 

MAXMAC 

» 

Parameter MAi'MACH 

4 

NUNES 

- 

Number of streamlines on blade 

5 

NSTNS 

- 

Number of stations on blade 

6 

REFSTG 

- 

Reference blade stagger angle 

7 

REFCRO 

- 

Reference blade chord 

8 

REFMAC 

- 

Reference Mach number 

9 

REFDEN 

- 

Reference density 

10 

REFVEL 

- 

Reference velocity 

11 

REFFL0 

- 

Reference flow angle 

12 

SLN 

- 

Streamline number 

13 

NSTNSX 

- 

Number of stations on streamline 

14 

STAGER 

- 

Blade stagger angle 

15 

CHORD 

- 

Blade chord 

16 

RADIUS 

- 

Radius of streamline 

17 

BSPACE 

- 

Blade spacing 

18 

MACH 

- 

Relative flow Mach number at blade leading edge 

19 

DEN 

- 

Gas density at blade leading edge 

20 

VEL 

- 

Relative flow velocuty at blade leading edge 

21 

FL0WA 


Relative flow angle at bladr leading edge 

22 

AHACH 

- 

Internal Mach number 

23 

REDF 

- 

Internal reduced frequency 

24 

BLSPC 

- 

Internal blade spacing 

25 

AMACHR 

- 

Internal reference Mach number 

26 

TS0NIC 

. 

Transonic indicator 
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ORIGINAL PAGS £3 
OF POOR QUALITY 


1 A MG XX - Open Core for C&r»ftrti*r 61 o-des 


•"}* 1 r - ,"■ I. 

Words 6 
te 

End of Record on ACRT 


F 

Column of AJJL 


2*NJ 




Free 


4 Buffers 

SKJ 

ATJL r ' 1 

MI® ' 
MFT — 


A Mg’2X - Open core for Section two 

f Record of ACPI 


Free 


D2JK 

dTjT 


3 Buffers 
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TAMG1L and TAMG2L - Common Blocks for Swept Turboprop Blade Method 
Words : 

1 IREF - Reference Streamline number 

2 MINMAC - Parameter MINMACH 

3 MAXMAC - Parameter MAXMACH 

4 NUNES - Number of streamlines on blade 

5 NSTNS - Number of stations on blai 

6 REFSTG - Reference blade stagger angle 

7 REFCRD - Reference blade chord 

8 REFMAC - Reference Mach number 

9 REFDEN - Reference density 

10 REFVEL - Reference velocity 

11 REFFL0 - Reference flow angle 

12 SLN - Streamline number 

13 NSTNSX - Number of stations on streamline 

14 STAGER - Blade stagger angle 

15 CH0RD - Blade chord 

16 DCBDZB - 3C/3Z 

17 BSPACE - Blade spacing 

18 MACH - Relative flow Mach number at blade leading edge 

19 DEN - Gas density at blade leading edge 

20 VEL - Relative flow velocity at blade leading edge 

21 SWEEP - Sweep angle of blade 

22 AMACH - Internal Mach number 

23 REDF - Internal reduced frequency 

24 BLSPC - Internal blade spacing 

25 AMACHR - Internal reference Mach number 

26 TS0NIC - Transonic indicator 
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TAMGXX - Open Core for Swept Turboprop Blades 



TAMG2X - Open Core for Section two 



3 , JL . 1 0 Diagnostic Messages 


System fatal messages 3001, 3002, 3003, 3007, 3008 and (10) 3061. User fatal 


messages 2264 and 2265. 
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J.A-.ll References 


Most of the equations and code for the Doublet Lattice Method were taken from 

(1) Glesing, O.P., Kalman, T.P., Rodden, W.P., "Application of the Doublet-Lattice 

Method and the Method of Images to Lifting-Surface/Body Interference," AFFDL-TR-71-5 
Part II, Vol . 1, April 1972. 


Most of the equations and code for the Mach Box Method were taken from 

(2) Donato, V.W., Huhn, C.R., Jr. , "Supersonic Unsteady Aerodynamics for Wings with 
Trailing Edge Control Surfaces and Folded Tips," AFFDL-TR-68-30, August 1968, 


Most of the equations and code for the Strip Theory Method were taken from 

(3) Albano, E., "Strip Theory Aerodynamic Influence Coefficients for Wings with Aero- 
dynamical ly Balanced Control Surfaces," Northrop Corporation, Moran 01 vision Report 
NOR 68-125, August 1968. 


Most of the equations and code for the Piston Theory Method were taken from 

(4) Rodden, W.P., Forkos, E.F., Malcom, H.A., and Kllszcwski, A.M., "Aerodynamic 

Influence Coefficients from Piston Theory; Analytical Development and Computational 
Procedure," Space Systems Division, United States Air Force Report No. TDR-169 
(3230-11 )TN-2, August 1962. 
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J)-2> FUNCTIONAL MODULE AMP (AERODYNAMIC MATRIX PROCESSOR) 
i ,3> .1 Entry Point : 

AMP 

3.3 .2 Purpose 

The purpose of this module Is to produce "modal" aerodynamic matrices. This requires the 
combination of matrices from four sources. 

1. The aerodynamic matrices for aerodynamic cells, produced by the Aerodynamic Matrix 
Generator (AMG) module. 

2. The Interpolation from the structure to the aerodynamic cells, produced by the Geometry 
Interpolator (GI) module. 

3. The modes of the structure, produced by the Real Eigenvalue Analysis (READ) module, and 
selected by GKAM. 

4. The matrix of downwashes due to "extra" points, which may be supplied by the user via 
module INPUTT2. These extra points In NASTRAN are used for control systems and other 
special effects. 

3, J .3 DMAP Calling Sequence 

AMP AJJL,SKJ,DUK,D2JK,GTKA,PHIDH,DUE,D2JE,US£TD,AER0/QHHL,QKHL,QHJL/V,N»N0UE/V,N,XQHHL/ 
V.N.GUSTAER0 $ 

3.3 .4 Input Data Blocks 

AJJL Aerodynamic Influence matrix list 

SKJ Integration matrix 

01 JK Real part of downwash matrix 

D2JK Complex part of downwash matrix 

GTKA Aerodynamic transformation matrix k-set to a-set 

PHIDH Transformation be*veen modal and physical coordinates 

D1JE Downwash factors due to extra points; real 

D2JE Downwash factors due to extra points; complex 

USETD Displacement sets definition - dynamics 
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AER0 Aerodynamic matrix generation data 
Notes: 

1. AJJL, SKJ. D1JK, D2JIC, GTKA, PHIDH, USETD, and AER0 may not be purged. 

2. D1JE and D2JE are used only If N0UE > 0. Even In this case they may be purged. 

3,3 .5 Output Data Blocks 

QHHL — Aerodynamic matrix list - h-set 

QKHL -- Aerodynamic transformation matrix between hand b sets 

QHJL — Aerodynamic transformation matrix between j and k sets 
Notes: 

1. QHHL, QKHL, and QHJL are matrix lists - one submatrix for each (m,k) pair. 

2. If QHHL, QKHL, and QHJL are present before the module begins (APPEND on restart) and 
XQHHL = • 1 , the old data needed Is read from these data blocks. 

3,3 -6 Parameters 

N0UE -- Integer, Input, no default. The number of extra points. 

XQHHL - Integer, Input/output, no default. If +1, the data found on appended data blocks 
must be discarded. If -1, It can be used. AMP sets XQHHL to -1 on exit. 

GUSTAER0 - Integer, Input, default = 0. If, and only If, GUSTAER0 < 0, AMP will compute QHJL. 

3,3 .7 Method 

There are several important features which must be kept In mind. 

1. In general, the Input and output matrices may depend upon the aerodynamic parameters k 
(reduced frequency) and m (Mach number). A set of matrices (called a list) are processed 
in one pass through the module. 

2. Special code will be Introduced for restart. This Is required for Doublet Lattice solu- 
tions where matrix solution time may be long. This will allow the addition (or deletion) 
of (m,k) pairs without redecomposlng the dcwnwash matrix. 

3. An output, Q^i relating aerodynamic pressures to modal coordinates may be required for 
use in a data reduction module. This output wiM not be used in Phase 1; hence It will 
be purged from the calling sequence. The matrix of generalized forces, may be 
purged. If only data reduction Is desired. 
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The flow chart Is shown In Figure 1. The basic loop Is to write out matrices for the list of 
(m,k) pairs found on the AER0 data block. The source of these matrices Is normally the input. 

In the case of a restart which Involves only changes In the (m,k) pclrs, special code Is provided 
to avoid recalculation; and the matrices are found on the output data block which Ip declared 
APPEND. This occurs when XQHHL p -1. 

3*3 .7.1 Subroutine AMPA 

Put the output m,k list In core. This list Is found In the second record of the AER0 data 
block. The Index I will be used to Index down this list o' pairs. IMAX Is the number of pairs. 

Check to see If the output data blocks QKFL, QHHL, and QHJL exist and are valid. If they do, 
then this is a restart. These must be copied onto scratch files. Their (m,k) lists are put In 
core. Build a scenario file which lists the (m,k) pair, the AJJL column associated with this (m,k) 
pair, and the corresponding QHHL column. 

7,3 i.7.2 Subroutine AMPB 

Calculate the Dj h matrices. The superscript (1) is for the real part and (2) Is for the 
Imaginary part. 


E* ai ] « Partition of [$ dh ] , 

1 

0) j 

1 . 
i 

K,,] - [OX,] 

[G^3 may be needed again later to calculate Q^ h . 


[Djjh-wj’hX,] 

(2) 

:#>] ■ n$>]X M ] 

(3) 

[O'”] - Werge 0»<’>] . 

(4) 

CD » Merge Cojf ) D<| > ] . 

(5) 

If the Input data blocks are purged, 05^*' is zero. Start a loop with I ° 0. 

Check the time 


left. 
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J,3 .7.3 Subroutine AMPC 

Calculate (or find) If It is needed. It will be needed If either (a) Is to be output. 


or (b) Q hh Is to be output and Is not found on the scratch file, 
output only when their output data blocks are purged. If can 

It from there; otherwise, it must be calculated. First, check to 
for the present k. If not, find it by 

The Q k() and are not to be 

be found on a scratch file, get 
see If D Jh (k) has been calculated 

[D jh ] • »<’>] ♦ 1 k[D<*>] 

. (6) 

and save for possible later use. Next, solve for Qj h> The algebra included here will be theory 
dependent. The header record of AJJL will specify aerodynamic groups (see Section 4.115.7,5). 
Retrieve the submatrix [Ajj] from AJJL. If there Is more than one group, 0j h must be unpacked 
Into row groups. For each group, solve for [Qj h ], then pack the groups. For Doublet Lattice 
method, and the Double Lattice method with slender bodies, 

^jh^group = ^ A jj^group^°jh^group 

(7) 

For other methods. 


^jh^group " ^jj-^jh^group 

(7a) 

3,3 .7.4 Subroutine AMPD 


Calculate (or find) [Q hh ] and [Q kh ] if they are needed. They will be needed unless the output 
data blocks are purged. If [Q^3 can be found on a scratch file, get it there, otherwise. It must 
be calculated. If it must be calculated [Q^] will t 2 available. To compute [Qj )h 3 

Where S k j is a matrix list for (m,k), 


= ^ S kj^Qjh^ • 

[Q kh 3 is copied onto QKHL. 

(8) 

tQ ih 3 = tG ki^ ( W * 

(9) 

Fq l 

[ w • " er3e [q ■ 

(10) 


where [q gh ] is zero. Note that this requires only an update of [Q^]'s trailer. 
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Check the time. If [Qj h ] and [Q^] were calculated (rather than found), then the time per 
calculation can be found. If the time per calculation Is known and it Is not enough (with a 10* 
margin), no more loops should be attempted. 



Figure 1. AMP Module Flow. 


Repeat the loop (subroutines AMPC and AMPD) for additional values of I until the job Is done. 









I 
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If GUSTAER0 > 0, the following equations are evaluated: 
Partition PHIOH (Real) 




t 

■ o 

1 

f-e— * ' 

d 

PHIOH 

^ | 

PHIAH 

' 


i 

i* * 


Multiply (TRANSP0SE) (Real) 


< — 1< — $ — h — ^ <5 — h — > 


T 

* • 

T A 

* 

T 

* 

a 

| 

GIKA 

o a 

1 

PHIAH 

j 

GKH 

1 

■ a 

1 

. , 

l 

> • 


Start loop on reduced (m,k) pairs (use all). 
Multiply (TRANSPOSE) (S Is complex) 


(ID 


• j — ) 


■ 

T 

• 

T 

■ 

SKJ(k) 

9 

GKH 

i 

— > j 

S(k) 



• 

i 

. 


(13) 


Partition into Groups - (1) = Doublet Lattice, (2) = non-Doublet Lattice 





? 

J 2 

i 


S(l) 

S(2) 


Solve each group Rj^ : 

a. Doublet Lattice group 


(14) 


Ajj(gn>up) R jh (l) = S(l) 
b. Non-Doublet Lattice group 

# Jh <2) ■ Ajj S(2) 


(15) 


(16) 
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Merge Results 


[ 

[ 

r 

r 

o 


V'> ' 


V 

• *" 
fv 

V ) _ 

' J 

J 


R jh 


Append 


QJHL 


V>] — *> [a HJL j 


Repeat the last five steps for (m,k) pairs. 


c 

0 

0 

K 

I 


1 
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3 , £ .7.5 Matrix List Data Blocks 


The matrix list data block Is a special data block used for NASTRAN aeroelastlc calculations. 
It Is used to store a series of matrices. The matrices In the list will depend qpon two para- 
meters. The format is similar to that jf a matrix. If there are NMK matrices, each with NR0W 
rows and NC0L columns, then it will be stored like a matrix with NR0W rows, and NMK times NC0L 
columns. The matrix for the first parameter pair Is stored in the first NC0L columns. The 
matrix for other parameter pairs Is then added on at the end, one at a time. 

A special header record Is written which contains the following Information: 


Word 


Value 


1-2 

3 

4 

5-(4+2NMK) 


The name 

NC0L, for the individual matrices 
NMK 

M(I ) , K{ I ) , 1=1, NMK 
Other information 


If a single matrix exists. It can be read as a normal NASTRAN matrix. It Is possible that the 
matrix list was not completed by the generating module. The number of columns (found In the 
trailer) divided by NC0L should be an Integer. This should be equal to NMK. If it is less than 
NMK, it is the actual number of matrices on the list. For the AJJL, there Is additional informa- 
tion in the heaoe> 


Word 


Value 


(6+2NMK) 

(7+2NMK) - 
(6+2NMK+3NGP) 


NGP, number of uncoupled aerodynamic groups 

KT(N) , NJ(N) , N=l, NK(N) to NGPT, the theory 
identifier and the number of U. degrees of 

freedom associated with this group. \ NJ = NC0L. 


The matrix AJJL might look like (1 1' the identifier for Doublet Lattice theory): 


Pi 
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NC0L 

° 7 

KT ( 1 ) ° 1 

NMK 

° 3 

KT(2) p 1 

NGP 

= 2 

NJ(1) » 5 



NJ(2) o 2 



The output QICHL might look like (for 3 modes): 




NC0L n 4 (number of modes 
+ e«tra points) 

NMK b 3 (the M(I), K(i) are 
found from the AER0 
data block) 


The output QHHL might look like: 




NC0L ° 4 
NMK n 3 
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3-3 .8 Subroutines 

Numerous utility subroutines are used by the functional phases as shown below. 


AMPA AMPB 

AMPC 

AMPD 

AMPS 

AMPF 

CYCT2B CALCV 

CYCT28 

CYCT2B 

QYCT2B 

CVCT2B 

SSG2B 

SSG2C 

SSG2B 

SSG2B 

SSG2B 

MERGED 

CFACTR 

5KPREC 

SSG2A 

CFACTR 

PARTN 

CFBS0R 


SKPREC 

CFBSJ3R 


FILSWI 



FILSWI 

3.3 .8.1 Subroutine Name: AMPA 

TRANP.1 

SSG2B 



SKPREC 

1 . ' Entry Point: AMPA 






2. Purpose: To provide a scenario for later phases and to prepare for use of the appended 
output files. 

3, Calling Sequence: CALL AMPA (AER0, QJHL. QHHL, AJJL, QHHL0, QJHL0, INDEX, IMAX, IANY) 

AER0, QJHL, QHHL, and AJJL are the GIN0 file numbers of their respective data blocks, 

QHHL0 and QJHL0 are the GIN0 file numbers of two scratch files to hold valid submatrices from 
QHHL and QJHL on restart. 

INDEX Is the GIN0 file number of the scenario data block. Its contents are as follows: 


Record No. 

Word 

Contents 

0 

1 

Header 

1 

1 

M column number 


2 

K column number 


3 

AJJL cotumn number 


4 

QHHL0 column number {0 Implies recompute) 


FUNCTIONAL MQDUME AMP {AERODYNAMIC MATRIX PROCESSOR) 


IMAX Is the total number of (m»k) pairs on output. 

I ANY Is the flag for necessity to compute at least 1 QJH or QHH (0 Implies compute — 1 
Implies all retrieved modes). 

4. Common Blocks 

/AMPC0M/NC0LJ »NSUB , XM , Xl< , AJ JC0L ,QHHC(2)L , NGP , NGPD ( 2 , 30 ) .MCBQHH { 7 ) .MCBQtU ( 7 ) . NC0LH , 
IDJH„MCBRJH(7) 

NC0LJ - Number of columns In a submatrix of AJJL 

NSUF* - Number of submatrices In AJJL 

XM - Current M value 

XK - Current K value 

AJJC0L - Current column number of AJJL 

QHHC0L - Current column number In QHHL0 (a zero value Implies compute a new QHH) 

NPG - Number of groups 

NGPD - Two words for each group - Theory ID (1 = D.L.) - Number of columns of AJJ 
belonging to this group 

MCBQHH - Matrix control block for qHHL 

MCBQKH - Matrix control block for QJHL 

NC0LH - Number of H points 

IDJH * Flag for change in k value In (m,k) pair 

MCBRJH - Matrix control block for QJHL 

J.3.8 . 2 Subroutine Name: AMPB 

1. Entry Point: AMPB 

2. Purpose: To compute GKI and the DJH1 and DJH2. 

3. Calling Sequence: CALL AMPB (PHIDH 1 GTKA,D1JK,D2JK,D1JE,D2JE,USETD,DJH1 „DJH2 t GKI ,SCR1 , 
SCR2.SCR3) 

All Inputs are GIN0 file numbers. 
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3.3 .8.3 Subroutine Name: AMPB1 

1. Entry Point: AMPB1 

2. Purpose: To build a partitioning vector which will add a given number of columns to 
another matrix. 

3. Calling Sequence: CALL AMPB1 {IPVECT.NC0L1 ,NC0L2) 


3 . 


B XiX s a tiU #ue*BMShaat aa&aasa 
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IPVECT - the GIN0 file number on which the partitioning matrix will be built. 

NC0L1 - the number of columns In first matrix. 

JC0L2 - the number of columns In second matrix (to add onto the first matrix). 

3 .} .8.4 Subroutine Name: AMPB2 

1. Entry Point: AMPB2 

2. Purpose: This routine Is a general driver for PARTN. 

3. Calling Sequence: CALL AMP82 (A e All .A12.A21 ,A22,RP,CP 0 N1 „N2) 

A, All s A12, A21 , A22, RP # and CP are the GIN0 file numbers of the matrices supplied to PARTN. 

i R p ] 

( ) Faii Aiil 

[A] = M 

( j |_A21 A22j 

If any partition Is not desired, set Its file name to zero. 

N1 and N2 are the number of rows of RP and CP respectively. These are used only If RP or 
CP ° 0. A., RP and CP must have matrix trailers. Trailers will be written on all existing 
outputs. 

3,3 .8.5 Subroutine Name: AMPC 

1. Entry Point: AMPC 

2. Purpose: To compute (or retrieve) QJH and to form QJHL (If not purged). 

3. Calling Sequence: CALL AMPC (DJH1 ,DJH2,DJH,AJJL„QJHL,QJH0 11 QJHUA,SCR1 ,SCR2„SCR3 1I SCR4, 
SCR5.SCR6) 

DJH1 ,DJH2 ,0JH , AJ JL and QJHL are the GIN0 file numbers (GFN) of their respective data blocks 
QJH0 Is the GFN of a data block containing old QJH's from restart 
QJHUA Is the GFN of a data block containing the current QJH 
SCR1 - SCR6 are the GIN0 file numbers of six scratch files 

1 -3 .8.6 Subroutine Name: AMPC1 

1. Entry Point: AMPC1 


PRECEDING EAGE BLANK NOT FILMED 
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2. Purpose: To copy columns of one open matrix to another matrix. 

3. Calling Sequence: CALL AMPC1 (INPUT .OUTPUT .NC0L.IZ.MCB) 

INPUT ° GIN0 file number of the Input matrix 
OUTPUT » GINJB file number of the output matrix 
NC0L ° the number of columns to copy 
IZ = open core 

MCB ° matrix control block for OUTPUT. 

/UNPAKX/ and /PACKX/ control AMPC1 . 

4. Design Requirements: Both matrices must be opened and properly positioned. No trailers 
are written. Both matrices are left open. 

3 , 3 > .8.7 Subroutine Name: AMPC2 

1. Entry Point: AMPC2 

2. Purpose: To copy each column of the INPUT file onto the bottom of each column of the 
OUTPUT file. 

3. Calling Sequence: CALL AMPC2 (INPUT, OUTPUT, SCR1) 

INPUT, OUTPUT, and SCR1 are the GIN0 file numbers of their respective data blocks. 

4. Method: On the first entry, INPUT and OUTPUT are switched. On subsequent entries 
OUTPUT and SCR1 are switched and read together, one column at a time to produce OUTPUT. 

3,3 -8.8 Subroutine Name: AMPD 

1. Entry Point: AMPD 

2. Purpose: To compute or retrieve QHH and to write QHHL. 

3. Calling Sequence: CALL AMPD (QJHUA ,QHH0 ,SKJ „GKI ,QHHL ,SCR 1 , SCR3 , SCR4 ) 

All Inputs are the GIN0 file numbers of their respective data blocks. 
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3,1 .8.9 Subroutine Name: AMPE 

1. Entry Point: AMPE 

2. Purpose: To compute GKH 

3. Calling Sequence: CALL AMPE( PHIDH ,GTKA ,GKH ,SCR1 ,SCR2,USETA) 

All Inputs are the GIN0 file numbers of their respective data blocks. 

4. Method: AMPE calls CALCV and SSG2A to partition PHIOH. It then calls SSG2B to compute GKH. 

3.3 .8.10 Subroutine Name: AMPF 

1. Entry Point: AMPF 

2. Purpose: To solve for QHJL 

3. Calling Sequence: CALL AMPF(SKJ, GKH, AJJL, QHJL, PLAN, IMAX.SCRl ,SCR2,SCR3 1 SCR4,SCR5,SCR6, 

SCR7,SCR8,SCR9,SCR10) 

SKH, GKH, AJJL, QHJL, PLAN, 5CR1-SCR10 are GIN0 file numbers of their respective data blocks. 
IMAX Is the number of m,k pairs. 

3.3 .9 Design Requirements 

1. AMP requires 14 scratch files. These files are used as follows: 
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NAME 

DATA BLOCK 

COMPUTED 8V 

USED BY 

SCR1 

Old QHHL (QHHD) 

AMPA 

A 0 0 

SCR2 

Old QKHL (QJH|) 

AMPA 

' AoC 

SCR3 

Index of wprk to do (INDEX) 

AMPA 

A. DRIVER 

SeR4 

DJHT 

AMPB 

BoC 

SCR5 

QJH2 

AMPB 

B,C 

SCAB 

6KI 

AMPB 

B,B 

SCR 7 

DJH 

AMPB 

CpC 

SCR8 

QJHUA 

AMPC 

C.0 

SCR9 

Scratch File 


8,C,0 

SCR TO 

Scratch File 


3 pC 0 P 

SCR 11 

Scratch File 


8„C,D 

SCR12 

Scratch File 


C„D 

SCR13 

Scratch File 


C 

SCRT4 

Scratch File 


C 


2. Open Core: 


ROUTINE 

OPEN CORE 

FUNCTION 

AMP 

AMPB2X 

Buffer 

AMPA 

AMPA IX 

See layout 

AMPB 

AMPB2X 

CALCV 

AMPB1 

AMPB1X 

Buffer 

AMPB2 

AMPB2X 

PARTN 

AMPC 

AMPC IX 

Buffer, CVCT2B, AMPC1 

AMPpZ 

AMPC IX 

Buffer 

AMPD 

AMPOIX 

Buffer, CVCT2B 

AMPE 

AMPEX 

Buffer, CVCT2B 

AMPF 

AMPFX 

Buffer, CVGT2§) , AMPC1 


PRECEDING EAGE BLANK NOT FILMED' 
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FUNCTIONAL MODULE AMP (AERODYNAMIC MATRIX PROCESSOR) 


Open core-AMPAlX Is laid out as follows: 


Contents Lenath Pointer 



3, The loop between AMPC and AMPD would require much overlaying. Thus, AMP currently Is a 
single overlay chain. 

5-3 .10 Diagnostic Messages 

The following messages may occur: 3045, 3001, 3002 f 3003, 3008, 3007 

3007 occurs when a theory Is used which AMP does not understand; 3045 occurs when Insufficient time 
remains to compute another m-k pair. 
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FUNCTIONAL MODULE FAS (FLUnER ANALVSIS - PHASE 2) 

Entry Point : FA2 
J„y .2 Purpose 

To cal lee t data for reduction and presentation for each loop through the configuration 
parameters. . 

3, y .3 DMAP Cal 11 no Sequence 

FA2 PHIH.CLAMA.FSAVE / PHIHL.CLAMAL.CASIYV JVG / V.N.JSTART / C.Y.VREM.O / C t V 9 PR11fT»YE3 $ 
■j. y .A Input Data Blochs 

PHIH = Complex eigenvectors - h set, modal formulations. 

CLAMA ° Complex eigenvalue output table. 

FSAVE - Flutter storage save table. 

Note: No Input data block may be purged. 

jf ,tf .5 Output Data Blocks 

PHIHL - Appended complex mode shapes - h set. 

CLAMAL » Appended complex eigenvalue output table. 

CASEYY - Appended case control data table. 

0V6 - Output aeroelastle curve requests (V-g or V-f). 

Notes: 

1. No output data block may be purged. 

2. All output data blocks are read (DMAP attribute APPEND) on subsequent calls (FL00P from 
FSAVE t 1 If the method Is K). 


3'ij .6 Parameters 

TSTART - Integer-lnput/output-no default value. On input TSTART is the CPU time at the 
start of the DMAP flutter loop. On output TSTART will be -1 If there 1§ in? 
sufficient time for another DMAP loop. 

VREF - Real -user Input; no default. V fl(J ^ will be scaled by VREF: 


V A n 
out 


V/V, 


ref 


PRINT - BCD-user 1nput°defauTt = YES.. If PRINT’ « N0. no flutter summary will be printed. 
For YES the wlngflutter summary will be printed. 

For YESB the blade summary will be printed. 
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2.*f .7 Method 

The primary purpose of module FA2 Is to gather data for reduction and presentation. The 
header record of PSAVE will contain the METHOD. The actions of FA2 are method’ dependant. 
hi ,7.1 K-rMetthod 

Th1$ module is near the end of a OMAP loop. Its output files PHIHL, CLAMAL, CASE VV 0V6 are 

appended for each entry. On the first pass, special code must be executed to Initiate the files. 

The complex eigenvalues X have been found by module CEAD. These should have been sorted by 
Im(X) Increasing, Only use the first "NVALUE" modes. The quantities that need to be computed a^ei 

V out ° M*)/v ref • 

((2.0) Re(X)/Im(X) If Im(X) t 0 
9 

( 0 If Im(X) » 0 , 

. f ■* k Im(X)/2Tib re ^. , 

V Mach 13 ^sound" 1 ^ V ref 

The values of the parameter FL00P, m, k, b pef and NVALUE are found In the file PSAVE. A printer 
output Is prepared. 

The PHIHL, CASEYY and CLAMAL data blocks are created by appending the PHIH . CASEVY and CLAMA 
data blocks. 

The CASEYY data block Is for modules SDR2 and PLOT. It must keep In step with the append vectors, 
m, k, p and FLO0P will be added to the LABEL. 

The 0VG data block Is appended each time through the L00P. This will be used to create V-g 
or V-f plots, m, k, p and FL00P will be added to thp LABEL. 

.7.2 PK-Method 

The values for X, I, and G are supplied by FA1 on FSAVE for all eigenvalues and al] MACH 
number RH0 pairs. FA2 collects all k values together and outputs each collection of N such eigen* 
values on a curve for V-g plotting. CASEYY, PHIHL and CLAMAL are not written. 

j.y .7.3 KE-Method 

Existing FSAVE records contain records of length 2 x N where 2 ° Real, Imag a V, N ° number 
of modes. 
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FUNCTIONAL MODULE FA2 (FLUTTEK ANALYSIS - PHASE 2) 


The records are sorted by (see Figure 1) 

m a mach number 
k ° reduced frequency 
p ° density. 

The output should be sorted by 

m o mach number 
p o density 
n = mode number. 

The records will be used for 0V6, and for formatted print. 


A special "sorting" algorithm will be used to order the roots. For the first k value 
In each loop, the roots are accepted In the order of ALLMAT. Define 

the 1th eigenvalue for the 
Pj n = nth reduced frequency k. 

In the above, 1 = 1, 2, 3, ... (number pf modes) 

n ° 1, 2, 3, ... (number of k values) 

Define the extrapolated value based upon previous n's to be 


P l!n ° P 1,(n-1) + (k n“ lt n-l ,(P 1 t (n-l)‘ p 1,(n-2) )/(k n-r l{ n-2 > 


n > 3 


where P^g will be chosen equal to ^ . Then, select for P^ n the root found In the nth loop 
closest to Pj® n . Delete that root and let P ? n be the one qf the remaining roots closest to 
P 2 6 n * C° nt * nue until all roots are exhausted. The measure of "closeness" of the complex num- 
bers Is the square of the magnitude of the difference. If P 1 and P g are two roots, 


P 1 a Re + 1 Im P 1 , 


then the square Is, 


rRe^Pg-P ])] 2 + Clm(P 2 »?i)3 a . 

All eigenvalues are put on the FSAVE data block to be passed to FA2 module. CASEVY, PHIHL and 
CLAMAL are not written. 
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K - method 
(FA2 In loop) 

Output in order received 

Point ® (m,k,p) triplet 

# of entries ° # of modes 

l( - method 
(no loop) 

Sorting required 

Point ° (m,p,mode) triplet 

0 of entries D 9 of k's 

PK - method 
(no loop) 

Transpose required 

Point ° (m,p,mode) triplet 

H of entries ° & of V's 


Come In loop and 
PHIH , CUMA slots 

(No change) 



Note : All must have 0VG Figure 1. 

J.V .8 Subroutines 

Utility routine CYCT2B Is called. 

3.^,9 Design Requirements 

Open core for FA2 Is at /FA2X/ . 

FA2 uses no scratch files. 

3,*} .10 Diagnostic Messages 

The following messages may occur: 3001, 3002, 3003, 3007, 3008 and 3045. Only 3p45 Is 
a user message. It Indicates that the DMAP loop was not completed by exhausting the '•onflgu 
ration parameters but rather by a tlme-to-go failure. 







FUNCTIONAL MODULE APDB (AERODYNAMIC POOL DISTRIBUTOR FOR BLADES) 


3-5 FUNCTIONAL MODULE APOB (AERODYNAMIC POOL DISTRIBUTOR FQR BLADES) 

3-S •) Entry Po1:it : APDB 

3.5 .2 Purpose 

Bulk data cards which control the solution of aerodynamic problems are 
processed and assembled Into various blocks For convenience and efficiency In 
the solution of the aerodynamic problem. APOB also generates the transfermatl on 
matrix [GkaJ T (GTI(A) and the partitioning vector PVECT. 

3iS *3 DMAP Calling Secueuce 

APDB EDT , USET, 8 GPDT , CSTM , EQEXIN, GM, G0/ AER0, ACPT, FLIST, G7KA, 

PVECT/ V, N, NK/ V, N, NJ/ V, Y, MINMACH/ V, Y. MAXMACH/ V, Y, JREF/ 
V, V, MTYPE/ V. N , NEIGV/ V. 1 , KINDEX v -1 $ 

3*5 '.4 Input Data Blocks 

EOT - Aerodynamic bulk data cards 

USET * Displacement set definition table 

BGPDT - Basic grid point definition table 

CSTM - Coordinate system transformation matrices 

EQEXIN - Equivalence between external points and scalar index values 

GM - Multipoint constraint transformation matrix 

G0 - Structural matrix partitioning transformation matrix 

Notes : 

1. EDT, USET, BGPDT and EQEXIN cannot be purged. 


2. CSTM may be purged if all points are in the basic system. 


FUNCTIONAL MOOULE APDB (AERODYNAMIC POOL DISTRIBUTOR FOR BLADES) 

3. GM and G(J may be purged IF there are no multipoint or no omitted 
points . 

2 * S .5 Output Data Blocks 

AER0 - Control Information for control of aerodynamic matrls generation 
and flutter analysis 

ACPT - Information pertaining to each Independent group of aerodynamic 
elements 


FL1ST • Contains 


FLFACT and FLUTTER cards copied from gpT 


GTKA - Aerodynamic transformation matrix - K set to a set 

PVECT - Cyclic modes partitioning vector fc matrix PHIA from module CVET2 

Notes : 

1. AER0, ACPT 0 FUST and GTKA cannot be purged. 

2. PVECT may be purged If there are no cyclic modes to be partitioned. 

3 .6 Parameters 

NK - Output - Integer - no default. Degrees of freedom In the NK 

displacement set. 

NJ - Output - Integer - no default. Degrees of freedom In the NJ 

displacement set. 


- Input - real - default <■ 0.8. This Is the maximum Mach. number 
below which the subsonic unsteady cascade theory Is valid. 

• Input - real - default 1.01. This is the minimum Macfe number 
above which the supersonic unsteady cascade theory Is valid. 
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FUNCTIONAL MODULE APOB (AERODYNAMIC POOL DISTRIBUTOR FOR BLADES) 


1REF ~ Input - Integer - default ° »1. This defines the streamline number 
of the reference stream surface. I RE F must equal an SEN Op 
a S7REAML2 card. The default value, -1, represents the stream 
surface at the blade tip. If IREF does not correspond to an S^N, 
then the default will be taken. 

MYTPE - Input - 8CD - default ° CPSINE. This controls which components of 
the cyclic modes are to be used in the modal formulation, fITVPE o 
SINE for sine components and MTVPE ° C0SJNE for cosing components, 

NEIGV - Input - BCD - no default. The number of eigenvalues found. 

Usually output by the READ module. 

KINDEX - Input - BCD - default ° «1. Harmonic index number used in cyclic 
analyses. 


3‘S Method 


Subroutine APOB is the main control program for this module, It allocates 
buffers, reads input files, and initialises output files. APOB creates the AER0, 

ACPT and FLIST tables and generates the PVECT partitioning vector. Subroutine APDBT or 
APDBg generates the GTKA transformation matrix, 
much like module SSG2, using the following matrix operations: 


They reduce [G* ] 




£-T 

6 kn 


gT 

KM 


' 


C«lJ ■ [«J T !“!J * t&D 


KN 

tej„3 - Ul 


KM 


KN 


Kf 


’KS 
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FUNCTIONAL MODULE APOB (AERODYNAMIC POOL DISTRIBUTOR FOR BLADES) 




4 


J K*S 




At each step where a matrix multiply Is Indicated, the multiply Is shipped If the 
result Is known to be zero (l.e., U n or Ug are null), 

3,J .8 Subroutines Called 

Utility routines BISL0C, CALCV, SSG2B, TRANSS and GMMAT? are called. 

3 ‘ Si. 8.1 Subroutine Name: AP0B1 


1. Entry Point: APDB1 

2. Purpose: To generate transformation tnatrlR [G^L for compressor blades 
Method 6. 

3. Calling Sequence: CALL APDB1 (IBUF1, IBUF2, NEXT, LEFT, NSTNS, NUNES, XgfGN, 
LC2TM, ACSTM, NODEX, NODE!, ISILC, XYZB), 

3.J” .8.2 Subroutine Name: APDB2 


1. Entry Point: APDB2 

2. Purpose: To generate transformation matrix [G^], for turboprop blades 
Method 7. 

3. Calling Sequence: CALL AP0B2 (IBUF1, IBUF2, NEXT, LEFT, NSTNS, NLINES, XSIGN, 
LCSTM, ACSTM, NODEX, NOOEI, ISILC, XYZB). 


3.5 .8.3 Subroutine Name: APDB2A 


1. Entry Point: AP0B2A 

2. Purpose: To generate basic to local transformation matrix for APD82. 

3. Calling Sequence: CALL APDB2A (NLINE, NL, 2CR1, NSTNS, Ml, SI. 2N, TBLT, TBLR). 
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3.5 .9 Design Requirements 

Open core Is located at /APDB23/. APOB uses five scratch files. 

3.5 .10 Diagnostic Messages 

System fatal messages 3001 , 3002# 3003, 3008 and 3037 may occur. The AP0B module 
also generates Its own messages that are not nyrnhered. These messages are self* 
explanatory. 
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SET1 

32 
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SETS 

SPLINE1 

32 

32 

©F POOR QUALITY 

SPLINES 

32 


MKAER01 

34 


MKAER02' 

34 


AIFACT 

35 


FLFACT 

36 


FLUTTER 

36 


AERO 

37 


CAER01 

37 


FMETHOOS 

38 


VREF 

39 


TF 

40 
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42 


CTYPE 

41 
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41 
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41 
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42 
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n 

, 
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42 
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42 
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42 
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58 
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59 
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60 
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File Name Restart Table 


File Name Bit Pos , 

File Name 

Bit Pos. 
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MODAL FLUTTER ANALYSIS OF AN ADVANCED TURBOPROPELLER 

A. Description 

The dynamic aeroelastic stability of a ten-bladed advanced turbopropeller 
at a given operating condition Is examined in three phases: 

Phase 1 generates the differential stiffness matrix when the propeller Is 
subjected to the centrifugal loads due to its rotation. The presence of any 
steady state airloads is neglected. 

Phase 2 calculates the "running" natural modes and frequencies using the 
total (elastic plus the differential) stiffness from, and the geometry at the 
end of, phase 1. This phase is checkpointed, and allows the user to select the 
structural modes to be included for flutter analysis. 

Phase 3 is a restart of Phase 2 and computes the flutter eigenvalues. V*g» 
V-f curves and plotted to examine stability. 

The blades of the propeller are assumed to be: 

a) identical in all respects, 

b) mounted on a relatively rigid hub, and hence 

c) structurally independent, and 

d) aerodynamical ly coupled via discrete interblade phase angles 

a = 2irn/10, n - 0, ±1, .... ±4, 5. 

Therefore, only one blade of the propeller is modelled as shown in Figure 1. 

B. Input 

1 . Parameters : 

Propeller 

Number of blades = 10 

Diameter at blade tip (grid point 4) = 23.8 in. 

Diameter at shank root (grid point 153) = 4.1 in. 

Chord! ength at tip = 1.28 in. 


4.1 


3.03 In. 


Chordlength at root 
Sweep angle at tip 
Sweep angle at root 
Young’s modulus 
Poisson's ratio 
Material density 
No structural damping Included in flutter calculations. 
Operating Point 


= 51 .0 deg. 

= -15.9 deg. 

■ 16.0 x IQ 6 lbf/1n 2 
= 0.35 

= 4.141 x 10" 4 lbf.sec 2 /1n 4 


Blade setting angle (at 8.98"R) with the 
plane of rotation 


= 69.0 deg. 
= 6800 rpm 
= 0.70 


= 9336 In/sec. 

« 9.763 x 10" 8 ",bf.sep 2 /1n 


Rotational speed 
Free stream Mach number 
Free stream velocity 
Free stream density 
2. Constraints: 

All degrees of freedom at the root of the shank are constrained to ?ero. 

C. Results 

Phase 1: Results are shown in Figures 2 and 3. 

Phase 2: Figures 4 through 6 illustrate the first three natural modes at 
6800 rpm. 

Phase 3: Typical flutter results are shown in Figures 7 and 8 wherein the 
first bending mode is seen to be unstable (g, u > 0). 

The results are in good agreement with the experimental observations, as 
reported i n Reference 1 . 

D. Reference 

1. Elchuri, V., and Smith, G. C. C., "NASTRAN Flutter Analysis of Advanced 
Turbopropel lers," Final Technical Report, NASA CR- 167926 , April 1982. 
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Figure 8. SR-5 FLUTTER ANALYSIS AT 6800 RPM ( 10 BLADES ) 

ROOT LOCUS OF FIRST BENDING MODE ( FREQ IN VACUUM = 187.9 HZ. 
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= 09 

GRID 

147 


- = 27 

2 o 44 

a 04 

GRID 

148 


oO 

2o44 

= 0 

GRID 

149 


o 27 

2o 44 

— = 04 

GRID 

150 


o47 

2o44 

— = 08 

GRID 

151 


-o47 

2 a 06 

= 09 

GRID 

152 


-o27 

2o06 

= 04 

GRID 

153 


oO 

2o 06 

= 0 

GRID 

154 


a 2 7 

2o 06 

-o04 

vjRID 

155 


= 47 

2a06 

— o09 

MAIL 

1 

1=6 

E7 

o 35 

=0004141 

P QUAD 2 

1 

1 

oO 12 



PQUAD2 

2 

1 

= 024 



PQUAD2 

3 

1 

o032 



PQUAD2 

4 

1 

o0 36 



PQUAD2 

5 

1 

= 030 




4.19 



I 

gp SORTED I U l K 0 A IT A ECHO 


0 1 

o o 2 

o o 

3 o o < 

PUUA02 

6 

1 

a 0 1 S 

PQUA02 

7 

1 

o014 

PQUA02 

e 

1 

o028 

P4JUAD2 

9 

1 

<>037 

PQUA02 

10 

1 

<=043 

PQUA02 

11 

1 

oQ36 

PQUAD2 

12 

1 

o022 

POUAD2 

13 

1 

o 0 16 

PQUAD2 

14 

1 

o032 

P QUAD 2 

15 

l 

o 04 8 

PUUAD2 

it> 

1 

o051 

PQUAQ2 

17 

1 

o042 

PQUAD2 

18 

1 

o023 

P3UAD2 

19 

1 

o018 

PQUAD2 

20 

1 

o034 

PQUAD2 

21 

1 

<,053 

P3UA02 

22 

1 

o058 

PQUA02 

23 

1 

0 046 

PQUA02 

24 

1 

o025 

PQUA02 

25 

1 

o021 

P3UA02 

26 

1 

o042 

P0UA02 

27 

1 

o 06 1 

PQUAD2 

28 

1 

a 066 

P3UAD2 

29 


o05 1 

PQUA02 

30 

1 

o027 

PQUA02 

31 

1 

<>024 

PQUA02 

32 

1 

a 049 

PQUAD2 

33 

1 

o070 

PQUAD2 

34 

1 

o07 3 

PQUAD2 

35 

i 

o057 

P3UAD2 

36 

1 

oC 30 

PQUA02 

37 

1 

o 028 

PQUAD2 

38 

1 

<.054 

PQUAD2 

39 

1 

o07d 

P0UAD2 

40 

1 

o082 

P3UAD2 

41 

1 

o065 

PQUAD2 

42 

1 

<=035 

PQUAD2 

43 

1 

<=031 

P0UA02 

44 

1 

o 06 1 

PQUA02 

45 

1 

o 08 8 

PQUA02 

46 

1 

<= 09 3 

PQUAD2 

47 

1 

o075 



4.20 


SQRFEO BULK D & T A ECHO 


I 


I 


I 


PQUAD2 
PQUAD2 
PQUAD2 
PQUA02 
PQUA02 
PQUAD2 
PQUAD2 
p P9UAD2 
PQUAD2 

* P9UA02 
PQUAD2 
P0UAO2 

* PQUA02 
PQUA02 

r PQUAQ2 
„ PQUAD2 
P QUAD 2 
- PQUAD2 
PQUA02 
41 PQUA02 
PQUAD2 
F PQUAD2 
t» PQUAD2 
P0UAO2 
" PQUAD2 
L P QUAD 2 
PQUAD2 
_ PQUAD2 
T P9UAD2 
PQUAD2 
PQUAD2 
t PQUAD2 
i PQUAD2 
PQUA02 

E PQUAD2 
PQUAD2 
PQUAD2 

E PWUAD2 
PQUAD2 
PQUAD2 
PQUAD2 
r PQUAD2 


o 2 

o o 

3 op ^ 

48 

1 

oG>39 

49 

1 

oQ38 

50 

1 

o 06 8 

51 

1 

o 09© 

52 

1 

o 103 

53 

1 

o083 

54 

1 

o046 

55 

1 

o041 

56 

1 

oQ76 

57 

1 

ollO 

58 

1 

o 1 18 

59 

1 

o 09 1 

6G 

1 

o04 7 

61 

1 

o043 

62 

1 

o083 

63 

1 

o 120 

64 

1 

o 129 

65 

1 

= 100 

66 

1 

= 044 

67 

I 

= 045 

68 

1 

= 090 

69 

1 

o 135 

70 

1 

o 138 

71 

1 

= 100 

72 

1 

o 048 

73 

1 

o 053 

74 

1 

= 106 

75 

1 

,152 

76 

i 

= 148 

77 

1 

o G99 

78 

1 

o 044 

79 

1 

o063 

80 

1 

o 123 

81 

i 

= 171 

82 

1 

= 157 

83 

1 

o 099 

84 

1 

= 046 

85 

1 

oQ7 L 

86 

1 

o 14 1 

87 

1 

= 206 

86 

1 

o 1 77 

89 

1 

= 112 


6 


9 


10 
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SORTED BULK 


DATA ECHO 


, L 

o o 2 

o o 

3 o o 4 

PQUA02 

90 

1 

o 048 

PQUAD2 

91 

1 

, 0 84 

P QUAD 2 

92 

1 

o 172 

PGUAQ2 

93 

1 

o232 

PUUAD2 

54 

1 

o 198 

PQUAD2 

55 

1 

o 135 

PQUAD2 

96 

1 

o 062 

PQUA02 

97 

1 

o 119 

PQUAD2 

58 

1 

o2Q6 

PQUAD2 

99 

1 

o 266 

PQUAD2 

LOO 

1 

o 230 

PQUAD2 

101 

1 

a 152 

PQUAD2 

102 

1 

,071 

PQUAD2 

103 

1 

o 16 1 

PQUAD2 

1 04 

1 

,237 

PQUAQ2 

105 

1 

o 347 

P0UAD2 

106 

1 

o 31 9 

P3UAD2 

107 

i 

o 167 

P0UAD2 

108 

1 

,075 

P3UA02 

1C9 

1 

,222 

P3UAD2 

no 

1 

,373 

PQUAD2 

121 

1 

o 242 

PQUAU2 

122 

1 

oQ89 

PQUAD2 

123 

1 

,44 1 

P3UAD2 

124 

1 

,830 

PQUAD2 

125 

1 

o 830 

P QUAD 2 

126 

1 

,441 

PQUAD2 

127 

1 

,441 

P QUAD 2 

128 

1 

,830 

P0UAD2 

129 

1 

o830 

PQUAD2 

130 

1 

,441 

P UR J A2 

in 

1 

,531 

PTRIA2 

112 

1 

,532 

PTRIA2 

113 

1 

,396 

PTRI A2 

114 

1 

,544 

PTR1A2 

115 

1 

,590 

PTRSA2 

116 

1 

,591 

PTRKA2 

117 

1 

,557 

PTR1A2 

118 

1 

,519 

PIRJ A2 

119 

1 

,396 

PTRIA2 

120 

1 

,377 

RFORCE 

1 

0 



9 o o 10 


1 13* 34 1 





NAS IRAN 


EXECUTIVE 


CONTROL 


DECK. 


ECHO 


NASTRAN SYSTEM(76)= 1 $ SWEPT TURBOPROP OPTION 

SC NASA p SR5PR0P 

APR AERO 

SOL 9 

TIME 12 $ IBM 370/3031 

$ 

CHKPNT YES 
4 

44 4$ $$$$$$$$«$$$$$$$»$$$$$$&&&$$$ 

4 BELL/NASA NASTRAN AERO SCL 9 ALTERS TO PLOT MODE SHAPES AND 4 
4 INCLUDE RE AND PK METHODS OF FLUTTER ANALYSIS 4 

444444S44444444444444444444$444$4$«444444444444444444444$$4 44 $$$44444$ 


4 ALTERS TO PLOT MODE SHAPE So 

4 

ALTER 10 3o 103 

SDR2 CASECCoCSTMpMPT oDIToEQEXINoSILoo oBGPDTeLAMAo oPHIGoEST 9g /„ 
OPIHSGo o oPPHIG/C oN oRE IG $ 

CFKPNT PHIGoPPHIG $ 

ALTER 10 5 

PL T SET PCEB o EQEXlNcEC T /PLTSETZ oPLTPARZeGPSETSZo ELSETSZ/ 

SoNoNS IL2/SoNoJUMPZ=-l 4 
PRTMSG PLTSETZ // 4 
CCNO PZZo JUMP 2 4 

PLOT PL TP ARZ o GP SETSZ oELSE TSZ oCASECC o BGPD To EQEXINoS IL o oPPHI Goo/ 
PLOTZ/VoNoNSILZ/VoNoLUSET/VoN 0 JUMP2/VoNoPLIFLGZ=-l/ 

SoNoPFILEZ = 0 & 

PRTMSG PLOTZ // 4 
LAEEL PZZ 4 
4 

4$444444444$444444$4444$444$4444$4 $$$$$$$$ 4444444444 $$$m$ $4444444444 
EXIT $ INCLUDE ONLY IF TERMINATION AFTER MODE PLOTS IS DESIREOo $ 

4$ 44444444444 444$444$4$$44 $$ 444$4$44$ $$$$ 4$4444$444$44$$44$44 444444444 
8 

4 ALTERS FOR KE AND PK METHODS OF FLUTTER ANALYSIS- 

$ 

ALTER 152ol53 

FA1 KHEo BHHoMHHoQHHL oCASECCoFLl ST /F SAVE oKXHH oBXHHo MXHH/ 

So No FLOOP/ S oN o I START /So No NOCE AD $ 

EGUIV KXhHoPFilH/NOCEAQ/B XHH oCLAMA/NOCEAD/ 

KXHHoPHIHL/NCCEAD/BXHHoCLAMAL/NOCEAD/ 

CASECCoCASEYV/NOCEAD 4 


4.24 


N IS I R H EXECUTIVE CONTROL DECK ECHO 

I CONO VORoNOCEAD $ 

ALTER 156 
LABEL VDR $ 

ENCALTER 

I ^ 

§ CEND 

I 

I 

I 



I 

1 

I 

I 

K 

I 

I 

I 


4.25 






N A S T R A N EXECUTIVE CONTROL PECK ECHO 


^ECHO OF FIRST CARO IN CHECKPOINT OICflONART TO @E PUNCHED OUT FOR THIS PROBLEM 


RESTART NASA 0 SR5f>ROP o VIS/BZ* 47702„ 


4.26 


SR5 ADVANCED TURBOPROP FLUTTER ANALYSIS 
BELL/NASA NASTRAN RF 9 AERO o MODES 

1C SLADE So 6800 RPM„ o7Q TUNNEL MACH NQ<? 




CASE CONTROL DECK ECHO 

CARO 

COUNT 

1 i 

2 TITLE ■ SR 5 ADVANCES) TURBOPROP FLUTTER ANALYSIS 

3 SUBTITLE = BijiMNASA NASTRAN Rp 9 AEROo MODES 

4 ' LABEL a 10 BLADES® 6800 RPM 0 e 70 TUNNEL MACH NOo 

5 $ 

6 $PC = L 

7 METHOD = l 

8 VECTOR « ALL 

9 $ 

10 OUTPUT IRLOTI 

U SET 1 a ALL 

12 PLOTTER NASTPLT »MODEL D 0 0 

13 PAPER SIZE @ 0 5 8V 11.0 

14 MAXIMUM DEFORMATION 0»5 

15 FIND SCALE <? ORIGIN 1„SET 1 

16 PTITLE = SOL 9 AERO 

IT CONTOUR VOISPLAC 

18 PLGT MODAL DEFORMATION CONTOUR ,SET loORIDIN 1»PEN 2 

19 BEGIN BULK 


{• *»* USER INFORMATION MESSAGE 207 0 BULK DATA NUT SORTED „XSORT WILL RE-ORDER DECKo 


I 


k 

V' 


!> 

< 

ii 

& 


I 


THE BLADE SHAPE DEFINED BY THE GRID DATA IS THE DEFORMED SHAPE OBTAINED 
AT THE END OF THE DIFFERENTIAL STIFFNESS RUN. 

DUE TO THE REQUEST IN THE EXECUTIVE CONTROL DECK TQ EXIT AFTER PL0TTING MQDE 
SHAPES „ IT IS AT THE USER'S OPTION TO INCLUDE OR EXCLUDE THE AERODYNAMIC DATA FQR 
FLUTTER ANALYSIS, IN THIS EXAMPLE THE AEROpYNAMIC DATA HAVE BEEN LEFT OUT. 


I 

K 


I 


K 

I 

I 

I 




S 

0 R 1 

6 

0 B 

U 

l * 

DATA 

o l 

0 o 2 

a o 3 

9 o 

4 

o o 

5 

o o 

6 P o 7 

C CRD <2 ft 

A 


o0 


pQ 


oO 

— O 0277 

♦C2RI 

o 9996 

-o 02 77 

Qo0 






CQUAD2 

1 

1 

1 


2 


9 

8 

CQUAD2 

2 

2 

2 


3 


10 

9 

CQUAD2 

3 

3 

3 


4 


u 

10 

CUUA02 

4 

4 

4 


5 


12 

11 

C QUAD 2 

5 

5 

5 


6 


13 

12 

CQUA02 

6 

6 

6 


7 


14 

13 

C6UA02 

7 

7 

8 


9 


16 

15 

CQUA02 

3 

8 

9 


l!) 


17 

16 

CQUA02 

9 

9 

10 


U 


18 

17 

C QUAD 2 

10 

10 

n 


12 


19 

18 

C0UA02 

A I 

11 

12 


13 


20 

19 

CQUA02 

12 

12 

13 


14 


21 

20 

CQUA02 

13 

13 

15 


16 


23 

22 

GQUAD2 

14 

14 

16 


17 


24 

23 

C QUAD 2 

15 

15 

17 


18 


25 

24 

CQUA02 

16 

16 

lb 


19 


26 

25 

CQUA02 

17 

17 

19 


20 


27 

26 

CQUA02 

13 

18 

20 


21 


28 

27 

G QUAD 2 

19 

19 

22 


23 


30 

29 

C WUAG2 

20 

20 

23 


24 


31 

30 

C QUAD 2 

21 

21 

24 


25 


32 

31 

C0UAO2 

22 

22 

25 


26 


33 

32 

CQUA02 

23 

23 

26 


27 


34 

33 

CQUA02 

24 

24 

27 


28 


35 

34 

C QUAD 2 

25 

25 

29 


3p 


37 

36 

C QUAD 2 

26 

26 

30 


31 


38 

37 

CQUAD2 

21 

21 

31 


32 


39 

38 

CQUA02 

23 

23 

32 


33 


40 

39 

GQUAD2 

29 

29 

33 


34 . 


41 

40 

CQUAD2 

30 

30 

34 


35 


42 

41 

C6UA02 

31 

31 

36 


3? 


44 

43 

CQUA02 

32 

32 

37 


38 


45 

44 

CQUAD2 

33 

33 

38 


39 


46 

45 

CQUA02 

34 

34 

39 


40 


47 

46 

C QUAD 2 

35 

35 

40 


41 


48 

47 

CQUA02 

36 

36 

41 


42 


49 

48 

CQUAD2 

37 

37 

43 


44 


51 

50 

C0UA02 

38 

38 

44 


45 


52 

51 

CQUA02 

39 

39 

45 


46 


53 

52 

CQUAD2 

40 

40 

4?>> 


47 


54 

53 


E C M P 


9 


.0 


4,28 


o 10 

+C2K1 




I 

I 

ft* SORTED SULK DATA ECHO- 


CQUAD2 

41 

41 

47 

46 

55 

54 

| CQUAD2 

42 

42 

48 

49 

56 

55 

C QUAD 2 

43 

43 

50 

51 

58 

57 

CQUAD2 

44 

44 

51 

52 

59 

58 

r CQUAD2 

45 

45 

52 

53 

60 

59 

i CQUA02 

46 

46 

53 

54 

61 

60 

C QUAD 2 

47 

47 

54 

55 

62 

61 

r C QUAD 2 

48 

48 

55 

56 

63 

62 

| CQUAD2 

49 

49 

57 

58 

65 

64 

CQUA02 

50 

50 

58 

59 

66 

65 

CQUAD2 

51 

51 

59 

60 

67 

66 

f 7 C QUAD 2 

52 

52 

60 

61 

68 

67 

a- C0UAD2 

53 

53 

61 

62 

69 

68 

CQUAD2 

54 

54 

62 

63 

70 

69 

r CQUA02 

55 

55 

64 

65 

72 

71 

| A COUAD2 

56 

56 

65 

66 

73 

72 

CQUAD2 

57 

5 2 

66 

67 

74 

73 

r CQUAD2 

58 

58 

67 

68 

75 

74 

l£ C QUAD 2 

59 

59 

68 

69 

76 

75 

Wb ‘ CQUAD2 

60 

60 

69 

70 

77 

76 

CQUAD2 

61 

61 

71 

72 

79 

78 

§ C QUAD 2 

62 

62 

72 

73 

00 

79 

in CQUAD2 

63 

63 

73 

74 

81 

80 

CQUAD2 

64 

64 

74 

75 

82 

81 

r CQUAD2 

65 

65 

75 

76 

83 

82 

r CQUAD2 

66 

66 

76 

77 

84 

83 

C QUAD 2 

67 

67 

78 

79 

66 

85 

SP CQUAD2 

68 

68 

79 

80 

6 7 

86 

| CQUA02 

69 

69 

80 

81 

88 

87 

CQUAD2 

70 

70 

81 

02 

89 

88 

CQUAD2 

71 

71 

82 

83 

90 

89 

f CQUAD2 

72 

72 

83 

84 

91 

90 

1, CQUAD2 

73 

73 

85 

86 

93 

92 

CQUAD2 

74 

74 

86 

87 

94 

93 

r CQUAD2 

75 

75 

67 

88 

95 

94 

I CQUAD2 

76 

76 

88 

89 

96 

95 

89 CQUAD2 

77 

77 

89 

90 

97 

96 

CQUAD2 

78 

78 

90 

91 

98 

97 

ffl CQUAD2 

79 

79 

92 

93 

100 

99 

* CQUAD2 

80 

80 

93 

94 

101 

100 

CQUAD2 

81 

81 

94 

95 

102 

101 

If CQUA0 2 

82 

82 

95 

96 

103 

102 


I 

I 
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S 

0 R T 

£ 

0 B 

U 

L K 

D 

A T A 

o 1 

o a 2 

o a 

3 

© © 

4 

<9 o 

5 

© O 

6 

0 © 

CQUAD2 

83 

83 


96 


97 


104 


103 

C QUADS 

84 

84 


97 


98 


105 


lp4 

CQUAD2 

85 

85 


99 


100 


107 


106 

C QUAD 2 

86 

86 


100 


101 


108 


107 

CQUA02 

8? 

87 


101 


102 


109 


106 

CQUA02 

88 

88 


102 


103 


no 


109 

CQUAD2 

89 

89 


103 


104 


ui 


no 

CGUAD2 

90 

90 


104 


10 s 


112 


m 

C QUAD 2 

91 

91 


106 


10 7 


114 


113 

CQUA02 

92 

92 


107 


106 


115 


114 

CQUAD2 

53 

93 


108 


139 


116 


115 

C QUAD 2 

94 

94 


109 


no 


117 


116 

CQUA02 

95 

95 


uo 


in 


118 


117 

CQUAD2 

96 

96 


ill 


112 


119 


118 

C3UA02 

97 

97 


113 


114 


121 


120 

CQUAD2 

58 

98 


114 


f*5 


122 


121 

COUAD2 

99 

99 


115 


J , 

116 


123 


122 

CQUAD2 

ICO 

100 


116 


117 


124 


123 

C QUAD 2 

101 

101 


117 


118 


125 


124 

C3UAD2 

102 

102 


118 


119 


126 


125 

C 3UAD2 

103 

103 


120 


121 


128 


127 

C3UAD2 

104 

104 


121 


122 


129 


128 

CQUAQ2 

1 05 

105 


122 


123 


130 


129 

CQUAD2 

106 

106 


123 


124 


131 


130 

CQUAD2 

107 

107 


124 


125 


132 


131 

C QUAD 2 

108 

108 


125 


126 


133 


13 2 

CQUA02 

109 

109 


127 


128 


135 


134 

CQUAD2 

no 

110 


128 


129 


136 


135 

CQUA02 

121 

121 


131 


132 


144 


143 

CQUAD2 

122 

122 


132 


133 


145 


144 

CQUAD2 

123 

123 


138 


139 


147 


146 

CQUAD2 

124 

124 


139 


140 


148 


147 

C QUAD 2 

125 

125 


140 


141 


149 


148 

CQUAD2 

126 

126 


141 


142 


150 


149 

CQUA02 

127 

127 


146 


147 


152 


151 

CQUA02 

128 

128 


147 


148 


153 


152 

CQUAD2 

129 

129 


148 


149 


154 


153 

CQUAD2 

130 

130 


149 


150 


155 


154 

CFR IA2 

ill 

111 


129 


138 


136 



CTHIA2 

112 

112 


129 


137 


138 



C FRIA2 

113 

113 


129 


130 


137 



CTRI A2 

114 

114 


137 


130 


140 




I 

I 

4.30 
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o o 
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SORTED BULK DATA ECHO 


C TRIA2 

115 

115 138 

137 

139 




CTRIA2 

116 

116 139 

137 

140 




CTRIA2 

117 

117 140 

130 

141 




C TRIA2 

118 

118 141 

13Q 

142 




CTRIA2 

119 

119 130 

131 

*42 




C TRIA2 

120 

120 142 

131 

143 




C VJQ1 N 

1 

155 






CYJ01N 

2 

151 






E ICR 

1 

PEER 


5 


+EI 


+ E1 

MAK 







GRD SE T 


1 


l 




GRID 

n 



2o419817E*00 

lo244079E+0i 

4 GO 

1 

#GD 

i 

-2 o 031 032 E *00 






GRID 

42 



2o457173E+00 

1 *2381436+01 

*G0 

2 

*GD 

2 

-2*0582496+00 






GRID 

43 



2a 5516Q2E+00 

1 a 22521 6E + 01 

*GD 

3 

4GD 

3 

-2* 1145276+00 






GRID 

44 



2 o 6836286+00 

U202870E+01 

4 GO 

4 

4 GO 

4 

-2o 173213E+Q0 






GRID 

45 



2* 865225E+00 

lo!72653E+01 

4 GO 

5 

♦ GD 

5 

-2o 216011 E+QO 






GRID 

46 



2 o 990463E+00 

1 o 1 58475E +0 1 

4GD 

6 

4 GO 

6 

-2*254392E+00 






GRID 

47 



3 a I26074E+QQ 

1 o 14001 IE +01 

*GU 

7 

*GD 

7 

-2*266630E+00 
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8o703647E~05 

2o440I56E+00 

♦ GO 

148 

*GU 

148 -io 4^6011 E-05 





GRID 

♦ 149 


2o?00679E-0i 

2o440040E*Q© 

♦ GD 

149 

♦ GD 

149 -4o 00 93 9 7 £-02 





! 

GRID 

*150 


4o 700607E-01 

2o439936E*00 

♦ GO 

iso ; 

♦ GD 

150 -8o 013493E-Q2 





1 

1 

GRID 

*151 


-4o700000E-01 

2o060000E*00 

♦ GO 

151 i : 

*G0 

151 Bo 999997E-02 





j 

GR I D 

*152 


- 29 TOOOOQE-OI 

2 o 060000E*00 

♦ GD 

152 

*GU 

152 4o000000E-02 





j 

GRID 

*153 


o0 

2o06Q000E*00 

♦GD 

153 1 

£ GO 

153 0oC 





i * 

GRI0 

*154 


2 0 700000E-01 

2o060000E*00 

♦ GO 

154 |*' 

*GD 

154 -4 0 000 0 00 E -02 





j - 

GRID 

*155 


4o 700000E-01 

2oQbOQQOE*QG 

*GD 

155 r 


*GD L 55 -fl« 9999976-02 

MAT1 i 1 o 6 E7 o 35 <,0004141 

PARAM C TYPE ROT 

PARAM KGGIN 1 

PARAM KiNDEX 0 

PARAM WSEGS 10 

PQUAD2 1 1 <.012 

PQUAD2 2 1 oQ24 

PQUA02 3 1 o 032 

PQUA02 4 1 o036 

PQUAD2 5 1 o030 

P QUAD 2 6 1 p018 

PQUA02 7 1 o0 14 

PQUAD2 8 1 o028 

PQUAD2 9 l o037 

PQUAD2 10 1 c043 


4.38 


o o 


2 


o 1 

PQUAD2 
PQUAD2 
PQUA02 
PQUAD2 
PQUAD2 
PQUA02 
PQUA02 
PQUA02 
PQUA02 
' PQUAD2 
PUUA02 
PQUAD2 
. P QUAD 2 
PQUAD2 
T PQUAD2 

* PQUA02 
PQUAD2 

- PQUAD2 
PQUA02 
PQUAD2 
PQUAD2 

t PQUA02 
: PQUA02 

PQUAD2 

- PQUAD2 
PQUAD2 

' PQUAD2 
_ PQUAD2 
; PQUAD2 
" PQUAD2 
PQUAD2 
r PQUAD2 
1 PQUA02 
PQUAD2 

- PQUAD2 
l PUUAD2 

* P QUAD 2 
P QUAD 2 

| PQUA02 
PQUAD2 
PQUA02 
PQUAD2 


U 1 

12 l 

13 1 

A4 l 

15 1 

16 1 

17 1 

18 1 

19 1 

20 1 

21 1 

22 1 

23 1 

24 i 

25 1 

26 1 

27 1 

28 1 

29 1 

30 1 

31 1 

32 1 

33 1 

34 1 

35 1 

36 1 

37 1 

38 1 

39 1 

4C 1 

41 1 

42 1 

43 1 

44 1 

45 1 

46 l 

47 1 

48 1 

49 1 

50 1 

51 1 

52 1 


o o 4 
o036 
<>022 
o 0 16 

*032 

oQ48 

o051 

o042 

0023 
oQ 18 
oQ34 
o053 
o 058 
o 046 
o025 
,021 
o042 

o 06 1 
o066 

005 1 
O 027 

0024 
o049. 
o070 
«>073 
©057 

0 030 
<>028 
« 054 
o078 

o082 

o 065 
o 035 

003 1 

006 1 

o088 

o093 
o075 
©039 
o 038 

o 06 8 

©098 
o 103 


o 1 

o a 2 

O o 

3 Q D 1 

PyUAD2 

53 

1 

0 Q83 

PQUA02 

54 

1 

.046 

PQUAQ2 

55 

1 

o041 

PQUAD2 

56 

1 

o 076 

PQUA02 

57 

1 

ollO 

PQUA02 

58 

1 

.118 

PQUAU2 

59 

1 

.091 

PQUAD2 

60 

1 

.047 

PQUA02 

61 

1 

.043 

PQUAD2 

62 

1 

.083 

P0UAD2 

63 

1 

o 1 20 

PQUA02 

64 

1 

.129 

PUUA02 

65 

l 

ol00 

PyUAD2 

66 

l 

o 044 

PWUAU2 

67 

1 

oQ45 

PyUAD2 

68 

1 

o Q90 

PUUAU2 

69 

1 

.135 

PWUA02 

70 

1 

.138 

pyUAD2 

71 

1 

.100 

PUUA02 

72 

1 

.048 

P QUAD 2 

73 

1 

.053 

PQUA02 

74 

1 

.106 

P QUAD 2 

75 

1 

.152 

PQUAD2 

76 

1 

.148 

pyUAD2 

77 

1 

.099 

PwUAD2 

78 

1 

o 044 

PQUAD2 

79 

1 

o063 

PWUAD2 

80 

1 

o 123 

PQUAD2 

81 

1 

o 171 

PUUA02 

62 

1 

.157 

PQUA02 

83 

1 

o 099 

P QUAD 2 

64 

A 

.046 

PQUAD2 

85 

1 

o071 

PQUAD2 

£6 

1 

.141 

PQUAD2 

87 

1 

.206 

PUUAD2 

68 

1 

.177 

PQUAD2 

85 

1 

.112 

PQUAD2 

90 

1 

.048 

PQUAD2 

91 

1 

.084 

PQUAD2 

92 

1 

.172 

PQUAD2 

93 

1 

.232 

PQUAD2 

54 

1 

.198 


SORTED BULK DATA (ECHO 


• 1 

a o 2 

o o 

3 

0 o 4 

o o 5 

o o 

PQUAD2 

95 

1 


= 135 



PGUA02 

96 

1 


o062 



PQUA02 

97 

1 


o 1 19 



PQUA02 

90 

1 


,206 



PUUA02 

99 

1 


= 266 



PQUAD2 

100 

l 


o 230 



PGUAD2 

101 

1 


o 152 



PQUAD2 

102 

1 


o07 1 



PQUAD2 

103 

1 


o 16 1 



PUUAD2 

104 

1 


o 237 



PaUA02 

105 

1 


o347 



PUUAD2 

106 

1 


o 3 19 



PGUAD2 

1C7 

1 


O 167 



PUUAD2 

1CB 

1 


o0?5 



PQUAD2 

1C9 

1 


= 222 



PQUAD2 

110 

1 


,373 



PQUAD2 

121 

1 


o242 



PQUA02 

122 

1 


o089 



PGUAD2 

123 

1 


,441 



PtrfUAD2 

124 

1 


,830 



PGUA02 

125 

1 


o 830 



PQUA02 

126 

1 


= 441 



PGUA02 

127 

1 


o44 1 



PUUA02 

128 

1 


,830 



PQUAD2 

129 

1 


o830 



PGUAD2 

130 

1 


o 441 



PTR1A2 

111 

1 


a 53 1 



PTRIA2 

112 

1 


= 532 



PTR I A2 

113 

1 


o396 



PTRIA2 

114 

1 


= 544 



PTRIA2 

115 

1 


= 590 



PTRIA2 

116 

1 


= 591 



PTRIA2 

117 

1 


= 557 



PTRIA2 

118 

1 


= 519 



PTKI A2 

119 

1 


o 396 



PTRIA2 

120 

1 


= 377 



RFOKCE 

1 

0 



113=34 

H 

O 

O 

SPC i 

1 

4 


I 

57 


SPC1 

1 

6 


7 

91 

98 

SPC1 

1 

123456 

151 

THRU 

155 

ENOOATA 








134 


145 



U S IR U EXECUTE V E CONTROL DECK ECHO 

NASTRAN SYSTEM (76)= 1 $ SWEPT TURBOPROP OPTION 

i D NASA d SR 5PROP 

APP AERO 

SOL 9 

% 


RESTART NASA 

b SR 5PR0P 

1 bOO/CO/COb 

00000 B 





1 0 

XVPS 

0 

FLAGS ® Ob 

REEL = 

1 9 

FILE 

s 

6 

2b 

REENTER 

AT 

OMAP SEQUENCE 

NUMBER 

7 




3b 

GPL 

V 

F LAGS = 0 o 

REEL = 

1 0 

FILE 

s 

7 

4b 

EQEXIN 

9 

flags ° Ob 

REEL = 

1 B 

FILE 

s 

8 

5b 

GPDT 

0 

FLAGS ° Ob 

REEL = 

1 B 

FILE 

SE 

9 

6 b 

CSTM 

9 

F LAG S c Ob 

REEL = 

1 B 

FILE 

= 

10 

7 o 

BGPOT 

0 

FLAGS = Ob 

REEL = 

1 V 

file 

a 

11 

8b 

S1L 

9 

FLAGS - Oo 

REEL - 

1 0 

FILE 

3 

12 

9 B 

XVPS 

9 

F LAGS a Ob 

REEL = 

1 0 

F ILE 

= 

13 

lOo 

REENTER 

AT 

DMAP SEQUENCE 

NUMBER 

10 




11. 

ECT 

9 

FLAGS = Ob 

REEL = 

1 B 

FILE 

= 

14 

12* 

XVPS 

9 

F LAG S 3 Ob 

REEL = 

1 B 

FILE 

= 

15 

13* 

REENTER 

AT 

OMAP SEQUENCE 

NUMBER 

12 




14o 

XVPS 

9 

F LAG S = Ob 

REEL = 

l B 

FILE 

= 

16 

15b 

GPTT 

9 

FLAGS = Ob 

REEL = 

Ob 

FILE 

S 

0 

16 o 

REENTER 

AT 

OMAP SEQUENCE 

NUMBER 

17 




17 B 

EST 

0 

FLAGS = Ob 

REEL = 

1 B 

FILE 


17 

18b 

GPECT 

9 

FLAGS = Ob 

REEL = 

l 0 

FILE 

= 

18 

19 o 

XVPS 

0 

F LAGS a 0 b 

REEL = 

1 0 

FILE 


19 

20b 

GEI 

9 

FLAGS = Cb 

REEL = 

Ob 

FILE 

s 

0 

21b 

OGPST 

9 

FLAGS = Ob 

REEL = 

Ob 

FILE 

= 

0 

22b 

REENTER 

AT 

OMAP SEQUENCE 

NUMBER 

29 




23b 

MELM 

9 

FLAGS * Ob 

REEL = 

1 6 

FILE 


20 

24 o 

MDICT 

9 

FLAGS = Ob 

REEL = 

1 0 

FILE 

= 

21 

"5b 

XVPS 

9 

FLAGS ® Ob 

REEL = 

1 B 

FILE 


22 

«, 8 B 

KELM 

0 

FLAGS = Ob 

REEL = 

0 0 

file 

= 

0 

27b 

KD1CT 

0 

FLAGS = 0„ 

REEL = 

0 B 

FILE 

= 

0 

28b 

REENTER 

AT 

OMAP SEQUENCE 

NUMBER 

36 




29b 

MGG 

0 

FLAGS = Ob 

REEL = 

1 B 

FILE 

= 

23 

30 b 

XVPS 

9 

FLAGS = Ob 

REEL = 

1 0 

FILE 

s 

24 

31b 

REENTER 

AT 

OMAP SEQUENCE 

NUMBER 

42 




32b 

KTOTAL 

0 

FLAG 5 = 4 p 

REEL = 

1 0 

file 

a 

25 

33b 

KGGX 

9 

FLAGS = 4 b 

REEL = 

1 0 

FILE 

= 

25 

34 b 

KGG 

9 

F LAGS — 4 o 

REEL = 

1 0 

FILE 

= 

25 

35* 

XVPS 

9 

FLAGS — Ob 

REEL = 

1 B 

FILE 

= 

26 

36 b 

REENTER 

AT 

OMAP SEQUENCE 

NUMBER 

54 





4 . 4 ? 


—I <*> 4ft 


DECK 




N A S T R A N EXECUTIVE 


CUNT 


ROC 


37 o 

CYCO 

9 

FLAGS 3 Os 

REEL = 

A s 

FILE 


27 

38 o 

XVPS 

0 

FLAG S 3 0 s 

REEL 3 

1 s 

FILE 

P 


39* 

REENTER 

AT 

DMAP SEQUENCE 

NUMBER 

63 



AO* 

KNN 

9 

FLAGS p 4s 

REEL = 

1 s 

FILE 

P 

25 

41® 

MGG 

9 

FLAG S 3 4 o 

REEL = 

I s 

FILE 

s 

23 

42 o 

MNN 

9 

F LAGS 3 4 s 

REEL 3 

i s 

FILE 

c 

23 

43* 

XVPS 

0 

FLAGS 3 Os 

REEL 3 

I s 

FILE 

o 

29 

44 s 

REENTER 

AT 

DMAP SEQUENCE 

NUMBER 

71 




45 s 

XVPS 

0 

FLAG S 3 0 s 

REEL 3 

1 s 

FILE 

c 

30 

46 s 

KFF 

9 

FLAGS = 0® 

REEL = 

0 s 

FILE 

* 

9 

47s 

MFF 

0 

FLAGS = Os 

REEL = 

Os 

FILE 

s 

0 

48s 

REENTER 

AT 

DMAP SEQUENCE 

NUMBER 

74 




49* 

KFF 

9 

FLAGS 3 Os 

REEL 3 

I s 

FILE 

3 

31 

50s 

KFS 

V 

FLAGS 3 0® 

REEL 3 

Is 

FILE 

S 

32 

5 1 s 

MFF 

9 

FLAGS 3 Os 

REEL 3 

1 S 

FILE 


33 

52s 

XVPS 

V 

FLAGS * Os 

REEL = 

1 s 

FILE 

= 

34 

53 s 

REENTER 

AT 

DMAP SEQUENCE 

NUMBER 

77 




54s 

KFF 

0 

FLAG S 3 4 s 

REEL 3 

1 0 

FILE 

- 

31 

55s 

KAA 

0 

FLAGS = 4® 

REEL 3 

1 s 

FILE 


31 

56s 

MAA 

9 

FLAGS = 4 o 

REEL 3 

1 s 

FILE 

= 

33 

57s 

MFF 

0 

FLAGS 3 4s 

REEL 3 

1 0 

FILE 

J= 

33 

58s 

XVPS 

0 

FLAGS 3 Os 

REEL = 

1 s 

FILE 

= 

35 

59s 

REENTER 

AT 

DMAP SEQUENCE 

NUMBER 

90 




60s 

KKK 

0 

FLAGS = Os 

REEL = 

1 ® 

FILE 

- 

36 

6 1 s 

MKK 

9 

FLAGS 3 Os 

REEL = 

1 s 

FILE 

P 

37 

62s 

XVPS 

9 

FLAGS 3 Os 

REEL 3 

1 s 

FILE 

= 

38 

63 * 

REENTER 

AT 

DMAP SEQUENCE 

NUMBER 

94 




64* 

LAMK 

9 

flag s = o s 

REEL 3 

1 s 

FILE 


39 

65* 

PHIK 

D 

FLAGS - Os 

REEL 3 

1 s 

FILE 

= 

40 

66s 

OEIGS 

0 

FLAGS = Os 

REEL 3 

1 s 

FILE 

s 

41 

67 s 

XVPS 

9 

FLAGS = Os 

REEL 3 

1 s 

FILE 

- 

42 

68 s 

MI 

9 

FLAGS 3 Os 

REEL 3 

0 s 

FILE 

a 

0 

69s 

KEEN TER 

AT 

DMAP SEQUENCE 

NUMBER 

101 




70s 

LAMA 

9 

FLAGS 3 Os 

REEL 3 

1 S 

FILE 

= 

43 

71s 

PHIA 

9 

FLAGS = Os 

REEL 3 

1 ® 

FILE 

= 

44 

72s 

XVPS 

9 

F LA G S 3 Os 

REEL 3 

1 s 

FILE 

3 

45 


ENC OF CFECKPOINT DICTIONARY 
IME 35 $ IBM 370/3031 

$$ $$$$$$>$$ $$>$&$$$$>$&$> $>$$$$$$> $$i $$$$$$$ 

$ BELL/NASA NASTRAN AERG SOL 9 ALTERS TO PLOT MODE SHAPES AND 


ECHO 


j 


i 

> 


1 

I 

] 

\ > 


r 


t 


i 

» 


i 

i 


) 

\ 


i 




N A S T R A N 


EXECUTIVE 


CONTROL 


0 E C K 


ECHO 


$ INCLUDE K£ AND PK ME THCDS OF FLUTTER ANALYSIS $ 

i $ $$$$$*$$ 4 . 4 $ $$$$$$$$ $$$$$$$&$$$$$$$$$$&$$$$>&& $$$$$&$$$$$ 

$ 

4 ALTERS TO PLOT MODE SHAPESo 

% 

ALTER 103ol03 

SUR2 CASECCpCSTMoMPTsDITffEOEXINoSI L o p o BGjPO ToLAMA 0 ppHIGpEST 0 0/ p » 
OPHIGpppPPHIG/CpNpREIG 4 
CHKPNT PH IG0PPHIG $ 

ALTER 105 

PLTSET PC08 ® EGEXINpECT/PLTSETZpPLTPAR^pGPSETSZpELSETSZ/ 

SpNoNSILZ/SoN, JUMPZ =-1 4 
PRTMSG PLTSETZ U $ 

COND PZZpJUMPZ $ 

PLOT PL TPARZj GPSETSZ 0ELSE TSZpCASECG oBGPOT pEQEX I N 0S IL 0 0 PPHI Go 0/ 
PLOTZ/V0N0NSILZ/V0N0LUSET i v S oN 9 JUMPZ/V 0N oPLTFLGZ^-l/ 

So N pPF ILEZ = 0 4 
PRTMSG PLOT Z // $ 

LABEL PZZ 4 
$ 

4 ALTERS FOR K£ AND PK METHODS OF FLUTTER ANALYSIS. 

4 

ALTER 1520153 

FA1 K HHo BHH0MHH0UHHL 0CASECC 0FLI ST/FSAVEPKXHH0BXHH0MXHH/ 

S0N0 FLOOP/S0N9T START ZSpNo NOCEAO 4 
EQUIV KXHH0PHI h/NOCEAD/BXHH pCLAMA/NGCE AD/ 

KXHH0PHIHL/NOCEAD/BXHH0CLAMAL/NQCEAO/ 

CASECC oCASE YY/NOCE AD $ 

COND VDRoNQCEAD 4 

ALTER 156 
LABEL VDR 4 
ENCALTER 
4 

CEND 
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1r5 advanced turboprop flutter analysis 

BELL/NASA NASTRAN RF 9 AERCo MODES 
So BLADES 0 6800 RPM 0 c 70 TUNNEL MACH N0 o 


UrtltAWAL PAGE US 

OF POOR QUALITY 



I 


I 


CASE CONTROL DECK ECHO 

CARD 

COUNT 

1 $ 

2 TITLE » SR 5 ADVANCED TURBOPROP FLUTTER ANALYSIS 

3 SUBTITLE « BELL/NASA NASTRAN RF 9 AEROo F LUTT£R 

4 LABEL ° 10 BLADES 6000 RPM 0 o T0 TUNNEL MACH NO, 

5 $ 

6 SPC O 1 

7 METHOD a 1 

8 F METHOD = 1 

9 $ 

10 OUTPUT 4 XYOUTJ 

n i 

12 PLOTTER NASTPLT D o 0 

13 XPAPER a 8 0 5 

14 YPAPER = IloO 

15 YAXIS * YES 

16 ^INTERCEPT = 7046«0 $ OPERATING VELOCITY 

17 XTAX5S = YES 

18 XBAXIS = YES 

19 CURVELINE SVMBCL = 6 

20 XDIVI SIGNS = 10 ] 

21 YT05 VI S1CNS = 10 

22 YBDS VISIONS = 10 

23 YTGRID LINES = YES j 

2k YBGRID LINES = Vt S } 

25 XTGRIO LINES s YES 

26 XbGRlU LINES = YES 

27 XtllLE = VELOCITY VSBAR <. IN/SEC o 0 0 » REF VS BAR = 7046 IN/SEC^ 

28 VI TITLE = DAMPING G 

29 VBTITLE = FREQUENCY l- 0 HZ 

30 TC UR VE — K = e. i.t>o2 o « — 006 go9 ol «£ 1 1 o5oSI G“0 oO * 

31 XYPLQT oXVPRI NT VG/ UG 0 FJ 0 24 G 0 FK 3(GoF ) 0 4lG 0 FJo SIGoFK 6 {GoF.| 

32 TC UR VE = Ka o i ( , o 2oo30o6oo90lo2i)lo59SIG = 36oO | 

33 XVPLOT oXYPRINT VG/ 7IG 0 F), SflGoFl* 95 Go F i 9 10 1 G 0 F ) 0 1 1 ( G„ F ) 0 124 G, F : 

34 TCURVE = Ka » 1 0 o 2 0 3 0 o 6 0 o 9 o 1 o 2 0 1 o 5 0 S I G = 7 2 oO i 

35 XYPLUT 0 XYPRI NT VG/13 4 G 0 ?! *14 5 G oF I 0 155G* F J 0 16 5 G, FJ 0 171 G v F J 0 1£4 G B F ; 

36 TCURVE = K=ol »«2oo3 s ,6m9o1«2i)1<.50SI6=IO8oO i 

37 XVPLOT 0 XYPRI NT VG/19 4G 0 H c20 4 G oFl o21 4 G 0 F i ,22 1 &« FI 0 2 35 Go F ) » 245 Go F ! 

38 TCURVE = K=ol0o20 O 30o6oo90io2olo5oS!G=144oO 

39 XVPLOT oXYPRINT VG /25 1 G 0 FI *26 4 G »F i *27 4G« F I 0 28 1 G 0 F 1 0 295 &* F 1 0 305 G„ F 

40 TCURVE = K = o 1 0 o 2 0 o 3 0 06 0 o 9 0 1 o 2 0 1 o 5 9 S I G = 1 80 o 0 

41 XYPLOT 0 XYPRI NT VG/31 4GoF) 932 5G e Ft 033(G«Fi o34(GoFl o35« G ? F)o 36( G 0 F 

42 TC UR VE a K= ol 9 o 2 oo3 9o60q9ol-j29 l o 5 0 $ IGa— 144 0 O I 
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5 ADVANCED TURBOPROP FLUTTER ANALVSI S 
ELL/NASA NAS TRAN RF 9 AERO o NCpES 

'O BLADES* 6800 K.'M* o 70 TUNNEL MACH NOo 


ORIQtMAL PAGE IS 
OF POOR QUALITY 


CARD 

COUNT 

43 

44 

45 

46 

47 

48 

49 

50 


CASE 


CCNTAOL DECK ECHO 


AVPLQT oXVPRJNT VG/37flG<,Fi 0 38 UGoFJ «39«GpFi o*0(GoF) o4A<! Go Pi o 424 
TCURVE = K= oi o o 2 0 G> 3 Do69o9oio2elo55S 3G=-I08 oO 

XVPLOT oKVPRI NT VG/43(5G 0 IFJ 0 44 UGoFJ »45« Go f D o46 f Go F» o 47H Go PJo <484 
TCURVE =K=oi 0o2co3oo6pp 9 0 1 0 2 o i o5 oS I G= -r7 2a 0 

KVPLOToKVPRXNT VG/498G 0 FI 0 ?0(G 9 F ) •5L{G*PI «52( Go F3 v 531 G»E )« 544 
T C UR VE =K ss oioo2oo3oo6(?o9olo2oAo5()SJ **36 9 0 

XVPLOT pXVPRINT VG/55«GoF!) P 56 SG oF A 0 57fl GoF i 0 58 * GoFA 0 5p« Go F i * 60« 
BEGIN BULK 


Go F I 
GoF f 
GoF 1 
GoF i 


<r -3 


4.46 




I 
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ORTO.L PAQE 18 
OF POOR QUALITY 


INPUT BULK. DATA DECK £ C H 0 


9 


10 


aerodynamic data pgr flutter analyses 


AERO 

0 

1 = 0 

2 = 905 

9 = 763£= 

S 




f FLFACT 

i 

1 = 0 







t FLFACT 

2 

0=0 

36=0 

72 = 0 

108 = 0 

144 = 0 

180 = 0 

-144=0 +FL21 

*FL2l 

-iOboO 

-72 = 

-36 o 






B" FLFACT 

3 

= 10 

= 20 

= 30 

= 6 

= 9 

1=2 

1 = 5 

fl FLUTTER 

1 

KE 

1 

2 

3 

L 

6 


MKAER02 

0= 0 

= 001 

0=0 

= 3 

0 = 0 

= 6 

0 = 0 

= 9 

MKAER02 

36 0 

= 001 

36 = 

= 3 

36 = 

= 6 

36 = 

=9 

j MKAER02 

72 o 

= 001 

72 = 

= 3 

72 = 

o 6 

72= 

= 9 

8s MKAERQ2 

IC8p 

= 001 

108 = 

= 3 

108 = 

= 6 

108 = 

= 9 

MKAERG2 

144 = 

= 001 

144 = 

= 3 

144 = 

= 6 

144 = 

P 9 

f MKAERC2 

100= 

= 001 

180 = 

= 3 

180 = 

= 6 

180 = 

= 9 

L MKAER02 

-144 = 

,001 

-144 = 

= 3 

-144 = 

o 6 

— 144 o 

= 9 

MKAERU2 

-lObo 

= 001 
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❖ GD 

51 

-2=0473316-01 


GRAD 

❖ 52 
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❖ GD 

52 

-3=1295406-01 


GRAD 
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❖ GD 

53 

-5=0075646-01 


GR I D 
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1=1036806*00 
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54 

-5= 8475206-01 
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❖ 55 
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1=0120886*01 
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1 = 1230936*01 
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1 = 0502406*01 
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9=625038E*00 
1=0994766*01 
1=08566 OE *01 
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9o746097E*00 
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59 
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♦ GD 

61 
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♦ GO 

62 
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62 
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63 
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♦64 
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64 
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69 

-2o 687562 E-01 
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-X©X35652E*00 

8 © X 51909£4 , 00 
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GRID 
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GRID 
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GRID 
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103 

4 0 400A90E-02 

GRID 

♦ 104 


♦ GO 

104 

■1 o 597679E-01 

GRID 

♦105 


♦ GO 

105 

■2o685744E-01 

GRID 

♦ 1C6 


♦ GO 

106 

8o 496489E-01 

GRID 

♦ 107 


♦ GO 

107 

7o 381414E-01 

GRID 

♦ 1C8 


♦ GO 

108 

4o 893306E-0 1 

GRID 
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♦ GO 

109 

1 o 591752E-01 

GRID 
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♦ GO 
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■6o 520444E-02 
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♦ GO 
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•1 o 757268E-01 

GRID 

♦112 


♦ GO 
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2o 750542E-0 1 

GRID 
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♦ GO 

113 

6o 919722E-01 

GRID 

♦ 114 


♦ GO 

114 

5o 999017E-01 

GRID 

♦ 115 


♦ GO 

115 

3©958974E-01 

GRID 
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♦ GO 

116 

lo095533E-01 

GRID 
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♦ GO 

117 

•8 D 331943E-02 

GRID 
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lo008383E*00 
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->5o963364£-01 

Ao954939E-01 

6 o 204426E-0i 

lo050502E«-00 

-2 <>496761 £<>00 

-2ol94051E*00 

-lo543487E*00 

-3o069078E-0l 
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6 o 5 11260E-01 

lo 101048E*00 

-2 ©2 73488E<>00 
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o 
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98 
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♦ GO 

99 

5©9O23UE*0O 

♦ GO 

100 
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♦GO 

101 
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♦ CD 
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4©999296E*G0 

♦GO 

103 
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♦ GO 

104 

4o 750490E <>00 

♦ GO 

105 

4o997375E*00 

♦ GD 

106 

4o997557£400 

♦ GD 

107 

4©903041E*00 

♦ GD 

108 

4©75189QE*00 

♦ GD 

109 

4 <>620979£ 400 

♦ GO 

no 

4o550488E<-00 

♦ GD 

in 

4©450042E*00 

♦ GD 
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4©304033E*00 

♦ GD 
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4o303668E+0C 
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♦ GD 
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118 


4.57 


b* WSHjjj&jat 


OF POOH OU “UTY 


o 1 

o o 

2 o o 3 a 

*GU 
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*GD 
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GRID 

*121 


*GD 

121 

4o447656E-01 

GRID 

*122 


*GD 

122 

2o925791E-01 

GRID 

*123 


*GD 

123 

4 o970558E-02 

GRID 

*124 


*GQ 

124 

•1 oll8595E-01 

GRID 

*125 


*G0 

125 -2o023238E-0l 

GRID 

*126 


*GD 

126 

■2a 828727E-01 

GRID 

*127 


*GD 
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3o625186E-01 

GRID 

*128 
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GRID 
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*GD 

130 

1«195536E-01 

GRI D 
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♦ GO 
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lo414182E-01 

GRID 

*132 


*GD 
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■2o2i 7990E-01 

GRID 

*133 


*GD 

133 

•2 o 995755 E— 0 1 

GRID 

*134 


♦ GO 

134 

2o 804042 E— 01 

GRID 

*135 


♦ GD 

135 

2c2 14295 E-01 

GRID 

*136 


*GD 
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1 o 307743E-01 

GRID 

*137 
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9o 034532 E-02 
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138 
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3 o 62391 1 E +00 

*G 0 
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i 
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4 
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s 
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126 1 
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i 
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*GD 
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*GD 
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*GD 
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GRID 
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f 

*GD 
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GRID 
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*GD 
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GRID 

*145 


L 
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r 
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147 
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